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Correlation of Polarized Light Phenomena With the 
Orientation of Some Metal Crystals 


C. J. Newton and H. C. Vacher 
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Introduction 

In continuation of an investigation |/ ?\*? of the 
utility of plane polarized light in metallography 
a study has been made of tin and aluminum crystals, 
and additional data have been obtained from monel, 
a 70 nickel-30 copper alloy. The previous work 
at the National Bureau of Standards by D. H 
Woodard [/] showed a nonuniformity of intensity 
of polarized light reflected from plastically deformed 
metal grains, suggesting an observable correlation 
between such intensity and crystallographic orienta- 
tion. The specific problem of such correlation was 
pursued in a study by H. C. Vacher [2], who showed 
a correlation between extinction positions and the 
projection azimuth of a cubic axis 

In the present study the relative maximum and 
minimum intensities of the reflected light were 
measured photometrically and correlated with crystal 
orientation. Woodrow, Mott, and Haines [3] used 
a photometric method but gave no data correlating 
intensities with orientation. Supplementary experi- 
ments were made to provide information concerning 
the cause of the observed optical anisotropy, which 
has been the subject of some uncertainty [4] 


2. Theoretical Aspects 


The theory of the reflection of plane polarized 
light from metallic specimens is one of considerable 
complexity [5]. This discussion will be limited to a 
qualitative description of the special case of a plane 
polarized beam incident normally upon the reflecting 
surface [6], which is very nearly the case when the 
metallographic microscope 1s used with low-powered 
objectives 

For isotropic metals (cubic structure), normally 
incident plane polarized light is reflected inherently 
with neither rotation of the plane of polarization nor 
introduction of elliptical polarization; therefore, 
brackets indicate 
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examination by rotation of the stage of the micro 
scope with crossed nicols exhibits continuous extine 
tion. A like situation exists when the optic axis of 
the ervstal is perpendicular to the surface of a ervstal 
belonging to the uniaxial classificat ion, W hic h includes 
the tetrogonal, hexagonal, and trigonal systems If 
the optic axis and the surface normal do not coincide 
it is necessary to consider the orientation of the plane 
they determine. The intersection of this plane with 
the surface determines a line ealled the principal 
direction, and the incident light an re 
solved Into components parallel to and perpendi ular 
to this direction Reflectivities and phase changes 
upon reflection differ for these two components The 
reflected light Is, In general, elliptically polarized, and 
the field observed through the crossed analyzer of the 
microscope bright There however, four 
special positions in a complete revolution of the stage 
where the plane of vibration of the normally incident 
beam is parallel to or perpendicular to the principal 
direction in the specimen surface At these special 
positions of the stage, the reflected light remains plane 
polarized in the original plane and is extinguished by 
the analy zer 

The theory of reflection from orthorhombic, mono 
clinic, and triclinic materials is of great complexity 
and will not here Fortunately these 
cases of lower symmetry rarely arise in metallurgical 
problems 

The confined to the 
effect on polarized light of the intrinsic optical char- 
acteristics a single flat surface In view of the 
possibility, however, that the optical anisotropy of 
metaliographic specimens may be caused by surface 
contours, the process of multiple reflection of plane 
polarized light must be considered To this end, 
experiments were conducted as will be deseribed 
later with models made with one, two, or three per- 
pendicular front-surfaced aluminized mirrors 

The polarization figure [7] is another phenomenon 
that can sometimes be observed with the polarizing 
microscope and, if distinet, can be used in crystal 
orientation problems. It resembles the interference 
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FicurReE 1 Polarization figures formed by reflection from mone 


figure [8] observed with convergent plane polarized 
light passing through transparent anisotropic mate- 
rials. Instead of arising from path differences 
caused by birefringence, the polarization figure, 
according to Cameron and Green |9], results from the 
rotation and ellipticity introduced in reflected polar- 
ized light when the incident beam is not normal to the 
reflecting surface. The procedure for observing the 
figures is given in section 3.3. When this procedure 
is followed, a bright field containing two dark lines, 
called the isogyres, will be observed. These arise 
from reflections at those special angles from the 
surface where, because of particular relations of the 
planes of the polarizer and analy zer, the angle and 
plane of incidence of the convergent beam and the 
optical characteristics of the surface itself, there is a 
cancelling out of the various factors that disturb the 
plane polarization in the reflected beam From those 
special angles alone, therefore, is the reflected light 
extinguished by the analyzer. As illustrated in 
figure 1, the isogvres generally appear as two branches 
of a hyperbola in the field, but a cruciform figure is 
formed at four special positions These are the 
same positions that give extinctions when the inci- 
dent light is quasi-normal as is the case withfa low- 
powered objective. The details, even if only quali- 
tative, of the cause of this phenomenon are rather 
involved. For a reasonably elementary explanation, 
the reader is referred to an article by E. N. Cameron 
and L. H. Green [9], who have studied these figures 
as they relate to the optical properties of ore minerals. 


3. Materials and Procedures 
3.1. Materials 


The greatest experimental difficulty encountered 
in this study was the securing of suitable metallo- 
graphic specimens with large grains and with sur- 
faces as free as possible from imperfections such as 
fine scratches, irregular pits or nonuniform films, 
which would mask the pertinent polarized-light 
phenomena. In order to avoid errors arising from 
differences in surface preparation, polycrystalline 
specimens were used. Many specimens were exam- 
ined, but the data presented in this report were 
obtained from four: One each of tin and aluminum 
and two of monel. The purity of the tin was 99.96 
percent and that of the aluminum, 99.99 percent. 
The monel was commercial grade, which had been 
heat treated to grow large crystals [2]. The prepara- 
tion of each specimen is described in section 4. 


3.2. Determination of Crystal Orientation 


A major preliminary activity was the determina- 
tion of the crystallographic orientation of each grain 
used in each specimen. This orientation was defined 
with respect to the normal to the surface of the speci- 
men and a reference mark on the surface defining 
zero azimuth. 

In the case of monel, these determinations were 
made by means of angular measurements of the 
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es of twinning planes in two mutually perpendic- 
surfaces |/0]. Monel is a face-centered cubic 
tal, and the twinning planes are (111) planes. 
th this information, and with the aid of stereo- 
hic plots and standard projections, it was possi- 
to establish the orientations of groups of twinned 
tals near the edges of the specimens 
pecimens of aluminum and tin were obtained 
grains large enough for the determination of 
itations by X-ray diffraction, using the back- 
ction Laue method with a Greninger net (77, 12]. 
s method was very simple to employ in the case 
he cubie aluminum, but considerable diffieulty 
encountered with the tetragonal tin. In prae- 
one pattern alone was not sufficient with the 
ter to allow the assignment of an unquestioned 
orientation to a grain. It was necessary to make 
ne or more check patterns with the X-ray-beam 
dent at an angle to the specimen surface so as to 
parallel to some low-index crystallographic direc- 
on. according to the tentative interpretation of the 
st exposure. If the resulting pattern agreed with 
the proper standard pattern made previously from 
another grain whose orientation was known with 
assurance, the orientation of the new grain was 
considered determined. 


3.3. Measurement of Reflected Polarized Light 


The microscope used in this study uses a Foster 
prism [{/3] as the polarizing-analyzing element. The 
objectives were nominally strain-free. This prism 
acts as permanently crossed nicols; hence, no effect 
of slightly uncrossing the analyzer was investigated 
The prepared specimen was placed on the rotating 
stage of this polarizing microscope and set in the 
zero azimuth position with reference to an edge, 
grain boundary, or scratch on the surface. A 5.6x 
objective lens and a 5x ocular were used. By means 
of stage adjusting screws, a particular grain was 
centered under the cross lines on the ground g 
the microscope camera. After satisfactory centering 
and magnification, the bellows was adjusted until 
the grain under study would cover the photocell 
aperture. The ground glass was removed and a 
slide bearing the photomultiplier tube, centered and 
limited by a %4o-in.-diameter aperture, was put in its 
The light flux passing into the sensitive tube 
caused a proportional deflection of the needle on the 
dial of a commercial electronic photometer. As the 
specimen was rotated, the intensity of the reflected 
light would rise and fall in an approximately sinu- 
soidal pattern. The usual procedure was to make 
three complete revolutions of the stage while record- 
ings were taken of the angular settings of the stage 
and the photometer current proportional to the light 
intensity at each of the four maxima and each of the 
four minima. 

[t was observed for any one surface that the mini- 
mum intensities from all grains were nearly the same, 
arising from the scattered light from surface imper- 
lections and from the microscope elements. All of 
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the maximum intensities from a grain were averaged 


together, and the difference taken between that 
average and the average minimum intensity read- 
ings from the grain. This intensity difference, AJ, 
is the quantity used in the reported results. Other 
intensity functions, such as the ratios of various in- 
tensities, were examined but were found to be less 
satisfactory. The angular for maximum 
intensity and minimum intensity (extinctions) were 
also averaged in a manner to allow for their distri- 
bution throughout the complete circle, yielding a 
grand average value of the first quadrant extinction 
angle, which was based upon all 24 angular readings. 
Through the use of results based on averages of large 
numbers of readings, personal errors were minimized. 
The extinction angle is defined as the angle from the 
fiducial line to the position of minimum reflection. 
If, moreover, polarization figures could be observed 
as they readily were with the monel specimens 
their crossing positions were also averaged into the 
value of the extinction angle for each grain. To 
make these figures visible, the incident light was 
made highly convergent by the use of a high-power- 
ed objective, such as 41x. Moreover, as the figure 
was formed at the back focal plane of the objective, 
the method of observation was modified by the re- 
moval of the eyepiece. Usually the figure could 
then be seen directly or through a pinhole evepiece; 
viewing was often improved by a lens system, such 
as a short-foeus telescope, which allowed focusing in 
the proper plane. 


settings 


4. Results 


4.1. Tin 

The tin specimen was prepared by melting the tin 
in a %-in. quartz ring resitng on a piece of polished 
graphite. It was then air-cooled. Approximately 
15 specimens were made, resulting in only one in- 
stance of a specimen (1144, fig. 2) whose surface 
contained numerous grains large enough (about 1 
mm in smallest dimension) for convenient determi- 
nation of orientation by X-ray diffraction. 

Considerable difficulty was encountered in prepar- 
ing the surface of this soft metal. It was never 
ground on papers but polished with only a diamond 
abrasive, 0 to 2 microns, on microcloth that has been 
boiled for 2 hours before being used, Then the 
specimen was etched in acidified ferric chloride 
reagent [/4]. Figures 3, a and b, show the same 
surface lightly and heavily etched. 

On the surface of the specimen there were 29 
grains large enough to give usable back-reflection 
Laue patterns. The. crystallographic orientation of 
the surface normal for each grain is shown in the 
stereographic triangle in figure 4. The stereographic 
angular coordinates of the optic axis (the c-axis, or 
[001] axis), the azimuth referred to a fiducial line in 
the surface and the colatitude, or zenith distance 
relative to the surface normal, are given in table 1. 
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Ficure 3, Surface of tin 1144 after etching with 5-percent ferric chloride solution. 


A, Light etch, x 50; B, deep etch, x BO. 
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lhe following symbols are used in the tables 

4/—Change of intensity of light falling on sensitive 
tube as grain is rotated from a position giving 
maximum intensity to an extinction position 
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OLATITUDE OF C—AXIS, DEGREES 


y Eextine tion angel The average angle obtained 
by a rotation of the microscope stave trom a 
reference line on the surface of the specimen 
to the position of the first extinetion 

o.— Azimuth of the e-axis (the optie axis) of tin, 
measured from the same fiducial line as that 
used in measuring extinetion angles 

? Azimuth of the projection of the cubie axis 

making the greatest angle with the surface 

norma | 


y= Adjusted extinction angle: y, plus such a mul- 
tiple of 90° as to give a value near ¢, or @,, as 
the case might be 

y Adjusted average of extinction angle and polat 
ization figure crossings when the latter were 
usable 

Ay Difference between the adjusted extinetion angle 


and the azimuth of the active axis; ¥ or ¥ 
minus ¢, or @,. as the case might be 

§6—Colatitude of the active axis: The angle between 
the active axis and the normal to the surface 

a=—Angle between surface normal and [110] diree- 
tion in Al 

8—Angle between surface normal and [111] dire« 
tion in Al 

y Angle between surface normal and [100] diree- 
tion in Al 


Table 1 also includes the photometer data and the 
extinction angles for the lightly etched specimen, 
with the differences between the c-axis azimuths and 
the extinction angles. The average difference was 
6.3°. However, if AJ=3 is arbitrarily taken as the 
value below which valid judgments of extinction 
angle cannot be made, four grains are excluded, and 
the average difference is 4.0 The results show 
that, with the exception of these four dim grains, 
there was excellent agreement of extinction angle 
with azimuth of the c-axis. Figure 5 is a graph of 











Figure 7. Surface of aluminum 1143 after certain treatments. 


vy; C, anodized, & 30; D and E, etched with modified Tucker's reagent, 
respectively 
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versus 6, 
the c-axis for this specimen of tin 
vn was for the equation 





} Al 15 sin* 6, 
vhich no theoretical justification was attempted 
vever, it exhibited reasonable SsvVinmetry and 
ivior at the end points, and the probable erro! 


the data points 
pared with those of points on this curve was less 


he value of the colatitudes of 


\leasurements were also taken from 
14 after the surface had been deeply etched 
age difference between extinction angle and 
S azimuth more than doubled, and the A/ 


tionship to the c-axis colatitude developed so 


specimen 


ch seattet that it could hardly be called more 
na trend 

4.2. Aluminum 
The aluminum specimen, 1143, was a coarse 


ined cast specimen that had been polished by 
nventional metallographic procedure. lts granular 
ucture is shown in figure 6. The intensity meas 
ements were after treating the polished 
rface in four surface 1 was the result of 
etching with Tucker’s reagent |/5]; surface 2 was 
ter repolishing, anodized in an electrolyte recom- 
: mended by Pearson, Machland, and Hay {/6); 
rface 3 was a repetition of the anodizing treatment; 
nd surface 4 was the result of a deep etch with a 
modified Tucker’s reagent [17]. These surfaces may 
be seen in figure 7. To show that the faces of the 
deep pits in this case were parallel to cubic planes, 
the positions of the specular reflections of ordinary 
ght from them were measured with an optical 
5 goniometer [/8]. The results checked satisfactorily 
The crystallographic orientation of the surface 
normal for the aluminum grains is shown in the 
stereographic triangle, figure 8. The stereographic 
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| URE Ss Crystallographic orientations of surfaces normal 


th respect to individual crystals in the polished surface of 
aluminum 1143. 


mbers 1 to 10 correspond to the number tifying the crystals in figure 6 
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rhe curve 


Figure 9 








coordinates, azimuth and colatitude, for the cubic 
axis making the greatest angle with surface normal, 
are given with photometer data from surface 4 In 
Table 3 presents the intensity and orienta- 
which is shown graphically in 


table 2 
tion relationship, 
figure 9 

For the etch-pit conditions, surfaces 1 and 4, the 


avreement between the extinetion angle and the 
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Relation of AI of aluminum grains to the difference 
in angles between surface normal and the nearest [110] and 


[111] ares 
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aizmuth of the cubie axis of greatest colatitude was 
ood The average difference for surface 1 was 4.2 
and for surface 4, 1.9 No doubt the intensitv was 
related in some wavy to the orientation of the surface 


with respect to the vrains crystallography, but no 


empirical relationship could be found that was 
satisfactory over the entire range of the limited 
amount of data at hand For convenience, one 
treatment of this data ts presented in tigure Y 


However, no real functional relationship Is implied 

The qualitative aspects of the reflection of plane 
polarized light from the aluminum with eteh pits 
correlated well with the multiple reflections observed 
from the pit models made with front-surfaced mirrors 
because the faces of the eteh pits are but slightly 
imperfect cubic planes 

\ few limited experiments were made with three 
mirror models. Onewas a single plane front-surfaced 
aluminized mirror with a stem normal to the plane 
One was a “cube corner,” three mutually perpendic- 
ular mirrors, with a stem making equal angles with 
the faces. The third consisted of two mirrors at 
right angles with a stem perpendicular to their line 
of intersection and making equal angles with each 
face. From experiments conducted with these mod- 
els, it appears that there is a useful analogy between 
the orientation of the plane of the reflected polarized 
light and the orientation of the image of a linear 
object multiply reflected by the same mirror system 

Reflection from the single-plane mirror, normal 
incidence, showed no rotation of the plane of polari- 
zation. Nor did reflection from the three mutually 
perpendicular plane mirrors, no matter what the 
position of the model as it was rotated about its 
stem, corresponding to a [111] direction. Study of 
the light reflected from the model made of two per- 
pendicular mirrors, however, was more complex 
When the line of intersection corresponding to a 
cubie axis) of the two surfaces was parallel to or at 
right angles to the plane of polarization of the inci- 
dent light, the plane was not changed upon reflection. 
When the line of intersection was at 45° to the plane 
of the incident light, the plane of the reflected light 
was rotated 90 Or, in general, for a rotation of 6° 
of the model about its stem, corresponding to a [110] 
direction, there was a 26° rotation of the plane. 


These results were implied in the equations pertain- 
ing to multiple reflection of polarized light published 
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by Olwen Jones in 1924 [/9]. This turning of the 
plane of polarization is analogous to the rotation of 
the Image of an object placed before such a rotating 
mirror, as illustrated in figure 10 

In observations with the etched specimen, the 
greatest AJ was observed from grains with thei 
surface normals near 1] 10] In accordance with 
reflection from the two-mirror model, the intensity 
of the light passing through the crossed analyzer 
rose and fell four times per revolution. In the metal 
lographic specimen, the grain with its surface normal 
nearest [111] had the smallest A/, that is, gave the 
smallest increase in light above the faint light at 
extinction. This agrees well with the three-mirror 
model, which did not rotate the plane of polarization 
in any position, and therefore always showed extine 
tion when viewed through a crossed analyzer. The 
results of the experiments with the mirror models 
would lead one further to expect a minimum A/ 
from grains with normals near [100]. This, however 
was not observed to be very pronounced. A possible 
explanation could be that pits of this type were sub 
ject, on the average, to a greater degree of imperfee 
tion because more cubic faces would be involved in 
each such pit than in the two simpler cases. 

The anodized films, surfaces 2 and 3, were reason 
ably light, just heavy enough to give smooth, nearly 
uniform cover. Polarized light, upon reflectior 
from these surfaces was markedly affected, both as 
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extinetions and to striking colors when examined 
hrough the sensitive tint plate. There was no appar- 
correlation, however, between this activity and 
ervystallographiec orientation of the underlying Degrees D 
etal \ film-orientation correlation was reported 6 
wever, by Pearson, Machland, and Hay |/6], em 3 mo rs " 
loving a heavier film than was used in this study : on 


4.3. Monel , ‘ 82 x 


f Che surlaces of the two monel specimens, 724 and 
SS. were polished mechanically and etched with 
two reagents: A, a evanide-persulate aqueous solu- 
on and B, monel contrast solution 20) Kach _ - ft seleee / i ei 
solutIOn Was used twice alternately on specimen 724: and precisior a es . y I seo urfaces on monel 
only the evanide solution was used on specimen 1138 
\s pointed out in section 3.2, orientations of grains 
a ised were determined by the twin-trace method; the 
positions of the surface normal relative to the 
vstallographic directions in the vrTrains are shown in 
stereographic triangles in figures 11 and 12 The 
stereographic coordinates, azimuth and colatitude, 
of the cubic axis of greatest colatitude are given in 
n table 4 for monel 724, along with photometer data 
ind extinction angles based both on photometer 
minima and polarization figure crossings. These 
data were taken from the surface with the second 
vanide etch, preparation 3. Table 5 compares the 
tensity of reflection and the precision of the cubic 
ixis azimuth-extinetion angle correlation for the 
) four surface preparations of this specimen. Table 6 
lists the grains in the four cases in order of decreasing 
A/. It is apparent that the several sets of measure- 
ents of reflection showed much variation, especially 
terms of intensity and its change, AJ. In all four 
ses there was reasonable agreement of extinction ‘ 
cle with the azimuth of the cubic axis of greatest ; 
ititude, but the average differences in the cases s ‘ . " 
nere the cy anide etch was used were less than half 
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TABLE 7 Photomete minima for grains of Vonel 1138 
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as large as those in the contrast etch cases. In all 


cases, however, there were large shifts in ranking the 
grains in terms of intensity, as can be seen in table 7. 
Hence, there was obviously no correlation of intensity 
with orientation, except the vague observation that 
the grain with surface normal nearest the [100] 
direction was always of smallest AJ, whereas the two 
grains with normals nearest [110] were always near 
the top of the list in AJ. 

The other monel specimen, 1138, was etched with 
the cyanide reagent only. The surface normals of 
the grains are located in the stereographic triangle 
in figure 12, and the coordinates of the cubic axis of 
greatest colatitude are given in table 7, along with 
the photometer and microscope data. These results 
are comparable with those from monel 724, showing 
good agreement of extinction position with the 
azimuth of the cubic axis of greatest colatitude, and 
only a hint of correlation of intensity with angles 
between the surface normal and the |110) direction 

It should be stated that the correlations were 
usually poor when the normals to the surface were 
far from the [110] direction. This was observed for 
several grains in a different area of specimen 1138, 
for which data are not presented here. The extine- 
tion angles had a fair degree of reproducibility but 
did not always correspond with the azimuths of the 
cubic axis of greatest colatitude. In these cases they 
often corresponded to the azimuths of the axis most 
nearly normal to the surface. This was in agreement 
with some of the observations reported in the paper 
by Vacher [2], where specific data on a few such 
grains with anomalous extinctions were presented. 
This disturbing behavior did not appear in the case 
of aluminum or tin, 


5. Discussion 


The practical value of the polarizing microscope 
for the determination of the orientation of crystal 
grains in metallographic specimens and the source 
or mechanism of the optical anisotropy will be dis- 
cussed in the light of the foregoing results. 


10 


Value of Polarized-Light 


Measurement: ip 
Determining Crystal Orientation 
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An examination of the results presented aby; 
reveals varying degrees of correlation between cer( yi 
characteristics of reflected polarized light and the 
crystallography of the grains in the reflecting meta] 
surface It is true, however, that a complete orien- 
tation determination of a grain in a metallographic 
specimen cannot be made alone from polarized. 
light observations such as these. The first difficulty 
encountered is the obtaining of a suitable specimen 
surface, that will not cause the pertinent phenomena 
to be masked by the extraneous optical effects of 
scratches, irregular pits, nonuniform films, cold- 
working, inadequate etching, ete. If a nearly 
perfect, lightly etched surface with polarizing 
activity is obtained, the observation with crossed 
polarizing elements of the four extinction angles as 
the specimen is rotated will give within a few degrees 
four possible positions 90° apart for the azimuth 
of the “active” axis. This is the optic axis of a 
uniaxial material, or the axis making the greatest 
angle with the normal for a material showing aniso- 
tropic effects due to etch pits with facets parallel 
to cubic planes; or stated another way, when the 
polarizingly active specimen is in an extinction 
position on the microscope stage, the plane con- 
taining the surface normal and the active axis will 
coincide with either the plane of polarization of the 
polarizer or that of the analyzer 

If the surface quality is such that polarization 
figures can be observed, this ambiguity is cut in half. 
The four settings of the specimen, observed cono- 
scopically, which give the cross figure, correspond to 
the extinction positions. As the specimen is rotated 
the isogyres spread apart and move toward the edges 
of the field in the quadrants containing the projection 
of the optic axis of a uniaxial material [9]. In this 
study the observation of polarization figures with 
uniaxial tin was found to be very difficult, but such 
observations were made quite successfully with 
etched monel, where the position of the cubic axis 
of greatest colatitude governed the behavior of the 
isogyres just as the optic axis is expected to do in 
theory. Locating the projection of the active axis 
in this way reduces its possible azimuth to two 
positions 180° apart. Polarization figures were not 
observed with aluminum either etched or anodized. 
In general, the degree of surface perfection required 
for the formation of these figures is even higher than 
that required for consistent extinction angle obser- 
vations. Apparently the presence of textured films, 
pits, ete., precluded their use with the aluminum and 
tinspecimens. The other quantity measured, change 
of intensity of polarized light as the specimen ts 
rotated, is less reliable in its relation to orientation, 
probably because of its very great dependence upon 
surface condition. Nevertheless, if we can generalize 
from the observations of tin, it appears that for 
uniaxial materials with carefully prepared surfaces, 
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order of grains ranked in terms of decreasing 
vill correspond very closely to the order of those 
ed in terms of decreasing colatitude of the optic 
approximating a fourth-power sine law 
in the other hand, the case for cubie materials is 
complex and apparently is subject to more 
irbing factors. This is shown by the great varia- 
s and inconsistencies of the intensities observed 
monel. The observations with etch pits on 
inum were somewhat encouraging, and lead 
to suppose that if suitable etchants that consist- 
create pits with known preferably 
c, are employed [78], empirical relationships 
ht be found for some isotropic materials. 
is apparent that the measurement of extinction 
e and intensity is not sufficient to permit one to 
iin a complete orientation determination. This 
eedure can be of great value, however, as an aid 
n other methods are being employed, as, for 
mple, the twin-trace method or the etch-pit 
ection method. Ambiguities often arise in twin- 
ce studies: and, with the method of reflections of 
nary light from etch pits, large deficiencies of 
ecision are often an impediment. The additional 
ormation supplied by polarized-light examination 
| be an aid in resolving the difficulties encountered 
such cases. In connection with the back-reflection 
Laue method, polarized-light data will facilitate the 
issignment of indices or the positioning of the 
pecimen so that more readily interpreted patterns 


faces, 


n be obtained. 


52. Cause of the Optical Anistropy of Metals 
Having Cubic Structure 


In section 1 brief reference is made to the question 
of the source of the polarization phenomena. Obser- 
itions reported here showed distinct cases of 
polarization activity due both to a surface film (the 


inodized surface) and to surface contour (the etch 


pits It is not surprising that an effect on the 
reflected polarized light can be due to either mechan- 
ism. The question would still remain, however, 


whether the effect of the film on the light was due to 
intrinsie anisotropy or to underlying metal-surface 
contour. It has been shown here that, when surface 
contour (etch pits) alone was involved, there was 
good correlation of the optical effects and the ery stal 
orientation. The. data from lightly anodized sur- 
faces, however, gave no such correlation. It is there- 
fore evident that the anisotropy of these light films 
was not related in any simple manner to the under- 
lving surface contour, if such contour was a function 
of the erystal orientation, as in the case of the etch 
pits. Hence, it seems probable that these films are 
inherently anisotropic but that their orientation is 
related in an as vet undefined manner to the orienta- 
tion of the underlying grains. 

lhe results from monel were very much less defi- 
nite than those from aluminum, but they showed a 
slight resemblance to the latter. Their lack of con- 
sistency can perhaps be attributed to the variation 


of degree of attack on cubic faces by the etching 
reagents 

When working with tin the conditions were some- 
what different. The fact that the metal itself was 
anisotropic was the controlling factor in its behavior 
Because of the softness of the metal, a very light 
etch improved the results by removing disturbing 
surface features such as polishing seratches, but the 
heavy etch worked against a correlation because of 
excessive irregular attack or an accumulation of an 
obscuring film of etch products 


6. Summary 


Metallographic specimens of tin, aluminum, and 
monel were prepared and the ery stallographic 
orientations of many of the grains in each specimen 
were determined The reflection of plane polarized 
light normally incident upon the specimens after 
various surface treatments was examined by means 
of a polarizing microscope and an electronic photom- 
eter. In some cases convergent light was used and 
polarization figures were observed. Very significant 
correlation was found in most cases between the 
extinction angle and the azimuth of the active axis 
Moreover, in the case of tin, the change of intensity 
of the light from a grain was approximately propor- 
tional to the fourth power of the sine of the angle 
between the optic axis and the surface normal. In 
the case of aluminum and monel, the relationship of 
intensity to orientation was not so clear. The results 
support the hypothesis that the source of the optical 
anisotropy of polarized-light reflection from cubic 
metals is an oriented- surface contour 


The authors acknowledge capable assistance of 
Ellen A. Buzzard in preparing the metallographic 
specimens and express their gratitude to Mrs 
Buzzard and Ruth E. Dowden for aid in the prepara- 
tion of many of the drawings and photographs used. 
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Separation of Iodide, Bromide, and Chloride From One 
Another and Their Subsequent Determination 


Thomas J. Murphy, W. Stanley Clabaugh, and Raleigh Gilchrist 


A method 1 


1 for eparatir 


First. the tocdicde , ) lized to elementary 
und t ed by distillatio 
ifecting the eblorick | e of a controll 
Aft il ! tillation, ehl 
Determinat or the ndividual halide 
t it ul 1 of 
\ t i ra i 


1. Introduction 


inv of the methods deseribed in the literature 
he separation of icicle bromide and chloride 
one another and for their subsequent deter 
tion vield unsatisfactory results In the separa 
of iodide only from a mixture of the three 


thre 


s lodide is usually 


} es published methods do report accurate 
separated by oxidizing it 
sulfate in a dilute solution 
the liberated 
‘ The separation of bromide from chloride is 
we difficult matter In 1928 Wales [2] studied 
nethods reported in the literature and found the 


nitrous acid or ferri 
ifturn acid | 


irl and volatilizing 


method of Winkler [3] to be the most promising 
\\ ler treated a mixture of bromide and chloride 
with sulfuric acid and permanganate and distilled 
o the bromine In 1931 Freeman } reported 
t} method to be unsatisfactory and recommended 

wdification, but in 1942 he [5] reported unsatis 
factory. results with his own revision He stated 
that in this method, as well as in the others studied 


the formation of 
which caused high results, or the forma 
low In the 
best method that he studied for the separation and 
determination of iodide, bromiuc and chloride, th 
recovery of bromide was only 95 percent 

Because of this lack of a reliable method, a study 
was undertaken and a method deve loped by which 
the three halogens can be quantitatively separated 
another Th method vields accurate 
the relative concentrations of 


thre principal eause of error was 
free chlorine 
tion of results 


bromate, which caused 


from one 
results regardless of 


the three halides, 


2. Discussion 


2.1. Separation of Iodide from Bromide and 
Chloride 

weakly acid solution it is possible tO ONDCIZEC 

to iodine with hydrogen peroxide without 

oxidizing either bromide or chloride. In solutions 


containing iodide, bromide, and chloride and hydro- 
y peroxide, the following are the significant re 
which can The free-energy values 


occur 





Ke 
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de. bromide. and chloride from one another 
odu ry ' ( peroxid i i 
Next 1 i; j i " 
con ira tric me ia l i 
le rema lual 
lepend t t f i If r 
i al | il 


in kiloealories per mole, are calculated from values 
reported by Latimer t) as follows 
| HO H 21 l 2H.0O; \f 3.284 
2 H.O I QO, + 2H 21 » S44 
; HO 2H 28 B 2H.O 31.785 
i HO tr io 2H 2b IS.655 
) H.O,+ 2H 4 ( 2H.O; 17.560 
ih Huo ( () oH A s2 S80 
7 HO Hoo O.: 25.220 
The large negative values of AF’ show that all of 
the reactions except reaction (2 tend to take place 


\ta pH of 5, reactions (1) and (2 


the same rate |7 


proceed nt about 
the net results of the reaction being 
peroxide As the 

takes place at oa 
and rate until at 
pli of 1 the rate of reduction of iodine is negligible 
If the iodine formed is 
reaction | 


the decomposition ol hydrogen 


reaction | 


2) at 


acidity is increased 


faster rate reaction a slowe! 
from the solution 

However as the 
to oxidize a small 
reaction 3 Reae 
than reaction ) 8], 


removed 
can vot completion 
acidity is increased it 1s possible 
amount of bromide lo bromine 

tion (4 
the acidity ts favored 
As the bromine and per 
oxide ts rapid in neutral or weakly acid solution and 
it follows that 
would be reduced 
rapidly . with the resultant decomposition of hvdrogen 
peroxide However, if the 
bromine reaction o removed 
“ither through 
reaction reaction | 


is much faster but as 


increased reaction > Is 


rate of reaction between 
the rate of oxidation of bromide slow 


anv bromine formed by reaction (3 


some of 
were 


reaction 7 
produced by 
iodine monobromide 
liberated iodine in 
bromide by reaction } if 
the bromine 
This actually Kor example, 
when 3 ¢ of sodium bromide in 100 ml of water was 
treated with 3 ml of hydrogen peroxide and 2 ml of 
and the solution boiled, bromine 
vapor could be seen distilling ovel However when 
the bromide content. was decreased to 1 ¢ of sodium 


as bromine or 
with the 
before conversion to 
would be possible to lose Some of 


proves to be the case 


phosphor ic’ acid 


bromide in 200 ml of water and the acidity adjusted 
toa pH of l with phosphoric ne id ho bromine could 
be seen distilling over, and a test of the distillate 
2 mg of bromide. This 
error would be negligible for large amounts of iodide 
in the presence of small amounts of bromide How 


showed it to contain only 0) 





ever, to determine small amounts of iodide in the 
presence of a large amount of bromide, the distillate 
obtained must be again treated with hydrogen per- 
oxide and phosphoric acid and redistilled When 
this is done no bromine is detected in the distillate 
When large amounts of both iodide and bromide 
occur together, bromine may be lost as iodine mono- 
bromide, if the acidity is not controlled carefully. 
When enough phosphoric acid ts added at the start 
so that the resultant pH will be 1 after the iodine 
is removed, considerable bromide is lost as iodine 
monobromide But when peroxide Is added first, 
the solution heated to boiling, and the acid added 
slowly us the iodine distills out, only a negligible 
amount of bromide is lost 

Reactions (5) and (6) both take place at about the 
same rate, with the resultant decomposition of hydro- 
gen peroxide [9] At low acidity, only a negligible 
amount of chlorine distills, even from fairly con- 
centrated chloride solutions. Thus under proper 
conditions it is possible to separate quantitatively 
iodide from bromide and chloride. 


2.2. Separation of Bromide from Chloride 


It has been found that nitric acid in suitable con- 
centration will quantitatively oxidize bromide to 
bromine without oxidizing chloride to chlorine 
Kebler [10] and Chamot and Mason [11] both used 
this reaction for the qualitative separation of bromide 
from chloride, but apparently no attempt was made 
to find out if the reaction could be made quantitative 
under proper conditions. Experiments showed that 
in diluted nitric acid (1-+-9)? the oxidation was com- 
plete but slow and that in diluted nitric acid (1+-2 
chlorine begins to distill. It was observed experi- 
mentally that satisfactory results were obtained when 
the range of nitric acid concentration was from (1-+-6) 
to (1+3 The nitric acid concentration used in the 
experimental work was (1+ 5 However, it was 
found that under these conditions a small amount of 
chloride is lost as chlorine. For example, when 20 
ml of nitric acid was added to 100 ml of water con- 
taining | g of potassium chloride and 25 ml of the 
solution distilled into hydrazine sulfate solution, the 
distillate was found to contain 0.5 mg of chloride 
In more dilute solutions of chloride and of acid the 
error is even less. This error would be negligible for 
large amounts of bromide in the presence of small 
amounts of chloride, but a second treatment and 
distillation is necessary in the case of large amounts 
of chloride in the presence of small amounts of 
bromide. This indicates that it is possible to sep- 
arate quantitatively bromide from chloride. 


3. Reagents Used 


Potassium lodide. The potassium iodide used was 
prepared by Mark Swanson, formerly of the Bureau. 
The method of preparation is not known. It was 


* Throughout this paper, whenever the use of an acid is prescribed with no 
indication of strength cr dilution, the acid to be used is the concentrated acid 
Dilutions are indicated by the volumes cf acid and water mixed to prepare a dilute 


reagent Diluted nitric (1+ means a diluted acid prepared by mixing 1 volume 
of the concentrated acid with r volumes of distilled water. 


found to contain 0.0003 percent of bromide 
0.0003 percent of chloride 


Sodium Bromide. A survey of available bron 
sults showed that they all contained chloride.  ‘{y 
sodium bromide used in these experiments was } 
pared from sodium carbonate and hydrobromie a: id 
as follows: 

An approximately 30-percent solution of sodium 
carbonate, prepared from sodium carbonate eop- 
taining 0.002 percent of chloride, was filtered throuch 
a glass filter of fine porosity. Sodium bicarbonate 
was precipitated from this solution by saturation 
with carbon dioxide. The mixture was cooled in ay 
ice bath, and the salt was caught on a porous glass 
filter, washed with cold water, and sucked dry. Thy 
sodium bicarbonate so obtained was redissolved and 
the solution acidified with specially prepared hydro- 
bromic acid, containing 0.0004 percent of chloricd 
to the end point of methyl red. The solution was 
boiled to expel carbon cioxide, evaporated to a small 
volume on the steam bath, cooled, filtered, and the 
resulting salt washed with 95-percent ethyl alcoho! 
This product was found to contain 0.0002 percent of 
iodide and 0.006 percent of chloride 

Potassium Chloride. The potassium chloride used 
was of ACS reagent quality and was found to contain 
0.0006 percent of iodide and 0.0002 percent of bro- 
mide. 

Phosphoric Acid. The phosphoric acid used was 
of ACS reagent quality and was found to contain 
less than 0.0005 percent of combined halides as 
chloride. 

Hydrogen Peroride. The hydrogen peroxide used 
was 30-percent peroxide of ACS reagent quality 
It was found to contain less than 0.0005 percent of 
combined halides as chloride. 

Nitric Acid. The nitric acid used was of ACS 
reagent quality. It was found to contain less than 
0.00001 percent of combined halides as chloride. 

Hydrazine Sulfate. The hydrazine sulfate was 
supplied commercially. It was found to contain 
less than 0.0005 percent of combined halides as 
chloride. 

Silver Nitrate. The silver nitrate used was of 
ACS reagent quality. It was found to contain 
less than 0.00005 percent of combined halides as 
chloride. 

Potassium Permanganate Reagent Solution \ 
survey of commercial samples showed them all to be 
high in chloride content. The following method of 
purification was used: 

Seventy-iive grams of potassium permanganate 
was dissolved in 1,250 ml of water, and 10 ml of nitric 
acid was added. The solution was heated to about 
5° C, and a stream of halide-free air was passed 
through the solution. At intervals the evolved gases 
were passed through a 2-percent solution of hydra- 
zine sulfate for % hour, and this solution was tested 
for chloride with silver nitrate. The treatment was 
continued until chlorine was no longer detected 
This solution, which contained about 0.9 g of 
potassium permanganate in 15 ml, was found to 

| contain less than 0.0001 percent of chloride. 
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Preparation and Standardization of 
Solutions 


Silver Nitrate on a Weight Basis (Approxi- 
mately 0.1 Molal) 


wut 70 ¢ of silver nitrate was dissolved in 4 
; of water. Two 50-ml portions and a 60-ml 
on were measured by means of a weight buret 
portions were concentrated on a steam bath to 
lume of about 10 ml, and 10 ml of 1 N hydro- 
ie acid was added to each The solutions were 
ted to 150 ml and stored in the dark overnight 
precipitates were collected in filtering crucibles of 
porosity, washed with very dilute hydrochloric 

14-999), dried at 120° C, and weighed The 
entration was found to be 0.010625 ¢ of silver 
ram of solution 


4.2. Silver Nitrate on a Volumetric Basis 
(Approximately 0.01 Normal) 
Ihe same procedure as above was followed except 
the 50-ml and 60-ml portions were accurately 
isured volumetrically instead of being weighed 
he solution was found to be 0.10256 N and was 
ited tenfold before use 


4.3. Iodide 


grams of the potassium 


i 


iodide was dis 


\\ eighed portions 


went, 
ved in 2 liters of distilled water 
the solution of about 45 ¢ to 60 @ were diluted to 
ml, and 0.1 molal silver nitrate was added by 
ins of a weight buret to within about 1 percent 
the equivale wee pont The titration was com- 
eted with 0.01 NV silver nitrate potentiometrically 
lhe concentration was found to be 0.0077132 ¢ of 
wlide per gram of solution 
4.4. 
l'wenty grams of the sodium bromide was dissolved 
liters of distilled water. The same procedure 
is used as for the standardization of the iodide 


Bromide 


olution The concentration was found to be 
0.0078208 ¢ of bromide per gram of solution 
4.5. Chloride 
Twenty grams of the potassium chloride was 


ssolved in 2 liters of distilled water. The same 
procedure was used as for the standardizations of 
the iodide and except that the 
solution was cooled in an ice bath during the titra- 
tion. The concentration was found to be 0.0047904 
g of chloride per gram of solution, 


bromide solut ions 


5. Procedure Used for Separating Iodide, 
Bromide, and Chloride 


5.1. Halides Present as Major Constituents 


The distilling apparatus was assembled as in 
ire | Twenty-five milliliters of a 2-percent 
ition of hydrazine sulfate was placed in flask G 

25 ml of a 0.5-percent solution of hydrazine 
ifate was placed in flask 1. About 4 in. of mercury 


s placed in gas washing bottle K and 50 ml of 





| N silver nitrate was placed in gas washing bottle L, 
When this system was connected to a vacuum line, 
the silver nitrate solution served to wash any halide 
from the air, and the mercury maintained the pres- 
sure on the system about 4 in. below atmospheric 
In this manner no halogen was lost through 
the glass joints, as air would have been drawn in 
through a loosely fitting joint 

By means of weight burets and the standardized 
halide solutions, synthetic mixtures of halide salts 
were obtained containing varving amounts of each 
halide. The amounts added were such that a total 
of about 1 ¢ of combined halide salts was present 

The standardized solutioas were added directly 
from weight burets to flask A to avoid error due to 
transfer. This apparatus and method of adding the 
known amounts of halide were used 
throughout all the separations, 
noted 

Separation of Iodide The synthetic halide mix- 
ture in flask A was diluted to about 200 ml. The 
system was closed and connected to a vacuum line 
The vacuum was regulated so that air slowly bubbled 
through the mercury. One milliliter of 30-percent 
hydrogen peroxide was added through delivery 
funnel B, and the solution was heated to boiling, 
Diluted phosphoric acid (14-99) was placed in the 
delivery funnel and showly dropped into the solution 
The addition of 


pressure 


solutions 


unless otherwise 


until iodine began to be evolved 
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phosphoric acid was stopped and the liberated iodine 
removed by distillation. More diluted phosphoric 
acid added and the iodine liberated again 
removed, This addition of phosphoric acid and 
removal of iodine was continued until iodine was no 
longer liberated from the solution. The amount of 
diluted phosphoric acid necessary will depend on the 
quantity of iodide present. Usually, 10 to 25 ml of 
the diluted acid is required. The distillation was 
continued until all of the iodine was distilled into 
flask G® It was necessary to keep the hydrazine 
sulfate solution in flask G warm because todine does 
not react readily with the cool solution. The rate 
at which the iodine distilled was kept slow, so that 
ontv about 25 ml to 30 ml of solution distilled 
About | ml of diluted phosphoric acid (1 +99) and 
0.2 ml of 30-percent hydrogen peroxide were added 
through the delivery funnel and the color of the 
solution noted If a color developed, the iodine 
formed removed by distillation and the treat- 
ment repeated until no color was developed.’ The 
solution in flash | was poured into flask G. Tubes 
F and H and flask I were rinsed and the washings 
caught directly in flask G. The flask and contents 
were reserved for the determination of todide 

Separation of Bromide. The apparatus was 
assembled as above, with another 3-neck flask for G 
Sixty milliliters of diluted nitrie acid (1+ 1) was 
added through delivery funnel B, and the solution 
wis refluxed until all the hydrogen peroxide was 
decomposed With a little practice, this can be 
determined by noting a change in the bubble forma- 
tion The bromine produced was slowly distilled, 
and the distillation continued for 15 minutes after 
the residual solution was colorless and no trace of 
bromine could be seen anywhere in the system 
The rate of distillation was kept slow, and the 
volume in flask A was never allowed to drop below 
120 ml. The solution in flask [ was poured into 
flask G. ‘Tubes F and H and flask I were rinsed and 
the washings caught directly in flask G. The flask 
and contents were reserved for the determination of 
the bromide 

Ne paration of Chloride As the todide and bromide 
originally present had already been removed, the 
only halide remaining in flask A was the chloride 
The flask and contents were reserved for the deter- 
mination of chloride 


5.2. 


was 


was 


re- 


Iodide as the Major Constituent, Bromide and 
Chloride as Minor Constituents 


Separation of lodide. The method of separation 
used was the same as for the separation of iodide 
under 5.1, except that the concentration and method 
of addition of the phosphoric acid were different. 
In this instance 10 ml of diluted phosphoric acid 
1+9) was added directly instead of dropping in the 
more diluted acid (1+ 99 

Separation of Bromide. The apparatus was as- 
sembled as in 5.1, except that tube H and flask I 





The side necks of flask A were gently heated with a toreh to dr r the 
iodine that condensed there 
‘On turning off the heat, it was necessary to open the stopcock on B and to 
disconnect the vacuum to prevent the sclution in G from sucking back 
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were omitted and that 0.5-percent hydrazine sul! 


solution was used in flask G. A piece of platin 
gauze and 60 ml of diluted nitrie acid (1+1) w » 
added to the solution in flask A. The purpose of + \o 
platinum gauze was to catalyze the decomposit 
of the hydrogen peroxide present. The solution ws 
refluxed until the peroxide was decomposed, and 
then for 30 minutes longer. The solution was slow |, 
distilled until no visual evidence of bromine remained 
or until 25 ml had distilled if no bromine was noted 
Tube F was rinsed and the washings caught direct 

in flask G. The flask and contents were reserved 
for the determination of bromide 

Separation of Chloride If more than 0.1 percent 
(1 mg) of chloride was present, the chloride was det: 
mined without further treatment If less than 0.) 
percent (1 mg) of chloride was present, the following 
method was used: 

The apparatus was assembled as in the separation 
of bromide. Fifteen milliliters of potassium perman 
ganate reagent solution was added to flask A through 
the delivery funnel. The solution was heated to 
boiling, refluxed for 30 minutes, 25 ml slowly dis- 
tilled into flask G, and this solution saved for the 
determination of chloride. This isolation of the 
small amount of chloride is advantageous because 
the chloride is thus concentrated in a small volume 
of solution that is free from substances such as nitri 
acid that might have a solvent effect on the silver 
chloride 


5.3. Bromide as the Major Constituent, Iodide and 
Chloride as Minor Constituents 


Separation of Iodide. The apparatus was assem- 
bled as in the separations of bromide and chloride in 
5.2. Using this modified apparatus, the iodide 
was separated as in 5.2 except that at the end of this 
separation a second apparatus similar to the first was 
set up in which flask G of the first apparatus became 
flask A’ of the second one. One-half gram of dipotas 
sium hydrogen phosphate, 2 ml of 30-percent hydro- 
gen peroxide, and 10 ml of diluted phosphoric acid 
(1+9) were added to this flask. The resulting solu- 
tion was heated to boiling, refluxed for 30 minutes, 
and 25 ml was slowly distilled into flask G’. The 
solution in this flask was saved for the determination 
of iodide. The solution in flask A’ was poured into 
flask A. 

This is, in effect, a second treatment and separa- 
tion of the iodide, necessary because of the small 
amount of bromide that comes over in the first sepa- 
ration. The role of the dipotassium hydrogen phos- 
phate is not clear, but halogens appear to distill from 
solutions that contain it faster than from those that 
do not 

Separation of Bromide. The method outlined for 
the separation of bromide in 5.1 was followed, except 
that a piece of platinum gauze was added to hasten 
the decomposition of the peroxide, and 30 ml of con 
centrated nitric acid was used instead of 60 ml of 
diluted nitric acid (1+-1). 

Separation of Chloride. The procedure outlined 
for the separation of chloride in 5.2 was followed 
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Chloride as Major Constituent, Iodide and 
Bromide as Minor Constituents 


Separation of Lodide The iodide was separated us 
5.3, except that a second treatment and separation 
re not hecessary 
Separation of Bromide The method outlined for 
separation of bromide in 5.2.was followed, except 
it at the end of this separation a second apparatus 
ilar to the first was set up in which flask G of the 
st apparatus became flask A’ of the second appara- 
s, One-half gram of dipotassium hydrogen phos- 
sate. 50 ml of water, and 15 ml of nitrie acid were 
ided to the solution in flask A’ 
dure followed as in the first separation. This 
cond treatment served to separate the bromide 
om the small amount of chloride that distilled inh 
w first separation After the second separation 

solution in flask A’ was poured into flask A 

Si paration of hlor ude The procedure outlined for 
he separation of chloride in 5.1 was followed 


and the same pro- 


5.5. Halides Present as Minor Constituents 


Se paratior of lodide The procedure outlined fol 
thy separation of iodide In O.3 was followed, except 
that a second treatment and separation were not 
necessary and that 0.5 ¢ of dipotassium hydrogen 
phosphate was added to flask A before the separa- 
tion 


Se parat ol of / 


Bromide The procedure outlined 


for the separation of bromide in 5.2 was followed 
Separation of Chloride The procedure outlined 
) 


for the se paration of chloride in 5.2 was followed 


6. Determination of the Separated Halides 


In the above procedures the se parated halides wer 
determined by two different methods, depending on 
the concentration of the halide. When the halide 
was In excess of 1 mg, a potentiometric titration 
method was used, and when present in less than this 
amount, a turbidimetric method was used 


6.1. Potentiometric Method 


The potentiometric method was based on a varia- 
tion of the method of Lange and Schwartz [12]. The 
apparatus consisted of the following items: Vacuum- 
tube voltmeter with full seale of 100 mv; silver 
clectrode: calomel] electrode: salt bridge, consisting 
of a glass tube fitted with fine porous frits at the 
ends and filled with a solution containing 1 part of 

veerine and | part of a saturated aqueous solution 

ammonium nitrate; and magnetic stirring appa- 
rautus 

The 3-neck flask containing the separated halide 
vas placed on the stage of the magnetic stirrer, and a 
leflon-coated stirring bar was placed in the flask 
lhe salt bridge was inserted into a side neck of the 
-neck flask so that it reached below the surface of 
the solution. The silver electrode was inserted 
through the other side neck in a similar manner 
The other end of the salt bridge was immersed in 
he calomel electrode, and the silver and calomel 


ectrodes were connected to the voltmeter When 





chloride was being determined, it was necessary to 


cool the flask in an ice bath during the titration 
During all determinations it was necessary to shield 
the flask from light. This was accomplished by 
enclosing the assembly in black paper 

Standardized silver nitrate was added from a 
weight buret until the end point was approached, 
usually to within | percent of the theoretical amount 
required. The titration was completed by using 
0.01 N silver nitrate. The change in electromotive 
force produced by the addition of known amounts of 
the silver nitrate solution was recorded, and the end 
point was determined by the method of plotting 
Ak Ai with re spect to the volume, of solution 
added From the weight of the first silver nitrate 
solution and the volume of the second, the amount 
of halide was calculated 


6.2. Turbidimetric Method 


When the halide content was small (less than 1 
mg), a turbidimetric method was used after suitable 
dilution The method was as follows 

The solution in flask (; was diluted sous to contam 
from 0.01 to 0.05 mg of halide in 50 ml of solution 
Controls containing varving amounts of halides (e. ¢ 
0.01, 0.02, 0.03 me of halide) and the same amount 
of hydrazine sulfate as in 50 ml of the diluted sample 
solution were prepared, One milliliter of nitrie acid 
and l mil of 0.1 \ silver nitrate were added to the 
50-ml portion of the sumple solution and to the 
controls The turbidities produced were compared 
visually in Nessler tubes 


7. Results 
Table 1 shows the results obtained for the separa 
tion and determination of known amounts of halides 
when at least 1 mg of each halide was present. Th 
separations were all performed according to the pro 
cedure outlined in section 5.1 The determinations 
were all by the potentiometric method 
Table 2 shows the results obtained for the separa 
tion and determination of known amounts of halides 
when one of the havlicle Ss Wis present us 2 Tabor cone 
stituent, with the other two halides as minor con 
stituents. The separations were all performed ac 
cording to the procedures outlined in sections 5.2 


5.3., and 5.4. The major constituent was not de 
termined because the results shown in table 1 prove 
that large amounts of halides can be recovered. The 
minor constituents (halides) were determined by the 
turbidimetric method 

Table 4 shows the results obtained for the isola 
tion and determination of small known amounts of 
the halides in the order of 1 mg or less of each. The 
separations were all done as in the procedure given 
in section 5.5. The determinations were all by the 
turbidimetric method 

The results given in tables 2 and 3 were corrected 
for any halide present in the reagents 

As shown in table 1, there was a standard deviation 
of 0.0002 ¢ for the iodide determinations, 0.0002 ¢ 
for the bromide determinations, and 0.0003 ¢ for 
the chloride determinations. Because all determi- 
nations reported in tables 2 and 3 were made by a 





constituent 


ls obtained when halide was a n ajor 


Panu | he 


I Hy 
‘ 
i } \ Preset t i | 
. TT ‘ ‘ iT 
: x44 “4 uy 271 ‘ ooo 
Kl , kh i ; TE 4 7 on 
7. un wire til O4 
, f ‘ 13s 000! f 000 
Kl ! Ke ‘ ' ‘ PT Te reir 4 000K 
) ” on ue 1508 o000 
Ki, wo i Kf - y ” wo) we io mol 
ia His Mal) sO TL 
KI iB Ki anne iam - NT nr O00 
; f w7s " nn 08 wU8 0000 
KI vals nm KOI wre wine LL ol we 
% % TT was wn rT 
| ‘ 
KI iN ; uN WINS ONS PTT 
, ‘ ‘ x 
KI, ® ' hk - : ‘ 
wis nis wun ‘ ] ue 
PARLI 2 Re alts obtair i her 0 OF the hatides were minor constituents 
ajo cor fil ent 
| } 
" 
\ I \ ! A 
Bn Ki 
I { 
‘ 
i ( Kl 
" K 
( 
iN ‘ tkKl 
« I " Tt Kt 
; 
" n " ( k 


visual turbidimetric method, standard deviations Panie 3. Results obtained wi 
were not calculated, as they would be meaningless 
For all practical purposes, the recovery of added 
halides in these determinations may be regarded as batho 
kK 


LOO percent ‘ 
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dewpoimt measurements made 


Relative Humidity-Temperature Relationships of Some 
Saturated Sait Solutions in the Temperature 
Range 0° to 50° C' 


Arnold Wexler and Saburo Hasegawa 


The relative humidity-temperature relationships have been determined in air i eq 
librium witl aturated salt solutior of lithium ehlorice Lif H.O:; magn mm chloride 
Mu 6HLO | n dichromate Na Cr.0,-2H_0 magne im trate Mei NO OLLO 

dium chloricd Nal ammonium fate NH g)SOg: pota im nitrate, KNO.: and 

tas im fate KSO, ver a temperature range of O° to 50° ¢ isme the lewpoint 

thod rhe relative humidity is a continuous funetion of temperature, and, except for 

Miium chloride, is monotonic The curve for sodium chloride increases from 74.9-pereent 
relative humidity at 0° C to a maximum of 75.6 percent at 30° C and then gradually ce 
crease to 74.7 pers t The maximum change in relative humidity with temperatur 
about 15-percent relative humidity as tl temperature increase from O° to 50° ¢ occur 
with saturated salt Nutions of sodium dichromate and magnesium nitrate 


1. Introduction 
Saturated salt solutions are very useful in produ 
r known relative humidities, principally for testing 
nd calibrating and hygro 
I iphs at temperatures above 0° C The saturated 
made slushy mixture with 
stilled water and chemically pure salt in a glass or 
nameled trav and a sealed metal or 
lass chamber. When equilibrium conditions, usu 
ly hastened by forced air circulation or stirring, are 
tained, the chamber space is at a constant relative 
pumiclity Some saturated salt solutions produce 
elative humidities that are roughly independent of 


selected hvgrometers 


ilt solution is up as a 


is enclosed in 


emperature 

The equilibrium values of relative humidity of the 
aturated solutions of several salts used in calibra 
on and testing are listed in NBS Circular 512 [1] 
nd are based partly on vapor-pressure data given 
n the International Critical Tables [2] and partly on 
at NBS The relative 
humidity temperature relationships ol these. as well 
s several other, saturated salt solutions have been 
edetermined and are the subject of this paper 

The measurement of vapor pressure or vapor pres 
sure lowering of a salt solution may be made in sev 
The differential 
he determination of the difference in vapor pressure 
wtween the solution and solvent. using a& sensilive 
lifferential manometer. In the dynamic method, the 
boiling point of the solution is determined under re 
The transpiration, or gas-saturation, 


al ways methods are based on 


wed pressure 
vethod involves the gravemetric measurement of the 
ater-vapor content of an inert gas saturated by 
massage over or through the salt solution, The dew 
romnt method consists in reducing the temperature ot 
mirror surface until condensation occurs 


raper 
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The dewpoint method was utilized in this investi 
With this method, 
readily be made in the presence of an air atmosphere 
under conditions similar to those occurring in the 
solutions for humidity-control 
Furthermore, other than using chemically 
pure salts of reagent grade and distilled water, no 
were taken Thus the 
and water was alr 


vation measurements could 


use of saturated salt 
purposes 
salts 


special precautions 


contained trace impurities, the 


saturated 


2. Description of Apparatus 


The experimental setup shown in figure 1 was used 
to determine the 
saturated salt solution 


Wolff flask, A, 


equilibrium vapor pressure of a 
The bottom 2 liter 
was filled to a depth of 2 inches 
with a saturated solution of a salt, B 
a slushy mixture with a pure reagent grade chemical 
and distilled water. A dewpoint apparatus, C 
inserted through one side neck of the flask: an 
stirrer, D 
the flask: and two copper-constantan thermocouples, 
E and KF, were inserted through the other side neck 
of the flask The Wolff flask 
liquid bath, O. Automatic temperature control of 
the bath was achieved by a system in which thermis 
bridge circuit, detected the temperature 
heater, H, supplied heat in 
unbalance sclow 


of a 
made up as 


was 
air 
was inserted through the central neck of 


was immersed in a 


tors, In a 
unbalance and a 
the 
a cooling system, |, was 
bath temperature 

The shaft of the air stirrer, D, was supported in a 
close-fitted bearing and driven, through a belt and 
pulley, by an external remotely positioned motor 

Thermocouple, I, measured the air temperature 
within the Wolff flask, thermocouple, F, measured 
the temperature of the saturated salt solution, and 
thermocouple, H, measured the temperature of the 
liquid bath. These thermocouples were made with 
cotton-covered Bakelite-insulated To insure 
that spurious would be 


portion lo room tempera 


ture, used to reduce the 


wire 


no electromotive force 
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introduced breve tise of electrolv tt aetion each 


junction was coated and protected with polystyrene 
The thermocouple leads were covered wit! 


eement 

grounded equipotential shields \ precision po 
tentiometer Was used fo mhensure the thermocot ple 
electromotive forese with nu precision ol better than 


1 microvolt 


A light was projected on the dewpoint mirror and 


the mirror was viewed, through a telescope, from 
position P, external to the liq | bath 
Details of the dewpoint apparatus are shown in 


} inches in diameter and 3 


figure 2. A brass cup, A 

inches deep, was attached to one end of a inch 
outside-diametetr copper tube, B inch in wall 
thickness. A dewpoint mirror, H, was soldered to 
the other end of the copper tube. This mirror 
L-inch-outside-diameter and inch thick, was 


machined from copper, with a recess equal to the 
wall thickness of the copper tube so that it would 
fit snuggly in the end of the tube. The reflective 
surface of the mirror was first ground and polished, 
then plated with chromium, and finally lapped to 
within one-half fringe of optical flatness. A l-mm- 
diameter hole was drilled in the center of the mirror, 
from the rear, to within about \, inch of the front 
reflective surface \ copper-constantan thermo- 
couple, C, B&S gage No. 30, was inserted into the 
hole and soldered to the mirror 

An insulated-wire heater, E, was wound around 
the copper tube close to the mirror end. This heater 
a variable transformer from a 60- 
Glass-wool insulation, G, 


was controlled by 
evele a-c powell 
was W rapped around the copper tube and retained in 
position by a Bakelite tube, F, ‘-ineh outside diam- 
eter A galvanometer mirror, I, attached to the 
Bakelite tube by a metal bracket, J, allowed con- 
venient viewing of the dewpoint mirror 

The addition of dry ice and aleohol to the brass 
A, figure 2, served as a cold source for cooling 


source 


eup, 


SS ES SS See 


the dewpoint mirror, H, through conduction along 
the copper tube, B. By decreasing the voltage 
applied to heater, E, the temperature of the dewpoint 
mirror could be lowered gradually until dew or frost 
was detected. By increasing the voltage applied to 
heater, FE, the temperature of the mirror could be 
raised slowly until all dew or frost disappeared. The 
dewpoint mirror temperature corresponding to the 
instant of appearance and disappearance was 
measured by thermocouple, + 


3. Experimental Procedure 


The procedure employed in making observations 
was to adjust the liquid-bath temperature to some 
desired value and then permit the Wolff flask, with 
its contents to come to temperature equilibrium 
After equilibruim had been established, an observer 
would manipulate the dewpoint heater control and 
simultaneously view the dewpoint mirror until he 
had obtained five successive appearances and dis- 
appearances of dew or frost on the dewpoint mirror 
A second observer would then repeat the process 
under the same operating conditions. Occasionally, 
the two observers would take repeat sets of observa- 
In this fashion, for any one salt, data were 


‘ 


tions. 
obtained at 10-deg. intervals from 0° to 50° ¢ 


Ld ie 
~ For purposes of this investigation, it was assumed 
that the dewpoint corresponded to the mean of the 
observed temperature at which dew first appeared 
and then disappeared. For a given observer, each 
dewpoint therefore corresponded, usually, to the 
average of a minimum of 10 temperature observa- 








4. Results 
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tions The dewpoint for a salt at a test tempera 
ture, was taken as the mean of the dewpoints ob 
tained by ene h observer 

The ambient temperature was assumed to be the 
mean of the air temperature and the salt temperature 
within the Wolff flask As these temperatures were 
relatively stable for appreciabl periods ot time 
only one such set of temperature readings was ob 
tained by each observe! lor his corresponding set ol 
dewpoint observations 

The relative humidity in equilibrium with the 
saturated salt solution ts viven hy 


RIT" 100 


where « is the vapor pressure ol pure water at thre 
dewpoint temperature, and ¢, is the vapor pressure 
of pure water at the ambient temperature. The 
vapor pressures given in the Smithsonian Tables [3] 
were used in this computation 

The mean relative humidity at each observed test 
temperature is presented in table | for eight different 
saturated salt solutions ranging in relative humidity 
from about 11 percent to 99 percent The data for 
each salt were plotted and the best smooth curve 
obtainable by eve was drawn through the plotted 
points. The faired values from this smooth curve 
are given in table 2. They represent the best 
estimates of the relative humidities obtained with 
these salts nt the st lected temperatures 








The plots of relative humidity versus temperature 


vield curves that are continuous and, generally, vary 


little with temperature ovet the temperature range 


These plots are shown as solid-line 
The maximum variation of 


of O° to 50° C 
curves ip figures 3 to 10 


relative humidity with temperature oceurs with 
saturated solutions of sodium dichromate and 
magnesium nitrate For these salts there is an 
absolute decrease in relative humidity of about 15 
percent as the temperature increases from 0° to 
50° (. The seatter of the mean experimental 


values above the smooth curves is given in table 3 
scatter of the data oecurs with sodium 
+ 0.6 percent relative humidity. For 
the seatter is equal to or less than 


The vreatest 
chloride and Is 
the other 

0.3 pereent relative humidity 
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5. Discussion 
5.1. General 
There are several mivantages to the use of the 


method in determining 
humility 
The experimental 
requiring only a temperature controlled liquid bath 


the equilibrium 
saturated salt 
relatively 


dewpoint 


relative over a solution 


setup can be simple, 
for maintaming the salt solution and its containing 
i dewpoint 
the 
atmos 


vessel at fixed temperature and 
The dewpoit can be 


other 


some 
detector measured ip 


presence ofan au ol noncondensible 


phere \ethods that depend on the direct measure 
ment of vapor pressure, either absolute or differen 
tial require careful elimination of all gases and 
Vapors except the water Vapor from above the salt 


solution ais well ns the elimination of any dissolved 
the solution because these would produce a 
error in the Vapor The absolute 
the Vapor 
creases With decreasing dewpoint for a given accuracy 
in the dewpoint measurement. For example, if the 
dewpoint is measured with an accuracy of 0.1 deg C, 


ruses if 


SCTIOUS pressure 


necuracy. of pressure measurement in 


the corresponding accuracy in vapor pressure is 
0.008 mm Hg at —20° C, 0.084 mm He at 0° C, 
0.109 mm Hg at 20° C, and 0.296 mm Hg at 40 °C 
The accuracy, in terms of percentage of vapor 
pressure, is OLS percent at 20° ©, 0.7 pereent at 
O° ©, 0.6 percent at 20° C, and 0.5 percent at 40° C 
Thus, the pereentage accuracy in) vapor-pressure 


determination by the dewpoint method, for a 
dewpoint accuracy, is roughly independent of dew- 
pomt or vapor pressure. The dewpoint method is 
an indireet method; it vields a temperature from 
which the vapor pressure is obtained by recourse to 
tables. One of the inherent limitations of the 
method lies in the basic accuracy of the tables. For 
most work, and particularly for this investigation, 
the standard tables are considered to have adequate 
Hecuracy 


viven 
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5.2. Errors 


In determining relative humidity by the dewpo: 
method, the principal errors arise in the meas) 
ment of the dewpoint temperature and the ambie»t 
temperature The precision of the temperat 
measurement, using calibrated single-junction copy 
constantan thermocouples and a calibrated precision 
potentiometer, was 0.0L deg C, but the accuracy y 
probably no better than 0.05 dee C As (1 the th 
mocouple for measuring the dewpoint was imbed 
only several hundredths of an inch below the surfac 
of the dewpoint mirror, (2) the body of the mirro, 
the rates of heating and cooling 
were uniform and slow, it seems reasonable to assume 
that any temperature differentials within or across 
the mirror were negligible, say, no more than 0.01 to 
0.02 deg C Because of the relatively low thermal 
conductivity of liquid water, the temperature differ 
ential across the thickness of a dew film was possibly 
of greater magnitude. In cooling, the mirror was 
colder than the exposed dew surface and in heating 
warmer. ‘Two tend to reduce the 
importance of this. First, by averaging the tem 
peratures at which dew is first observed to appear 
and then disappear, the errors due to the differentials 
across the dew, because they are of opposite sign 
tend to cancel Second, by keeping the thickness 
and quantity of dew small, the differentials can be 
reduced to a minimum 

The manually operated and visually 
dewpoint hygrometer required a certain degree of 
skill for successful After some experience in 
its use, it was possible to control the rate of heating 
and cooling so that the mirror temperature could 
be followed easily with a potentiometer, the appear 
ance and disappearance of dew could be detected 
readily, and the thickness of the dew deposit kept 
as small us desirable An observer could repeat a 
series of, say, five dewpoint determinations, at any 
one time, with an average deviation of a singh 
determination from the mean of 0.04 deg C, with an 
average deviation of the mean of 0.01 C, and with 
an average range of 0.13 degC. The average differ 
ence in temperature between the appearance and dis 
appearance of dew was 0.5 deg C, whereas the maxi 
mum difference ever observed was 2 degC. The dew 
point appeared to be independent of the various 
observed differences in temperature between the 
appearance and disappearance of dew 

The dewpoint (the average, usually, of five repeat 
determinations) obtained by any one observer had 
an average deviation from the mean of two or more 
observers of 0.07 deg C The average deviation of 
the mean dewpoint of all observers for a salt solution 
at any temperature was 0.04 degC. The mean range 
in dewpoint for all of the test conditions for two or 
more observers was 0.16 deg C. 

Because the dew was visually detected, it is likely 
that different observers used different criteria fo! 
defining the appearance or disappearance of dew 


was copper, and (3 


it was factors 


observed 


ise 





















may account to some extent for the spread in 
ypoint from one observer to another 

Below 0° C, either frost or dew 
| y be deposited on the 
' In all cusses 


{ ns 93° C) the deposit first detected was dew 
he mirror temperature was not permitted to fall 
far below the dewpoint, then dew remained on 
mirror throughout the determination If the 
ror temperature dropped appreciably below the 
ypoint, the dew changed to frost 
The ambient temperature was taken as the mean 
perature of the salt solution and au These 
temperatures differed from other 
nuse of the heat sources and sinks associated with 
he setup The fan used for stirring the air within 
chamber introduced heat through agitation 
Lo duction along the shaft added or abstracted heat, 
whether the 


supercooled watel 
surface of the dewpoint 
which included dewpoints aus 








er two each 


; vending 








on room temperature was 
ove or below that of the test chambet The 
ny ' ypoint hvgrometer served as a serious heat sink 
8 by the very nature of its operation, it has to 
maintained below test-chamber temperature 
hi Although the insulation surrounding the conduction 
n | aided appreciably ino reducing the magnituck 
al f the heat loss there was still enough loss to produce: 
ils temperature differential between the salt: solution 
n air. The test-chamber temperature may there 
5S ve have an average uncertainty, because of this 
bi ferential, of 0.09 deg C 
; The root-mean-square uncertainty in the dewpoint 
ed : ermination is 0.07 deg C and in the ambient (test 
of hamber) temperature ts 0 lO deg C The correspond 
in incertainty in relative humidity varies from a 
i maximum of 1.2. pereent for the high-humidits 
Id rated salt solutions to 0.2 percent for the low 
I imiditv saturated salt solutions 
od 
pl 
5.3. Comparison of Results With Those of Other 
Ly lavestigators 
rs 
in The NBS results are compared with those of othe: 
th nvestigators in figures 3 to 10, in which the NBS 
' esults are shown as solid-line curves. The results 
Is of some workers were reported for a limited tempera 
1 ' ture range or, occasionally, for only one temperature 
\ Wherever necessary, results that were presented in 
us terms of vapor pressure were converted into relative 
he humidity 
lt is apparent from these figures that there is a 
at certain degree of scatter among the values reported 
ud hy othe researchers aus well as some disagreement 
r between the NBS results and those values. In 
of neral, these other results fall within a band of 
yn values that is, on an average, roughly within + 1‘- 
ng percent relative humidity of the NBS results 
ar 


\ 5.4. Use of Saturated Salt Solutions 


When saturated salt solutions are emploved for 
midity control, experience has shown that certain 
must be observed in order that the 


| eautions 
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used without the need of 


theoretical values may be 


measurement It is necessary to enclose the satu 
rated salt solution in a sealed chamber The cham 
ber and the fixtures therein must be made of non 


hyvgroscopi materials, preferably metal or glass else 
the time required for 
achieved may be very 


humidity equilibrium to be 


vreat, sometimes of the ordet 
The chamber 
should be 
equilibrium It 


weeks salt solution, and 


ambient ai 


of days o1 
brought to 

the salt 
occupy as large a surface area as possible and for 


temperature 


is desirable for solution to 


some means of air ventilation or circulation to be 
provided In the latter regard, if at all) possible, 
the motor that drives the fan or blower should be 


external to the chamber Otherwise, the heat «is 
sipated by the motor will gradually raise the internal 
chamber t and introduce uneer- 
tainty in the equilibrium relative humidity The 


time required for humidity equilibrium to be reached 


nperatvure some 


with saturated salt) solutions depends on several 
factors 1) the ratio of free surface area of the 
solution to chamber volume ») the amount of an 
stirring, and (3) the presence of hygroscopic mate 


rials 
salt solution may in rease the rate at which equilib 
riuth is achieved 


It is conjectured that agitating the saturated 


As ideal conditions are rarely obtained in practice 
it is probable that the theoretical values of relative 
humidity are reached In 
saturated salt should not be expected to 
the relative humidity than 
l-percent relative humidity of the theoretical values 


seldom general lst 


solutions 
control 


to closet about 
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Disintegration Rate of Carbon-14 





R. S. Caswell, J. M. Brabant,' and A. Schwebei 


The energy en ion rates of © samy 
ionization chamber From the energev emiss 
through knowledge of the average beta-r 
earlier data on the isotopic abundance, a v: 

obtained 


In view of the large discrepancies existing between 
ous determinations of disintegration rate (and 
nsequently of the half-life) of carbon-14 [1]? a 
isurement by an independent method has seemed 
sirable. In the present work an extrapolation 
umber measurement of the energy emission of a 
sample is combined with the average energy per 
sintegration from beta-ray spectrometer measure- 
nts to \ ield the disintegration rate The half-life 
then determined by using previously reported 
isurements of isotopic abundance [1] The dis- 
tegration found by the relation 2, —n&, 
here 2 is the energy produced per second per gram 
of material, n is the number of disintegrations per 


rate 1s 


cond per gram, and / is the average energy per 
sintegration. Using the number of C™ atoms per 
from mass-spectrometer measurement) the 
O.693n 


ram, “, 
half-life is determined by the relation 7% 

The type of extrapolation chamber used here has 
previously been used in combination with 4-pi beta 
counting to determine the average energy of beta-ray 
spectra [2, 3]. The excellent agreement between 
average energies determined in this way and average 
energies calculated from spectrometer data or from 
theory (using experimental values for 
demonstrates that the extrapolation chamber 
accuracy, Nuclides 
previously studied range in energy from Ca® (4 
0.075 Mev) to Y*® (£=0.895 Mev). 

The theory of the extrapolation chamber has been 
discussed elsewhere [2, 4]. The energy production 
rate for an air-cavity water-electrode chamber is 
given by the Bragg-Gray cavity theorem [5] 
l JnWarpm, Where J, is the number of ion pairs 
formed per gram of air per second, p, is the mass 
stopping power of the water relative to air, and W,,, 
the average energy required to produce an lon pair 
inair. In the present experiment, one electrode was 
aluminum and the other was a dilute, thoroughly 
mixed water solution of C™ (as Na.CO,). In this 
cause } a 2S Wate BR, where B Is called the “hack- 
scattering correction factor.’ B has been evaluated 
experimentally by comparison with a water back- 
cattering electrode (solidified with agar), and by 
xtrapolation chamber studies of the variation of 
ickseattering with atomic number [2]. B is found to 


hbeta-decay 
/ 


does measure , a to 


ax 


or¢ vod 


of California, Berkeley, Calif 
ndicate the literature reference t the en 


w at University 


in brackets 


le have breve Th if d \ na eX i ila 

i rate tl lisintegration rates are determined 
av energy emitted per disintegration Fron 
ilue for the half-life of C™ of 5,900 250 vear 


be remarkably iIndepe ndent of beta energy 
by less than 2% from P and is taken as | 
for CC". The relative stopping power of water to air, 
Pm, iS taken as 1.13, which is the ratio of the number 
of electrons per gram of the two media (1,11 
rected for the difference in stopping power caused by 
the effective ionization potentials of water and air 
Failla and Rossi [6] have reported a value of 
$2.5 ev per ion pair for S*®, which should be very 
close to the value for C™ because both nuclides have 
the same spectrum shape and nearly the same beta 
energy. This value agrees with that calculated by 
Wang [7] from the formula of Gerbes by averaging 
over the energy of the from the initial 
energy until brought to For this report, we 
take W,,,=32.5 ev per ion pair. A recent review by 
Binks [8] of a large number of measurements of W,,, 
indicates 33 ev per ion pair may be preferable. If so, 
the results of this experiment should be changed 
accordingly 

The average beta energy of C'™ was numerically 
calculated to be 49.7 kev, using 155 kev for the maxi- 
mum energy [9] and considering the beta spectrum 
as “‘allowed’’ [10] 

The amount of C™ in the water solution, n, 
determined by comparison with standard ampoules 
prepared by Manov [1,11] and on which 
spectrometer measurements have been made i four 
laboratories. The standard sample is taken as having 
3.13210" atoms of C™ per milliliter of solution, 


varving 


125 


tos 


cor- 


electrons 


rest 


was 


Mlwass 


> Om 
The value of n for the experimental sample is deter- 
mined by evolving CO, from both standard and 
experimental solutions and observing the relative 
ionization currents in a CO,-filled ionization chamber. 

Three runs were made, giving values of n of 12.55 
ue/ml, 4.53 we/ml, and 1.49 ue/ml, respectively. The 
corresponding values of n were 1.2510" atoms/ml, 
1.5310" atoms/ml, and 1.47 10" atoms/ml, respee- 
tively. These runs yield values for the half-life 
of C™ of 5,900 years, 5,940 vears, and 5,840 years, 
respectively. The third run somewhat lower 
accuracy than the first two because of the low activity 
of the sample 

In terms of the disintegration rate values assigned 
to the standard ampoules by Manov and Curtiss 
[11] based on CO,-CS, gas counting, the disintegra- 
tion rates, #, for the present samples should have 
been 13.8 pe/ml, 5.01 pe/ml, and 1.62 ye/ml, respec- 
tively. These values are systematically about 9 
percent higher than the values obtained in the present 
work. This comparison (which is independent of 


is of 
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isotopic abundance measurements) shows a dis 
crepancy of 9 percent 

This discrepancy is of the same order as the un 
certainties in the gas-counting method itself as 


shown by disintegration-rate mtercomparisons [1] 


and indirectly by the half-life determinations (tabl 
| The recent work of Crane [12] suggests the 
possibility of multiple counts in CO,-CS, counters 


clue to producti m of pulses by both electrons and 
hevative 


Table | shows a 


comparison of recent half-life 


values obtained by gas counting, calorimetri 
mensurement, and the present method. “nn” was 
determined by mass spectrometer measurement in 
all cuses except for the calorimeter In tae ealori 
metric measurement Wwe have used our valu 19.7 
ke for the average energy of the beta spectrum 
Probable errors in the present measurements are 
taken as: WH } percent; J | pereent, B, 
LO percent; p | pereent; 7 1.6 percent 


giving an over-all probable error in 


/ | pereent 
the half-life of about 250 vears 


rawr 
Hawk "i 1 

‘) \I ~ 

te Ni ~ 

{ \l “” 

' ‘ 
Ihe present value Is ih best auvreement with the 


calorimeter 
energy 


Inenstrements 
emission rather than a direct disintegration 
determination, The present not in 
good agreement with either group of gas-counting 


5500 vears and 6.400 vears 


which also depend upon 


rate value is 


measurements (about 


In view of the excellent beta-ray spectrometer data 


on C' it appears very unlikely that any uncertainty 
in / can account for the difference. No large error 
should) be present in (VW pm PB) because this 


quantitv. has been checked in’ the 


independently 


has an average energy about 1.5 tir o< 


average measurements 


previous energy 
ample, Ca 
that of C'™ and also has an “allowed” 
With the same value of (W p R) used here 
disintegration rate of a Ca® source determined |), 
extrapolation chamber coincided with a value dk i 
mined by 4-2 counting within | percent 

In conclusion, a reinvestigation of disintegra 
rate measurements by several independent: methods 
appears desirable. Considerable emphasis shouk 
placed on calorimeter measurements of disintes 
tion rate because they involve a minimum of un 
taintv. These measurements should be done o 
sample of high specific activity in conjunction \ 
abundance measurements. Further studies 
of this nature are under way in the Radioactiy 
Section at the National Bureau of Standards 


spect! 





Isotope 


We are indebted to the University of Kentucky a 
to T. I Davenport of the National Bureau of Stand 
ards for the loan of the equipment used 
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James B. Saunders 
A procedur rtesting large optical surfaces with relative nall optical andards 
. lescribed Simplified formulas are used to apply a statistical method { obtai increased 
ecisiot \ practical example is used to illustrate the procedure for testing surfaces that 
navy be a ! 1 to have revolution vi \ 
1. Introduction \ procedure is deseribed in which an ordinary 
sized standard test plate can be used to measure the 
Vith the increased use of massive optical parts shape of surfaces that are much larger than the stand 
z is need for a practi al method of testing lare ard Use is made of statistical methods for HMproy 
iLsurfaces with relatively small standards, such | ing the computed results. Simplified formulas will 
available in the average optical shop The con be developed that mav be applied by opticians, alter 
tional test of observing fringes between the a short training period, even though they may not 
ndard and unknown at various positions on thre be able to follow the mathematical derivations and 


we 
values either 


ial tests or new methods 
of 


Massive 
the 


1) methods testing 
tical for 


plates ol 


opt 


SIZCS 











enables the optician to form a 
of the nature of the surface 
extensions 


opty al 
al elements 
required for 





esting of Large Optical Surfaces With Small Test Plates 


rough esti 
For 
of the present conven 
Present 
hot 
as standard 
conventional 


nore ore 


must be 1 sed 


surfaces are 


methods are seldom available 
\ test for shape of very larg optical surfaces, by 








5 ins of interference of light, may be made with a 
(5 ndard test plate flat on spherical that is much 
: iller than the surface to be tested. The prinei 
involved are not new and may be considered 
mentary but those who perform the tests on 
) eal surfaces during the polishing are not as a 
sl ficiently versed in mathematics to apply 
BA he principles outlined in the present paper to the 
st advantage They can, however apply observed 
6, . 
tn to a set of simple formulas and compute th 
nsuing results 
H 
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theory The application of the general formulas to 
nN tvpical set of data will then be made to illustrate 
the testing of a surface by means of a smaller standard 
optical flat The method is not limited to plane 
surfaces. The sphericity of very large spherical o1 
aspherical surfaces may be measured with standard 
test plates ol approximately the same radius, if the 
viewing system permits the use of normal incident 
light. Figure 1, (A, B, C) shows simple optical 
arrangements for testing plane convex, and concave 
surtaces, respectively 


The basic principle hina also be applied to the 


testing of off-axis curved, or nonsvmmetrical 
surfaces, with modifications of the formulas How 
ever, as these are encountered only rarely the 


present papel will be limited to surfaces of revolution 

We will consider a Fizeau viewer (ig. 1) that has 
an aperture equal to or than the thre 
Spherical surfaces are referred to 
the spherical master with which they are tested, and 
plane or approximately plane surfaces 
to planes The and 
identical for plane and spherical 
discussion will deal with 
mately 


greatet a ol 


ure 
standard surface 
are referred 


equations menusurement are 


upPPronrt 


surfaces 
ot an 
standard 


the testing 


plane surface against a optical 


flat The unknown surface is assumed to be a 
figure of revolution about an axis normal to and 
passing through a known point on if In general 


this point will be at or near the center, if the surface 


is circular 


2. Experimental Procedures 
If the standard flat is adjusted normally to the 
collimated beam of light and the unknown is placed 
close to and approximately parallel to the standard 
interference fringes. may be with mono 
chromatic light of known Measure 
ments made on the fringe pattern permit a compu 
tation of the shape of the unknown relative to a 
plane (sphere for spherical surfaces over the visible 
region of interference. If adjustments permitted, 
the entire area of the unknown could be covered by 


obset \ ed 
W avelength 





moving the standard step by step without rotation 


im a plane parallel to the standard plane The 
avvregate of the resultant fringe patterns when 
properly assembled to form a composite pattern 


of the whole surface, would be quite similar to what 
would have been obtained with a standard surface 
that covered the whole of the unknown 

As it impossible to make the above 


assumed adjustments and 


is virtually 
corrections for rotation 
the 
standard 


changes im unknown, relative to 


the original surface of the 


separation of 
made 
applving the 
always 


must be 
changes 
required The 
subject to errors of observation 
of the standard 
are required to measure the 


by measuring these 


corrections corrections are 
If several successive 
the un 
shape along a 
Consequently 


The “Method 


and y ields 


positionings relative to 


known 


ciameter these errors accumulate 


it is desirable to reduce 
of Averages which is 


these errors 


easy to apply 


simple working formulas, is used since it provides 
adequate preersion 
In veneral large optical surfaces are vround ana 


polished by machines that produce figures of revolu 
about a point, usually located at the 
center of the surface If the 
revolution,’ its departure from a straight line that is 
tangent to it at the center of revolution is a measure 
of its departure from its tangent plane and conse- 
quently determines the shape of the entire surface 
other chosen plane The axis of 
2) is chosen as the intersection of a 


tion known 


surface is one of 


relative to any 


abscissas (see fig 
plane through the axis of rotation of the surface with 
the plane that is tangent to the surface at its center 
The unit of abscissas is chosen as the separation of 
along a 
diameter is 


reference marks chosen 
surface The 
The axis of ordinates is the 
axis of revolution of the surface. The unit of ordi- 
nates is one-half the wavelength of the light used 
The following steps are the chosen procedure for 
acquiring the data necessary to compute the shape 


spare ed 
the 


the axis of abscissas 


equally 


diameter of chosen 


of the surface. The master flat is placed on, and 
concentric with, the surface (positions A and A’ 
fig. 2 The magnitude of the air wedge between 


the two surfaces ts adjusted to produce a satisfactory 
number of fringes in the field (see fig. 3 The 
direction of this wedge is adjusted to make the fringes 
approximately normal to the line along which the 
chosen reference pots lie Except for plane sul 
faces this wedge varies from one reference point to 
the next. Consequently, if linear interpolation is 
used in estimating fractions of fringes at the reference 
the error in estimation varies with the frac- 
tion observed For most observers this error is a 
minimum when the fraction is 0.0 or 0.5, that is, 
when the reference point falls on the center of a dark 
or a bright fringe. Best accuracy is obtained by 
reading abscissas corresponding to the center of all 
fringes—both dark and bright-—and from these data, 
by nonlinear interpolation, compute the relative 


points, 


1. B. Scarborough, Numerical mathematical analysis, p. 446 (Johns Hopkit 
Press, Baltimore, Md., 10% 
?W. A. Granville, Elements of differential and integral calculus, p. 264 (Ginn 


& Co., New York, N. Y., 1911 
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represent the several position { the tandard test plat The correspor ne 
traight lines below represent sections through the standard surfa { referene 


orders: at the chosen reference points, The proce- 
dure will be explained in greater detail with the held 
of the fringe configuration shown in figure 3. 

The fringes (relative orders of interference) in 
position A are evaluated at the several equally spaced 


reference points by the method indicated above 
The standard is then moved to position B, which 
overlaps an appropriate amount of the area covered 


in its first position, A. Again the fringe readings at 
all reference points, covered by the standard in this 
new position, are evaluated. The standard is then 
moved to position ( and the corresponding fringe 
readings noted This procedure Is repeated until 
the surface covered by the master extends to the 
edge of the unknown surface of revolution. The 
positions represented by primed letters are a second 
independent set, useful in checking the precision of 
the method. 

In general, the fringe values at the reference points 
will not be integral. Consequently, the fractional! 
parts must be obtained by interpolation and, in ma! 
ginal cases, by extrapolation. Linear interpolation 
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xtrapolation is accurate only when the fringes are 
ily spaced along the straight line on which the 


rence points are chosen If they are equally 

meed, the surface is plane and no further test is 

essary. In general, the fringes will be curved and 

line nlinear interpolation should be used if highest 


iracy is to be obtained Of the several methods 
of performing nonlinear interpolation, the graphical 
method ts simplest and, for these purposes, adequate 
ce- i lhe details of this are best shown by a description ol 
eld This will be given in 
iter section, where the shape of the surface of a 
in iss disk will be computed from photographs of 


application to actual data 





ed terference fringes that are obtained when the glass 
ve disk is tested interferometrically against a standard 
ich | 

red 

bie 3. Derivation of Formulas 

11S 

en \ set of formulas are now derived for use in the 
wre bove-mentioned computation. In figure 4. the 
itil ved line represents the surface to be measured 
the chosen reference plane is represented by the 
‘he aALS of abscissas. The chosen reference points are 
nd cated by circles on the curved line The several 


of wsitions of the master plane relative to the coordi 
system are represented by straight lines at 





nts ious angles The coordinates of all chosen refer- 
nal points, which are marked on the surface of the 
al inknown, are (R, Y,). The observed fringe (a rela- 
on order of interference) at (R, Y,) is designated 
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fandard test 


plate, 


A, for the first position of the standard, B, for the 
for the third 
position, and so forth. The slopes and intersections 
of the straight lines with the axis of ordinates are 
represented by V, and Ap, respectively, for the A 
position, V/, and J}, for the B V7, and (, 
for the © position, and so forth 

The slope-intercept form of the equation for the 
straight line representing the standard in the .Y posi 
tion (Y=A, B, ¢ 


second position of the standard, CR 


position, 


, ete Is 


Vet (Xe—A VR 1) 
In position A, for example, the ordinates of the line 
are Ay )’; A VU,R, where Ap is the separa 
tion of the two optical The absolute 
values for Y, are unknown, but differences in Np», 
for any given position of the standard, are directly 
Consequently, the integral part of the 
smallest VY, will be subtracted from all \V,’s for pur 
computation. The quantity (V,p—Ny) is 
unaffected by this operation 

We will evaluate M/y and Yy by applying the 
method of averages (footnote 1) to the data. The 
grouping of observations will be made in such manner 
that there is no overlapping of the two groups of 
points except when an odd number of reference points 
are to be applied to an evaluation. When the num- 
ber of reference points in the two groups 1s odd, two 
observations are considered to be made at each of the 
points, and the two observations made at the central 


surfaces 


observ able 


OSC of 
I 
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’ i™ pats 
' ie ee OUT i 
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RADI 
hy ne 4 | j 
pravnnat anne divided between the two vt tip The iwo 
roups then have an equal number of observations 
and all observations are thereby assigned equal 
wemlht When the total number of reference pots 
iseven, there no overlapping of groups for an equal 


number of points to be allocated to each of the two 
rroups for equal we iwhting of datum 
footnote | 


In applving the method of averages 


we require that the algebraic sum of errors be zero 
for each of the two groups of abservations This 
requirement is represented by the following equa 
Thos 
4 . | 
zeta \ My SUR+ Ny NX } 


3. (Yat Na=MxD'R+iNxXo, 9 


where summations in the first of this pair of equa 


tions are from the lowest value of PR Ri) to the 
largest value of R(— FA, in this group and summa 
tions in the second equation are from the lowest 


value of R /?,) in the second group to the largest 
value of RU Rp im the second group, Equality of 
PR and Ry will result when the number of 
points is odd, whereas they will differ by unity when 
the total number of reference points is The 
number of reference points used in each evaluation 


is Ny 


Solving for 


reference 


even 


Vy ana \ 


from the pair of equations 


) we get 
ih PG OT ae 2 
_ ye 
Vx P } 3 
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Kquations (1 s), and (4 together with the 
assumption that the surface to be measured has 


revolution computation of ) 


for all chosen reference points and consequently th 


sVinmhetry, permit a 
shape of the surface 

The foregoing analysis permits a statistical evalua 
tion of AZ, and A all values of V Gi. e., A, / 
tf except the quantity «1 An in 4 
represents an error in the position of the referenc: 
plane and since this has no effect on the comput: 
shape of the surface, the observed A 
be free from error 


for 
erro! 


is assumed to 


As the surface is assumed to have revolution 
svinmetry, 
Ve=) 2 
From eq | Yp=M,R+(A Ag and = ) 


) 


M,R+(A,—A_e On eliminating and 


from these three equations, we have 


Vi, 


Ve=Ap—i(Ae+A_o 


As 1, does not enter into the computation, it need 


not be evaluated. 


4. Illustrative Example 


Equations (1), (3), (4), and (6) form a set of fund 
mental equations from which all values for }', m 
| be computed 
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» clarify any possible lack of understanding in 


development and application of these formulas 
vill apply them to a set of observations and will 
this 
2 represents the periph 


rule the resultant shape ofa surface In 
the large circle in figure 


! 


fa 3t.l-em-diameter disk that is to be measured 
vans of a standard whose diameter is 8.2 em 
the 


clisk represent the several successive positions 


even circles centered on one diameter of 
The eight circular fringe patterns 
n figure 3 (forming two mdependent 
ed and unprimed photo 
hing the fringe patterns formed by light reflected 
ally from the Lop surface of the large disk ana 


standard 





sets 


were obtaimed by 


tundard flat when in the several positions 
ited m= figure The differences in’ values 
ordinates (), - obtained from the two 
of photographs is a measure of the accuracy 
nable These differences nre due to errors of 
rvation and to an error in the choice of the 
point. If the unknown is ground and polished 
pindle that is not concentric with i the uXIS 
ation of the unknown surface will not be 
ed at the chosen origin of coordinates Hlow 


from the 
respec tively 


the average value of Vy, and )Y_, 
primed and unprimed 
free from the 


more pre cist 


ets of data 


be almost error of centering 
values for VV, 


Re ference marks 


order to obtain 
ollowing procedure is followed 
h then nssociated RP 


ed on the glass surface along the chosen reference 


values, or abseissas, are 
" These reference points are indicated by black 
centered in the small circles shown in figure 2 

tandard is placed in the desired position and the 

wedge adjusted so that a desirable number (5 to 
The present 
fringes approximately 
references line Considerable 
is usually tolerated 
I fringe pattern is photographed This is 

a vated for all desired of the standard 
\ fine rht line, using a needle point and a straight 


straig 
‘ is drawn the 





0) of fringes cross the reference line 
hes 
have the 
the 


hor prefers to 
thy 1) pendieular lo 
urture from this, however 

positions 
ueCTOSS photograph through 
reference points on if The positions of the 
the points marked by 
cking the photograph with a needle. This is to 
their positions, If no 
tortion is introduced in the photographie repro 
ction, the separation of all adjacent points will be 
qual. The intersections of all fringes, both dark 
‘ integral interference) and = light 
‘ ntegral plus half orders), with the reference line 
ire marked also by pricking the photograph. The 

ection of increasing interference is 

‘ iscertained from tests and notations made when the 
ives are photographed The relative orders of 
terference, beginning with 1, are indicated by 
nserting numbers on the dark fringes. In figure 3 
these numbers appear each of the fringe 
patterns. The lesser visible numbers, adjacent to 
the reference points, represent the 


rhe’ Thadhe 


nters of references are 


rmit accurate readimes on 


orders of 


orders of 


“ol I 
eed ' 


abo e 


abscissas (or 
Ml — 


nicrometer eyepiece may be used to read the data directly and thus d 


vith photography 













O17 4G 





P-values) of the several reference points that were 


marked on the unknown surface 
In general, the seale of the photographic reproduc 

tion will differ from that of the chosen abscissa scale 
For interpolation purposes any convenient seale may 
be used to measure the relative positions of potnts 
on the photograph \ mill 
meter marks was used for this compilation We will 
define the seale used as the P-seale and the readings 
from itas P-values. The /?-values corresponding to 
all reference potnts anal fringes are read off the centi 
which 
The corresponding 
numbers associated with the fiduciary ref 


centimeter seale with 


meter scale and inserted in column | of table 1, 
represents a tvypteal data sheet 
R values 
erenee pots ana appearing in the photographs 
inserted in column 2: the le 
lo integral 


ure 
values ‘ orresponding 


and half-integral orders of interference 
are inserted in their respective places in column 3 
If an observed 7?-value, corresponding to an P-value 
falls pre isely on oa fringe reference 


two identical 7?-values are 


pont centered 
ona fringe inserted as for 
PR 2 and riven here 

The object of this listing and the following compu 
1y-Values 


ly 5 mm the case 
tutions i to obtam more precis orders of 
interference) than can be obtained by direct interpo 


lation or reading from the photographs 


A direct plot of the 7?-values versus sA,-values 
first entered m column 38, table 1 (the blank spaces 
to be fille d later shown as small Dylan k dots in hivure 


>, may be made A smooth curve is then drawn 


through these points The values for A, in column 
3, Which correspond to P-values in column 2, are 
located on the curve The corresponding ordimates 
represent the interpolated 1, values desired The 
Values with asterisks, shown in table 1, are more pre 


cise Values obtained by an interpolation process de- 
scribed inp the following 

The values for A, from the 
direct plot ure nol as precise ws the cine tly observed 
values from which the plot was made 


paragraphs 


determined above 


Nlore pre Ist 
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P (CENTIME TERS) 
nonlinea 


Direet 


inte polation of fringe orders 


Fiovune 5 and residual plot of data fo 


values may be obtained from a residual plot.! To 
obtain data for the residual plot, a column of values 
(table 1, column 4), computed with a slide rule, is 
added to the table of observed quantities. This 
column is a product of corresponding /’-values and 
an appropriately chosen constant. This constant 
represents the slope of a straight line that crudely 
approximates the plotted points. It is indicated in 
figure 5. The approximate slope of this line is 0.7. 
A relatively large error in the choice of the slope 
contributes little or no error in the final results. 
Values of O.7P are accordingly computed to two 
decimal places. The sums of corresponding values 
from columns 3 and 4, except where vacancies exist 
in column 3, are entered in column 5. These values 
are plotted, on an appropriately chosen scale, against 
corresponding /?-values, and a smooth curve that 
best fits the points is drawn in by inspection (residual 
plot, fig. 5). One may now locate the points, shown 
as circles, on this curve, whose abscissas are those 
observed for the chosen reference points. The cor- 
responding ordinates are read off and inserted in the 
(0.7P+-A,) column. These values are marked with 
an asterisk for distinction. The corresponding 
values of 0.7P are subtracted from (A,+0.7P), and 
the resultant Ag values, shown also with asterisks, 
are inserted in the vacant places of the A, column. 
These are the desired orders of interference at the 
several reference points. Values for Vy(.V=—B,C.D, 

) are obtained from all photographs in this 
manner, 


‘H, M. Goodwin, Elements of the precision of measurements and graphical 
methods, p. @ (MeCiraw-TTlill Book Co., Now York, N. Y., 1920 
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Computation sheet for statistical evaluation of 


: 9 

PaBLe 2. 

The average of corresponding values in columns 3 and 
receptable result 


represents t 


/ ix Vy / bz 
‘ &. 23 
TD 
“) 
! 4.42 
4. 33 0. 00 8.71 Om 
2. 41 ow 7.13 on 
1. a2 2 7h 2 
Ow “> a ¥ 
4 uy OS 3 Lo 
1, 57 | $4 1. 58 
‘ 2) 1. 4 y .| 
ie 0.77 4.13 
s a n~ ia 
" 67 » 
( 0 6.74 
] 54 5S 
2 10. 30 aud 
13 12. 37 11. Se 
4 14.51 13. 0% 
l 16, Se 2 


In order to determine the shape of the surfac 
from the fringe readings and their corresponding 
R-values, these quantities are now 
applied to the formulas of eq (3), (4), and (6). Table 
2 is a typical computation sheet. From eq (6) on 
computes ),=3.33— (8.23 4-0.39 0.98, and like- 
wise for all values of )', in the range covered by th 
standard in ; 


abscissas, or 


position A (column 3 of table 2 In 
order to proceed to the evaluation of V, for other 
values of FR, eq (3) and (4) are used Values for Be 
found in the same manner as described previously 
for Ag, are inserted in column 4, table 2 This per- 
mits a summation to be made of all terms shown in 
eq (3) and (4) and consequently an evaluation of 
M, and By. As the number of reference points 
covered by the standard in each of the two positions 
A and #, is odd (i. e., from R=0 through R=4, o1 
5 points), they are each assumed to represent 2 
observations, making the total number of observa- 
tions 10, which may now be divided into 2 equal 
groups. With this in view, it will be seen that 
. 
>  R=2X0+2X1+1X2=4, for all 
7 
summations shown in eq (3) and (4 The resultant 
values for M, and Fy are, respectively, — 1.580 and 
8.67. Shlide-rule computation is adequate except 
for final precision. In general, the computed para- 
meter B, will differ slightly from the observed value 
found in table 2. On substituting these values for 
My and Xy in eq (1), and using values for R and Py, 
in columns | and 4, values for )'y may be computed 
from R=5 through 8. These operations are now 
repeated, using consecutively the C and D data, to 
compute additional values for VY, to the edge of the 
surface. These values are inserted in column 3, 
table 2. A similar treatment of the primed data 
yields an independent set of value for Y_, shown in 
column 5, table 2, which should agree approximately 
with the above set. An average of Vp and Y_, 
values for each value of R is the accepted value 
This averaging operation tends to decrease the effects 
of error in choice of the axis of rotation of the un- 
known surface as well as observational errors. 


and = likewise 


Wasuineron, December 3, 1953. 
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Bacher and Goudsmit Theory of Complex Spectra 
R. E. Trees 


The Bacher and Goudsmit theory expresses the energy of a state of an atom in term 


of a linear combination of energies of states of other ions of the same atom This theory 
is tested in the spectra with d- and electrons and general conclusions drawn about its 
accuracy 

In deriving their theory, Bacher and Goudsmit formally introduced many-partick 
nteractions Chey showed (to a limited extent) that this was equivalent to a second 
procedure in which the zero ord r wave functions were assumed to be inde pe ndent of the 
ionization and successive orders of perturbations were included 

Recently has been shown that a linear theory, which is a combination of the Slater 
theory and an empirically derived correction, predicts the separations of terms in any one 


configuration with about the same accuracy 


as the Bacher and Goudsmit theory: in thi 


theory, the wave functions are allowed to depend on the degree of ionization It is show! 
that the Bacher and Goudsmit theory can be interpreted in a third way that is consistent 
with this linear theory 

In the light of this latter interpretation, Racah has suggested an empirical modificatior 
of the Bacher and Goudsmit theory that often leads to better results: this sugwestion has 


heen it corpora ed in the present paper 


1. Introduction 


Bacher and Goudsmit have shown that it is pos- 
bie to express the energy levels of an atom cor- 
ectly to second order by taking linear combinations 
of the energy levels of the ions of that atom.! They 
ave explicit formulas for these linear combinations 
for atoms containing s- and Pp electrons, and tested 
the theory in the spectra of oxygen, nitrogen, and 
The theory usually predicted the absolute 
erm values with an error less than 2500 K 
rations of terms in a given configuration were often 

ven very much better, sometimes within a few 
hundred kaysers. 

In section 2 the theory of Bacher and Goudsmit is 
derived briefly from a new viewpoint. Bacher and 
Goudsmit derived their formulas by formally intro- 
ducing many-body interactions; no explicit form was 
riven for these interactions. They were able to 
show that in the first and second approximations this 
was equivalent to using zero-order wave functions 
that were independent of the degree of ionization, 
and, with these, making first- and second-order per- 
turbation calculations of the energy, respectively. 
In the present paper, a third derivation is given, in 
which the zero-order wave functions of the Slater 
theory are allowed to vary regularly with the degree 
of ionization; an explicit two-body interaction is also 
introduced. This interpretation is justified by the 
observation that there is a smooth variation of the 
one-particle and two-particle parameters in a linear 
theory. The more specific and easily understood 
meaning of this interpretation in terms of a linear 


carbon 
Sepa- 


R. F. Bacher and 8, Goudsmit, Phys. Rev. 46, 948 (1934); referred to in the 
resent paper as “BG 

*Kayser”’ is the name adopted by the Joint Commission for Spectroscopy for 

mit of wave number, cm~', and is abbreviated “K 


R. E. Trees and M. M. Harvey, J. Research N BS 49, 397 (1952) RP2378. In 
tion | of R P2378 the literature is reviewed, and the methods of analyzing spec- 
s with the Slater theory, for including effects of configuration interaction, and 
naking allowance for polarization energy with the L(L+1) correction are also 
itlined and applied. In referring to a “linear theory’’, the L(L+1) correction 
combined with the Slater theory, and the radial integrals are regarded as 
rameters 
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theory is given in section 5; it is equivalent to assum- 
ing that the parameters in different ions of the same 
atom are given by series expansions in powers of the 
degree of ionization, It is shown that the BG nth 
approximation accounts for terms up to and including 
the nth power for the one-particle parameters and 
terms to the (n—1)th power for the two-particle 
parameters 

In section 3 a few formulas applicable to terms in 
d", d"s, and d"s* configurations are obtained by using 
tables of coefficients of fractional parentage and 
recoupling transformations given by Racah. In 
section 4 an application of these formulas is made in 
the spectra of titanium, vanadium, chromium, and 
zirconium, It is found that separations of terms in 
the same configuration are given very closely (the 
error is less than 400 K Separations of terms in 
different configurations may be given very poorly, 
however, as an error of 6000 K is found in titanium 
(Z=22). The latter type of error is found to 
decrease very rapidly in spectra of more highly 
ionized atoms: in chromium spectra (Z=24) the 
error is only 600 K, one-tenth the error found for the 
related titanium spectra. The error is also found to 
decrease for spectra with higher nuclear charge as it 
is only 2000 K in zirconium (Z=40). 

Racah has pointed out that a power-series expan- 
sion, of the type described previously, converges 
more rapidly if spectra are chosen having the same 
degree of ionization. A comparison with experiment, 
using this empirical modification of the original 
theory of Bacher and Goudsmit, is also given in 
section 4. 

As indicated previously, the Bacher and Goudsmit 
theory can be interpreted in terms of a linear theory 
in which the radial parameters are given by series 
expansions in the degree of ionization. Checking 
these expansions is the easiest way to check all the 
data against the Bacher and Goudsmit theory. The 


‘GG. Racah, Phys. Rev. 62, 438 (1942); 63, 367 (1943), These papers are referred 
to herein as LI and III, respectively 
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conclusions ren hed from the direct apply ation of thre 
theorv already cited above have been vertfied by 
using all available data with this method The de 
tailed caleulation nemled for thi vertieation are 
omitted for bore vily 

In section 6 the effeet of weak configuration inter 
netion mn which cone electron is exetted 1 considered 
and om ection 7 the effect of irrationality of the 
energy cigenvalues is analyzed, The purpose is to 
ilustrate in more detail the direct appheation of the 
theory of BG and to indicate where a direct applica 
tion might lead to results that differ from those 


obtainable with a linear theory Such differences are 
probably tow mall to be detected withoun the necurney 
obtaimable at present with any theory but more com 


plets compart ons of theory with the presently avail 


able experimental data are needed to verify thi 
conelusion 
2. Bacher and Goudsmit Theory 
2.1. General Theory 
ln this seetion is riven a new derivation of the 
theory in which it is assumed that the interaction 


parameters vary smoothly, The Hamiltonian of 
the n-electron system is taken as 

H-MWeAH (1) 
where 

of 1 “Ae 
H=-> ( —— p ) 

and 

H,= >> - 


/7, and //, are examples of what are called “type 
one” and “type-two” respectively Hi, isa 
sum of funetions, each depending on the coordinates 
of a single electron, and //, is a sum of functions for 
two electrons A particular selection 7? of r of the 
n electrons will be indicated by 2? The sum of the 
Hamiltonians //(7’,) of the r-electron systems for all 
possible choices 7? is a symmetric function of the 
and the n-eleetron Hamiltonian /7/ is 
expressible identically as a linear combination of such 
two or 


terms 


ne lectrons. 


vymimetric sums for more values of 7 


H= SOSH P): (2) 
D 


The conditions that this be an identity are the follow 


Se W)C) and VEC" YK S)=(G) @ 


In general, (") is the binomial coefficient when | 





The values of tl} 
indeterminate of 


qr ad 4 


for bo O 
eflicients ¢ 
than two nonzero ¢ 


und zero 
there nre 
The method ol ¢ 


sil ly cusses will hy 


will be 
values 
ing best values of the CC, in 
sidered ly | mw 
The eigenfunction wy) of the n-electron syst, 

a defini 

ane P Is Spree ihe | by ; ol 7 
/-values; the electrons are regarded as distinguish 
even though some may be equivalent The elect ns 


in the seler tion ,? have allowed ementunetions y 


assumed to have in its specification 


fiyvuration/,, / 4 





of the r-eleetron Hamiltonians //(7 the residual 
n r) electrons have allowed eigenfunctions ¥ of 
the " ; electron Hamiltonian The /-value of 


the electrons are included im the Spec ifications i 
and y, along with S- and L-values for the combine: 
system and any other quantum numbers needed to 
specify the coupling The funetion formed by ve 
tor coupling electron eigenfunction to 
the r-electron eigenfunction is indicated by ¥(a,p 

\ linear combination of thes> functions can be formed 
that will be identical with w(y); this was first pomter 


the 7] 7 


out by Bacher and Goudsmit, 
v\4 S35 (ap,7 h ¥(a,B,y). | 
af 
a,8,y hy) is called the coefficient of fractional parent 
age The caleulation of the coefficients of fractions 


parentage is a mathematical problem that has been 
considered in detail by Racah He has tabulated 
these coefficients for configurations of equivalent p 
and d-eleetrons in III; their use will be illustrated in 
section 3 

A partial energy /(a@) of the choice 7’, is introduced 
that is associated with the r-electrons while they ar 
part of the n-electron system and combined to form 
the state a This energy cannot be observed, and 
must be approximated from the observed energy of 
the r-electron ion of the atom in the state a, which is 
defined as WW (a The average cnergy of the choice 
P, when present in the original n-electron configura 
tion y is determined with the help of (4 

Eiy:P.) > (aB.vh Kia ° 0 

op 


The corresponding approximate energy can be cal 
culated from the observed energies in the r-eleetron 
On, 

V(y;P,) > (aB,y}y) Wie). (6 
op 


The energy of the n-electron system is then given 
” 


identically by an expression similar to (2 


E(y)= SoC SSE (y;P). 7 
D j 


’ Ther more than one tert 1 given (SL)-value i 
configuration considered (the configuration may be that contained in a, 8 or 
In such cases the function obtained by a legitimate coupling procedure will t 
be an eigenfuncticn of the Hamiltonian, even though it bas an allowed symmet 
f it mixes with other allowed eigenfunctions in that configuration having 
(SL)-valu In such cases a purely formal difficulty arises becarms 
and the theery given by Bac! 
amplified by more or less arbitrary additions, In 
hese omit 


may be sllowed wit! 


wne 


characterization by parentage is indeterminate, 
and Goudsmit must be 
tions 2, 3, 
cases are considered in 


snd 4 of the present paper these cases are omitted 
soctions 5 and 7 














































ponding approrimat 3 caleulated 


wm the approximate average 


Wiy=MOST 


The energies VU and WI od 
1} used in Bay for both expres 
has Wia), Wla,B), ete 
i4) in BO’ are 
terms, but on 


above are the nine as 


21008 


(Quantities 
on the meht side of eq r 





sometimes enermes of observed 


, ; eetral reneral they are linear com 


obs rved 


mtions of the eneren in the present 


stron they would bye r placed 1 W 5 r | y r’, 
indicated by W 


hand side of 5) and (26) of BG 


nD Direetly observed eneret nre 
n the right 


Nhe difference between the exact value (7) and the 
proximate value (S) is estimated by separating 
Ie’s and W into their type-one and type-two 
and by assuming that the following expar rons 
n be mac 
Thies 
hen Wile / vi 1 aly) ! 
) 
Wola Efa)) 1 h } h ? 
| 
The coefficient a;(y) in (9) will be eliminated from 
prt ; the difference between ual (S) af 
na , 
mel ‘ / 
ited : 24! , , (— )( | ) 0 10) 


climinated if 


)(,)=0 Ul 


values defines an 


the co flicient h, 2’ will be 


al 7 } ( 


unre | 





how i ol f 


- 


upproximation of 


y Ol the same order as the largest value of ¢ for which 
his 0) is satisfied (it is assumed that (10) and (11) are 
mee atisfied for all smaller values of ¢ and that relations 
ira are satisfied For mstance, na second approxi 
mation, the type-one energies are expanded im a 

enes contlamning powers up to and meluding the 

5 squares of (7 ? whereas the Lvpe-two energies are 
expanded only to first power in (n—? Formulas 


al 3), (10), and (11) are not independent. They may 
by combined into a definition that is easily shown to 
be equivalent to that just given, by 
coeflicients C, to satisfy the 


tions in order to define 


requiring the 
following (m 1) equa 


an mth approximation 
~ S°C,(0 ry )() 5 t33 t=0,1 ,m) 


In the derivation of BG, many-body 
are formally introduced, According to the present 
notation, this Is equivalent lo considering hot only 
type-one and type-two terms as being present in (1 
but also type three, type-four, et terms This 
leads immediately, by generalizing (3), to the defini 
- approximation given by Bacher and 


Interactions 


tion of the m' 





the following expansion ts utilized for powers (0 of 
an integer p:’ 

p= Sa. (O(") 14) 
where 
Insertin 14) in (12) leads to the result 
Sa ( No! ‘0 (= |] m) 
os « JS 

Iva) 

Equation (13) is equivalent to the formula 
> (’ 4 I k=1,2,. " I 3a) 
| } h 





Croudsmut Accordmg to thei definition thre ( 
values satisfy the following (n 1) equations 

Ss ) . 

DC) } )-( ' ) j | 13) 
To show the definition 12) and (13) are equivalent 


It follows therefore 
binominal coeflicient relation 


by repeated application of the 


Ia) 
satisfied 


that the left sil of 


whenever (13) is 


and the use of Loe 
vanishes, as if should 
When ¢--0 in (12 
when ke 1 
defining equations is demonstrated 


the equation is the same as (13) 
so the equivalence of the two sets of 
\ special Mi 


approXimatllon is defined hy the following set of ¢ 


th 


values 
} vii ] 
( ( ) | ] 4 Mn | } 
” 7 | 
( 0 } m 7 / l j 
15) 
If there are only two nonzero values of (C., and 
these are C, and ¢ then they are defined by (15 


and relation (S) has the expli it form 


SW y;P)—(n— DSW 


} 4 


VW y yp 1) 


Relations (8) and (12) of BG are spe ial cases of 
16), for which » has the values 4 and 5 respectively 
alternative, it can be 


s an assumed that C) and 
( are the nonzero coefficients, and then the fol 
\ ‘ Gi, Me iby F. J, 
1 A { xs 


lowing formula is obtained by solving (12) or (13 


the specification y will be omitted for brevity 


S*W(P,). (17) 


SoWiP > > 


——/ 


2.2. The Second Approximation 


If ¢ ( and © are taen as the nonzero coefli- 
cients of the second approximation the following 
formula is obtained from (15 

Wn) SMP) -a—- 29S WP, 4 
! , 
,( 20 SOW IS 
Relations (9) and (13) of BG are special cases of this 


which » has the values 4 and 5 
perturbation theory 


formula for respec 
tively In appendix lL of BG 
was used to show that IS) ise juivalent to a ealeu 
lation that is correct up to and including the second 
reference to “the theory of 
in the present work applies 


ordet In reneral 
Bac her and CGioudsmit”’ 
chiefly to this second approximation 

The complete set of one-electron functions used 
in appendix | of BG) was that appropriate to the 
ionized state of the atom with only a single electron 
shell The Ssatne zero-ordet 
functions were used for the atom in all other states 
and the second-order effeet was evalu 
This Is & poor zero order upproxti 
mation tor spectra,’ so that 
effect.’ in the sense used by Bacher and Goudsmit, 
is not synonymous with polarization energy The 
latter energy is usually defined as that part of the 
observed energy that cannot be explained by a cal 
the one-electron functions are 
Hartree-Fock method for each 
No investigation has been made 


outside of a closed 
of ronization 
ated with them 


sone second order 


culation in) which 
determined by the 
atate of mnization 


to determine the variation in accuracy of a second- 


order calculation with change of the zero-order 
eigenfunctions, but Bacher and Goudsmit often 
obtained poor values in their calculations because 


the higher-order effects became more important as 
the functions worse approxi- 
thations 

The influence of polarization on the separation of 
terms in a single configuration containing d- and 
selectrons may be represented by a correction pro- 
portional to L(+1)." Raeah has shown that this 
correction may be derived by assuming a special 
type of linearity for the polarization energy, and he 
has outlined a proof of this assumption of linearity 
The proof is based on the assumption that in the 
configuration d@" the L(14-1) correction arises from 


zero-ordet became 


10 to ©) perceat in sticeessive ions 
when they are evaluated by least squares; they would have the same value if the 
rero-order functions were the same in all ions A similar problem is considered 

See also FE. P. Wigner, Phys. Rev. 94, 77 


Similarly defined parameters change by 


by A. Rahman, Physica 19, 377 (108 
(ia 
*R. BE. Trees, Phys. Rev. 84, 1080 (1051) 


*cG). Racah, Phys. Rev. 85, 381, (1942), and private communications, 





interactions with configurations d" 


*P and d®-2Ul’ + vay 
electrons of the original configur; nN 
The linear behavior then follows from 54 
Alternatively, it can be demonstrate: 
comparing (33a) of ITL with (33c) of LLL, and as 
ing that the perturbing configurations are = s: {] 
ciently far away so that the energy denominator. jy 
the second-order perturbation can be assumed 
stant The proof shows that the (14-1) corres n 
is accounted for initially in the first approximation 
so that it is on the same footing as the Slater em 
Ap interaction in which only one electron is excited 
is included fully for the first time in the second op 
proximation, as shown by (53) and (59) of BG 
Any deviation from linear behavior that is produced 
assumed negligible 


have two 
excited 
of BG 





by such an interaction must be 
if the L«/ | 
polarization 
In section 4 it is shown that the Bacher and Goud- 
smut with the ob 
served separation for the d® s 
figuration in both titanium and zirconium, the errors 
being 357 K and 2938 K, respectively Bacher and 
Goudsmit noted that the separation of terms in t 
configuration well predicted in 
containing 2s- and 2p-electrons It 
that the error of the Bacher and Goudsmit theory in 
predicting separations of terms in the same con 
figuration will generally be less than 400 K and that 
is therefore about the same as that ob 


correction is used to account or 


theorv leads to close avreement 


two terms mn 


con 


e 


she are spectra 


Is conclude ad 


the accuracy 
tained with a linear theory 

A large part of the polarization energy is the same 
for all terms of a given configuration. The energy 
common to all terms can be investigated with self 


consistent field and other variational methods, but 
itis not easy to get a reliable pieture of the behavior 
from the few calculations that have been mad 
One approximation, used fairly often with minor 


variations, ascribes an energy of 5000 to 10000 K 
either to each electron or to each pair of equivalent 
electrons; interactions between pairs of nonequiva 
lent electrons are neglected in the latter approxima 
tion. The energy is independent of Z but depends 
a little on the types of electrons that are paired 
This part of the polarization energy is very large 
and the departures from the described behavior 
could be appreciable 

The Bacher and Goudsmit theory can be tested in 
this respect by calculating the separations of terms 
in different configurations. In section 4 it is shown 





s noted in BG that 


X VILLI of BG, ecaleulation It 
wreement better than 200 to 300 K must be regarded as fortuitous so tha 
ultimate significant accuracy is regarded as about the same as obtained with the 

L+1) correction Many separations of terms in the same configuration ar: 
predicted with errors greater than 1000 K in BG, and this must be attributed 
the neglect of configuration interaction, or to the fact that the theory is less a 
curate for spectra with 2s- and 2p-electrons 

F. Seitz, Modern theory of solids, p. 2477 (McGraw-Hill Book Co., New 

York, N. Y., 140 Seitz evaluates energy per electron rather than per pair 
equivalent electrons, but his remarks on page 248 justify the latter interpretatic 
Self-consistent field calculations of O, O*, and O**, including exchange and « 
flwuration interaction (PD. R. Hartree, W. Hartree, B. Swirles, Phil. Trans. Re 
Soc, [A] 238, 220 (1930)), support an interpretation based on pairs, a value of 7000 
K being associated with each pair of equivalent 2p-electrons. Although a stro 
dependence on multiplicity is eapested (page 2%), this is shown only to a mir 
extent in the oxygen calculations; if present, it would be accounted for in 
previous consideration of the separations of terms in the same configuration. A 
additional peint that should be noted is that comparison of the self-consiste 
field calculations in oxygen (noted above), with and withcut configuration ints 
action, indicates that configuration interaction is important 


lables XVIT and 
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hat there is an error of 6000 K in predicting the 
, nergy difference between a term in the 3d* con 
‘ uration and one in the 3¢° 48 configuration of tita 
ium, and a corresponding error of 2000 K in similar 
reonium spectrum with 4d- and 5s-electrons, and 
o that this error decreases rapidly as ionization 
creases mm the spectra with 3d- and ts-clectrons 
ause the spectra of titantum, vanadium, and 
hromium show errors of 6000, 1800, and 600 K 
pectively Bacher and Goudsmit found an error 
: about 2OOO K in the spectra with 2s- and 2p 
: ectrons: the decrease of the error with increase ot 
: was slight in these spectra 
sox Che spectra considered here contain four electrons 
a mutsicl of closed shells and ure the easiest spec tra to 
nalvze with relation (18 \ more complete ex 
loration of the theory would consider also spectra 
th five or more electrons 
2.3. Higher Approximations 
nel A third approximation may be obtained from (15 
1e n which ©,, (2, Cy, and Cy are nonzero; 
ira 
yo Wyn= SW (P)y—n—4) SS WP) 
ns a WP 
‘ ) 
= “SW 19 
I 
ut 
a Relation (14) of BG is a special case of this formula 
™ for which n=—5 Bacher and Goudsmit considered 
= this approximation equivalent to a calculation by 
K perturbation theory correct to thrd order, but they 
nt vere not able to carry out the complex calculation 
PM that would be necessary to demonstrate this rigor 
a ously. Separations of terms in spectra of oxygen 


little better in the 
In section 5 the 


ds ind nitrogen were calculated ATT 
HW, than in the W, approximation 











" relative accuracy of Wy, and W, is discussed from 

a mother viewpoint (formulas (56) and (60)), and it ts 

concluded that higher approximations can only ac 

in count for a small part of the erro already noted in 
na the 1 approximation 

vn Any interpretation used to justify the method of 

Bacher and Goudsmit leads to the conclusion that 

is more accurate the more similar the binding of 

the eleetrons (and the structure of the configurations 

for the different energy levels utilized in the formulas 

XIV ni I ! 1 hh 1 « 1 
. - read T 7 bw a] ' 1 ul 1 

{ti th nw en wit 

by W I by t t ' that tl 

F i n il im tat XVIII K " ‘ val number 

‘ ent was worse 1 with W. for 

nfigur m xy w tlLan 1000 K 

X1\ BG), whereas tl Wy apt x t indicated term values in the 

t A figuration with a larger error, about 1400 K (table XV of BG Bacher 

4 Cioudsmit note that the 2s? 2p* configuration calculated better in the UW 

simation, regarding the sclectrons as | the closed shell, than in the 

prox ’ n The relative err fal lute term values will naturally 

‘or th er approximation as the absolute value of the energy to the 





For example, in calculating term values for /" all 
configurations specifying a choice of r-elections would 
belong to 7’. On the other hand, in calculating term 
values for /"~' s the selection of r-elections is specified 
by either the /’ or the /’ s configuration, so that 
there are twice as many different configurations to 
According to this argument, a given 
order of approximation would probably predict terms 
in the 7" configuration best, terms in /"~' s next best, 
and terms in /"~? s? least accurately This argument 
has some confirmation which was pointed out) by 
Bacher and Goudsmit."' Bacher and Goudsmit also 
have noted that the method is inadequate for con 
figurations such as 2s*2p*4s, in which the total energy 
of the excited electron is less than the 
of the inne such 
separations of terms in the same configurations would 
Similiarly, the method gives 
negative ions.’ 


be consider ed 


electrostatic 


energy ones In even the 


CAUSES, 


not be given accurately 
poor predictions for 


3. Formulas for Spectra With d- 
and s-Electrons 


3.1. The 1" Configuration 


When the configuration specified in ¥ is /' 
6) has the form 


, equation 


WitySL;/ + > l(a SL, l"~"(a2S2L2) 


SLUM"ySLYWlaS Lh 20) 


If this relation is used in (17) for Wi . the 


. yi fi ~ 
over P introduces a factor ( ) ( ) n, and (17) 
/ /i 


Stim 


be« Omes 


Wi ("yS1 “—— DS) (aS, Ly ISL ySL 
‘Was, L)—-—s WO. (21) 
“ - 


The first term on the right defines the separation of 
terms in the / 
as a relation viven by Racah for calculating niatrix 
elements of 
tween pairs of 
bv utilizing (20 


configuration: it has the same form 


it reneral scalar operator ty defined be- 


interacting electrons Similiarly, 


with 7 ” relation (16) becomes 


n(n | . 
Wi'ySL >» POS,L,)1 oS ol 
4 Lia I 
SLV"ySLY WPS n(n YW 29 
‘Table XIV of BG Absolute ter i { the 2 nfiguratior yeen 
are predicted within a few hundred Kays« t t iv it tt t t fleant 
wreement obtainable Errors for 2 sare about 1000 K, and they are about 
nH) WK for 2s? 2p 
In footnote 7 of BG it is noted that tthe siculated electron aff en 
e000 K and that the observed valu SOOO Ak Phe discontinuity bir ng 
energy of negative ions has been discussed by Ta-You Wu in Ph Kev. 89, 629 


*G. Racah, Phys. Rev 








riven 


clements of a sealar 


ame as an alternati formula 


This 1 thre 
by Racah for calculating 
operator 

lt is noted in Ill that the formula corr 
to (22 very long ¢ 
time only ith 1 few }>: 
direct use of "dy impractical except 


mais 


ponding 


need Hleulations and is of practs il 
riicular cause Althoueh the 


when | 


ori the trivial case where ros unit a step-by-step 


rein with | rad cle 


conerally practy il \\ 


procedure tn i 
ith thi 


mount at om tine ! 


procedure onl table of the coetherent 

SLILSLYV \/ are needed These coeflicient 
ire defined b O) TEL and (11) TEL, and the coetherents 
have been tabulated in TEL for p rie ceontieu 
ration 


approXimation IS 
1) term by usin 
tuble 


\ ne example thre econd 


evaluated mm detail for ad 


tep procedure By usin 


will be 
”) mon tepebyy 
Ill of TIL, relation (20 


hoe cvautiie 


i 1) Sion I’ 7 1On 


The on the meht are formally evalu 


table Ll of II 


ewe ee 7 on 


obtained directly 


Alternately “) could have beer 

from 1) with t > by first enleulatin the coeth 
erent of fractional parentage \,/ SN. / 
1) }) with 2 ol itl Phe tep-by tep 


procedure used above | imopler because the quantum 
number a \,/ ado thot 
ith, Ihe 4 


mite CLS 


( ) to obtain the desired formula 


enter thre Coicuintion expire 


valu rr! ty nr now in erted 


rl the un over 7’ rep wed by factors 


Wid? D) = GOW ete Inout I 


vood one to use in testing the 


All term Values used 


Formula (27) is a 
Bacher and Goudsmit theory, 
in the formula should be unperturbed 


known experimentally Phe 


thre youre 
following 


and 
likely to be 
formulas, obtained in ao simular 
generally useful: the that are 
perturbed by configuration interaction 


two 


not so terms most 


likely to be 





are starred 


f 


: ' 
Formula th all other tern would require an 
tion to the theory olf Bacher anid Croudsmiut L hie 

f duplieated terms wm and o; 


paurentuges of 
on the particular 
herens mm the formulas for the te 


defined by p 


Waottiad depend pectrum 


eon “lieration 
iven above the coellerents were 
} 


mathematical consideration 





3.2. The i '? © 


hguration 


when thre 


i ‘g 


ruration, ¢ 


evaluation of peciiieation 


contains carried out 


the /"°'’ conti in be 


imularly if tables of coeflicients of fractional parent 
avatlable.* In 


nee tor nonequivalent cleetronus are 
tables are not available and the co 


be calculated from tables for 


veneral such 
eflicients 
lent 
between different coupling 
"2 of Il Again a step-by-step procedure is simplest 
so that only selections 7? with | electrons need be 
considered, The parent term in the 
assumed, for simplicity, to be in an 
The 


electrons has an 


Thtuist equiva 
electrons with the help of transformations 
Ss hemes given in section 


conhgura 


tion again is 
exactly 


Laing 


! 
specified state single election con 


the equivalent average 
WSLUSN L'il which is the same as 
the ron We/ yS/ The 


specified by the 


energ\ wil 
the 
only 
[*-2] 


parentage 


eneryey ith 
new the 
configuration The 
will have the form 


obser \ ed 


feature is selection 


coetlicients of fractional 


l’ aS hLyUvS°L’) ISL’ YLUS'L’ 


Y2(aS,L,ISLW" ySI 


STAT) i 8’ 18. 10 (8° 2!" is’ 1’ 0) 


On the right, the first factor is the coefficient of 
fractional parentage for equivalent electrons; the 
lables for « i ha } ven by Mi ko hy Rev 
1, 871 (198s), With? , the normalization of the coefficients, the 
! ver 7’ yu | t t ike t NY it iced by factors 
( ” ) + wer Pa t ’ ! ! calculatio 
the tabl ! vllow wha to be the etd t Racal's tables, t 
elicient wtional ‘ ace have been nermalized to unity for each selec 
t ithe | t paper 
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cond factor is a transtor tion coethicirent evaluated 
») of III 
yy, / \/ LIS . af YF 
S+1)(28 / >] 


W(Sh 48” S'S, ) WOLVU Ag 


e W-funetions in this expr on are 
i i) of 1] ome Vulue ot the \\ funetion | 
» Deen tabulated 


Che particular case where in s-cleetron ! 
one we consider in detail In this ea 

ilt of ipplyving a ral CcLUce to a tart 

ple formulation ‘To evialtinte thre nveruce 

rvies tor a term S/ S’/ 

a) Evaluate Wa N/ nm term of the 

erved energvu iI SN, / i already ae ermbed 

b) If the spin S; of the ton term in does not 
tive nme value SS a he oral terme wut ley 


kderation mn then add an electron to th 


term. maki the re tant pin S equal to NS 
pin of thre parent erm in = chiar 
| uw, / to WN, | \/ 
If thee pin S; has the value S then add 

ectron to the ton, an pit tt latter into two 

' each bavi one of the two spins that are 
owed. 1. ¢ ive S the twe ius Sand (2S ‘ 
The part havin the resultant pin SI iven a 
oelherent 2S | whereas the other part has a 
eHherent of umitv minu th ‘ change 

N, / into 2S | iI uN’ ] \/ | 
S i uN] S S)/ 
This formulation could also be used when /’ #0 to 
ke care of the spin Another subdivision would 
then be made among the /-values, using only that 
part of (S31) that depends on the orbital momentum 

\s an example, the second approximation (1S) wall 
be evaluated for a @CE Kk term Step (a) 1 
ready carnmed out im 7 o that according Ito 
tep by 


The procedure is repeated for the terms on the right 


and the results obtained are substituted back in (32 
to vet the average energ, for the selection 

W (PCR )s °F sds Wids ?D 35 
This result could also be written by inspection 
The average energies (32 33), and (25) are inserted 
into (18 together with the trivial values for 
Wiysd®), Weysd), and Wy the sums over 7’ intro- 


rt ORNI Ous 





] 


) when 7? Is 


the J contivuration 


duce factors ( 


conhwuration 


iN i 7) 


Formulas 1) and 5) are valid even when con 


hivruration interaction 1 present 


a. Tae "6 


Configuration 


The final configuration considered 1 


the / and the 1? configuration are assumed to have 
SZ-values that are exactly pecified kor the selec 
tion 7’ that contains only one of the /’-electrons, 
the coetlicient ol fractional parentage l i ply 1 


transformation coefficient given by (6 


P-2(4SL)US8S,L,US8'' 1 SLAIN LOS’ 1 
SL.I72(S’ L’)S8"' I SL’(S. L038" i 
Ss Ss 1)(2/ / 
W (S3S’3:SS)OWUL'UV LL’ SS) 


This leads to a simple formulation for the /"~?s* con 
figuration d) The parentage of an electron in a 
term 1" @s’¢ySZ is divided in the ratio of the multi- 
plicities between the possible parent terms generally 


two 





For the selection 7’, that contains both /’ elec 


trons, the coefficient of fractional parentage becomes, 
S’L)S”’L’") 


(I? “aS, LVAS’L)SL ISUL’W YSL)I’ 


(/" (aS, L)ISLW YSL) 


SI LUSDS’L(S’L’) |S LS’L(S.LIS’L"). (39) 


The transformation coefficient is given by (5) of III 
and has a form similar to (31 In a term /"-? s*ySL, 
however, the closed shell of two s-electrons can be 
disregarded for this selection, so that (39) reduces to 
the procedure already given for a configuration of 
equivalent electrons 


To evaluate the average energy (6) for a d@s? °F 
term the formulation (d) is used 
Wid? Fd? 8 2/sWid? « *F P/3 Wd? s I 10 
iW d 8 I eds W /s ) (41 
Relation (41) is obtained by inspection; using (d) 
again with (41) leads to the average energy for a 
‘ hoice ds 
Wid? s? °F sds) —3/4W(ds*D) 4 1/4Wids '"D). (42 
The result (42) could be obtained also from (40) by 
use of formulations (b) and (e 
Wid? s *F-ds) — Wids*D 
Wid? s?F sds) = 1/4W (ds *D) 4+.3/4Wids 'D). (43 


When these values are substituted for the observed 
energies W on the right of (40) the result (42) is 
obtained: this checks the CONSISLERCY of the proce 
dures and illustrates the use of (c¢ The values (40 
$1), and (42), together with trivial values for W(y:d 


Wiy:), Wiysd 


mn ; : n—2 
(18). The number of selections 7’, of type / is ( ); 
/ 


and Wy:s) are then inserted into 


n—2 
I~ s ean be seleeted in 2 ( ; ) Wave and / ‘ 
in ("> ) ways 
Wid? 2 SF) *—4/3W (da? s *F 23 Wd? s ?*F)* 
2W (ds? ?D)*—3Wids*D 
Wids 'D)*— W(d? °F 
Ws? 'S)*+2W(d)+2W(s 14 


The following formula is obtained in the same wa\ 


W(d? 8? 8P)* = 4/3 W (de 8 *P) 4+2/3W (ae 8 *P)* 


+-2W ids? ?D)* 3Wids*D 
Wids 'D)*—W(2'G 
Ws IS) *-4 2UW id ou (8 5 
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The formulas for the d? s* configuration alway. op. 
tain terms that are likely to be perturbed by op. 


The s?'S term is usual! yp. 


figuration interaction 
of the formu ie 


known, and the applicability 
limited by this also 


4. Comparison With Experiment 


In all the formulas derived above, the observed 
energies, W, are referred to the core with all 
electrons removed. Because ionization pote 
are often not accurate, these absolute energies may 
not be well known. To check the theory of Bache; 
and Goudsmit only differences between terms hay, 
lonization potentials cancel in thy 
the 


low esi 


our 


als 


been calculated 
differences, and energies W can then be taken a 
relative-energy levels with respect to the 
term as zero; these are given in a Bureau Circular 
In particular, the following two differences would by 
unaffected by configuration interaction: 


D, = W,(d* 8 ®P)— W,(d8 8 °F) = Wd *P)— Wa 
W(d? s *P Wid? s *F Wd? ®P) + W(d? 5F 
16 
D, = W,(d* *D) — W,(d* 8 °F) =6/5W (da? *P) 


5Wid? s *P) 12/5Wie? s 4F 


-O/5W (ds *F 3 


6/5W(d? *P)—9/5W(? BB) +E 


W(d)— W(s) (47 
In table 1 the calculated and observed values of 
D, and DP, are compared In the spectra of titanium 
and zirconium sufficient data are available to com- 
pare the /,-values. The calculated and observed 
values are in close agreement, the calculated valu 
exceeding the observed value by 357 K in titanium 
and by 293 K in zirconium. Separations of terms in 
the same configuration are thus given closely by 
the theory. By assuming that /), will be given with 
the same error in more highly ionized spectra, it is 
possible to calculate positions for some terms that 
are not known experimentally, and use the values 
obtained to calculate DY), in these spectra.“ Accord- 
ing to table 1 the difference between the two terms in 
different configurations is predicted very badly by 
17) for titanium, the error in DP), being 6068 Kk 
However, the error is much less for the spectra of 
the higher ionization of vanadium (error 1841 K 
and it is still less in chromium (error 635 K). The 
error is also less in the homologous sequences of 
zirconium (error 2296 K), where the nuclear charge is 
much greater (Z=40) than in titanium (Z=22 
but the degree of ionization is the same. The 


C. E. Moore, 
(1982 

The @CP)s *P term is unknown in V mand Cr tv; the calculated positions 
these terms are given in table 1. If the errors in D, were to decrease rapidly * 
the errors in Dy do, then these calculated positions would be about 357 K too h 
and the calculated values of D, would be too small by about three-fifths of t 
error, or 220 K rhis is a small error compared to the total errors found in 
calculated values of D). However, it is noted later that the errors in D; proba 
will not decrease 
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improvement from titanium to vanadium to chro- 
mium is especially remarkable in view of the fact that 
the errors in 7), form a rapidly decreasing geometric 
sequence, Whereas /), itself starts at a large positive 
value and decreases to a twice as large negative value 

Racah “ has pointed out that the energies would 
be expected to vary more smoothly if the degree of 
ionization (rather than the atomic number) were the 
same in the different spectra. The results obtained 
by use of this empirical modification of the Bacher 
and Goudsmit theory are given in table 2. Except 
for the first spectra, the errors in calculated values 
of D, are on the order of 100 K; part of the good 
agreement must be fortuitous, however 

\ partial explanation of the poor agreement ob- 
tained in the first spectra may be that the d* *P (of 
Se rand Y 1) is perturbed by interaction with the 
term dsd’f*P (in the spectra with 4d- and 5s-electrons, 
this has the added effect of making the calculated 
position for the d? °F too high, as seen from the small 
separation of the d?°P and d@? °F in Sr 1 This effect 
can be estimated in the spectra containing 3d- and 
ts-electrons by calculating the position of d *P so 
that the error in 7), will be 170 K as it is in the second 
spectra. When this calculated position (40048 K) 
is used, the error in J), is still very large 5769 K 
Racah has noted that the observed position of d®? D 
in K tis about 5117 K lower than expected from the 
behavior of the difference between baricenters of the 
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d"~' s and d" configurations in Ca 1, Se 1, Ti 1, and 
V 1, and this may explain the poor agreement ob- 


tained in first spectra 


5. Relations Between Radial Parameters 


All formulas of the Bacher and Goudsmiut theory 
ean be used also to calculate the matrix elements of 
sealar operators ; this has been noted Ih speciite Cases 
in the discussions of formulas (21) and (22), and it 
was also pointed outin BG If the term values could 
be expressed in terms of such matrix elements, and 
Same for all 


if the zero-order functions were the 
contiigurations then the formulas derived from 16 
to (19) (such as (46) and (47)) would be identities if 


the H's were replaced by matrix elements of these 
operators 

The matrix elements of the operators are 
combinations of a limited number of radial integrals 
For the Slater theory, these combinations are well 
most of the formulas necessary for the pres 
The obser ved sepa 


linear 


known 
ent paper are contained in I] 
rations of terms in the same contiguration are pre 
dicted well by the Slater with an accu 
racy equal to that obtainable with the Bacher and 
an empiri al correction pro 
account for the 
Radial 
fune 


theory (i. e 


CGoudsmit theory) uf 
portional to LOL l 
polarization energy, 

delined by 


is added to 
see footnotes 8 and 9 
integrals zero-order one-electron 
tions having the same n- and /-values must 
sidered independent of each other when in different 
however: 
directly 


be con 


configurations to get this good agreement 

if the Bacher and Goudsmit theory were 
applicable all such integrals would have to have the 
same value. The derivation of the theory, given in 
shows how much of the observed variation 
integrals will be accounted for by 
CGoudsmit theory in riven 


section 2 
of these radial 
the Bacher and 
approximation 
To illustrate, formulas for the W's that enter into 
(46) are given here according to the notation used 


in section 2 


any 


’ 
lbtald s 


bhatd’ s 


W [ ‘i! | / lia d 

Wid J id lid L5Bid bald 

Wie stP 1(d? s TBid: s (dd? s lOal(d s 

Wid? s ‘I Ald? s SBid? s »t),(d? s 

Wid? *P lid TBid l0a(d 

Wid? * lid SBid {8 
i LiL+ | ‘ la " | t 1 

aprerat ' > The value S for a t tier 

! I 

isall(L +1 alii+t In the present work, this, rather, thar ply ali L+1 

is acided to each term; this change in the form of the correctk compensated by 

a change in the value of the constant parameter A rhe mod ‘ has been 

made so that the relationships obtained later between the 4-walue n different 


configurations will have a simpler fort 





, then the 
is an identity if the /’s and a’s in different confi . 


If these values are substituted in (46 


tions are assumed equal. Relation (46) will sti be 
satisfied if the following two less restrictive cond nea 
hold 

Bid’ s Bid Bid? s Bid 

a ds ald a d a d ) 
It is sufficient for /—2/3a@ to satisfy a relation of 


similar form for (46) to be satisfied, but if separa is 
of all pairs of terms in the d@* s configuration are to be 
correct, then both relations must hold and in addition 


the following two 


Formulas of this tvpe can be derived easily 

Inspection, and there ts no need to considet fractional 
To obtain the second equation In (50 
note that in (1S) the combination of a 
/- and an s-electron (needed to define G 


- it occurs two wavs in each of three 


parentage 

for instance 
occurs three 
wavs in Wild's 
possible selections / that contain the s-eleetron (in 
the d? s configuration 
the three possible selections 7’ 

s-clectron (in the ds configuration 
relation for the parameter @, is concerned, this dire 
interpretation of (18) leads immediately to (50 
multiplied by a factor of three As B,C, 
ure all defined by the combination of Iwo d-electrons 
a similar derivation leads to the same form of equa- 
tion for all three of parameters, and = this 
equation is the one represented by (49) and the first 
In this wav it can be shown that 
CGoudsmit approximation 
will predict correctly the 


- it occurs only once in each of 
that contain the 
As far as the 


and a 


these 


equation of (50 
the Bacher and 
relation (18) with n=4) 
separations of terms in the same configuration if the 
following relations hold between the parameters 


second 


Bd‘) =2Bi(d Bid 5} 
Bid x Bid ? Bid? s Bid 52 
Bid s 2B(d* s) Bid 53 


Three similar relations must be satisfied by the C's 


apd three more by the a’s. Also, it is necessary 


G.(d N 2G,(d* N (,(ds) o4 

The derivation from (18) fails to bring out the fact 
that the two-electron parameters must have a linear 
variation if the formulas are to be valid. In general, a 
two-electron parameter defined between electrons 
with orbital angular momentum / and /’ in the 
configuration y will be specified by (y)G@(W’). Ae 
cording to the derivation with the second series of (9 
only a first-power variation with the number ol 





























the Bacher 
In particular 
parameters #, Cor @ must 


electrons thoy can be accounted for by 
d Goudsmit second approximation 
en. anv one of the ; 


itisfv a relation of the form 


s eusv to verifv that (51) to (53) are satisfied 
hen a formula of this kind is applicable 


ar spectra rather 


possible to creneralize > to include all 


than just those of the 1OnS 


the same atom, mn two thre and fow electron 
ectra If a7 is the degree of ionization this veneral- 
on will have the form Gidd i+ rl 
( iD The difference 1), defined by 1G ana 
int to 15/ lOa@ according to the linen theory 


houutd niso be ( Xpre ssible by i formula ot this U\ pe 
i eal has fitted 4 experimentally observed [?.] 


lerences mm spectra with Te nnd tx-clectrons with 
mean error of 340 Kk by use of the formula 

D, = 2679 +853 13j4h sLSSq—11lg nh 

rl formula is adjusted to spectra with eight 

different tonizations, and without a quadratt term 

the mean erro would bye Inereas a lo about 5?0 kK 

According to t| Is formula thre differences iD eal 


culated trom ti will bye too great on the 
ly about 220 Kk \Modifving thre 


spectra with the same lonization are 


nverave 

theory so that 
utilized tends 

This could also 

approximation 

lhe original series (9) will often have to be satisfied 


to eliminate this svstematic error 


be accomplished by using the i 


for a given y, by more than one value of a, and this 
leads to the conclusion that the ratios of similarl 
defined radial parameters in different spectra of a 
given element should bi In particular, the 
constant in configurations 


eq al 
following ratio must be 


with d- and s-electrons 


constant (57) 


This is not really a necessary consequence of the 
theory deve loped in section 2, however, as the series 
9) can be reinterpreted in terms of radial parameters 
that independently ) 

When irrational roots are considered, it is shown in 


satisfy relations such as (55 
section 7 that (57) is a necessary condition for con- 
sistency of this interpretation, in terms of a linear 
theory, with a direct application of the method of BG 

In the same way that the differences such as (46 
were replaced by the formulas (51 to (54 the 
differences similar to 17) can be replaced by rela- 
tionships between the A values. The Bacher and 
Goudsmit second approximation will give the sepa- 
rations of terms in different configurations correctly 
if the following two relationships between the A’s 
are satisfied 


A(d*)— A(@ 8s) =3A(d 3A(d s)—3A(d*) 





A(s) +-2G)(d s 


By redefining the sl-values suitably?! the dependence 


on G, may be omitted from (59); equations (5S) and 


59) mav therefore be regarded as two relations 
between the differences A(d"* s Aig?" In 
the second approximation, the parameter .1 must 


relation, and these differences then 
Ssatisly (yl adrath equations having the general form 


satisfy a cubs 


A(a’* # A(d s' D+ kk Fr + Gkr 3} Th 


0) 


The i dependenes s omitted beentuse il ndds nothing 
extra when & has only the two values 0 and 1 it is 
easy to vernfy that (58) and (59) are satisfied when 


a formula like this applicable This formula can 
be generalized to apply to all similar spectra ina 
sl for 55) af the “uA . 
suitably defined im the spectra with more than four 
electrons If fourth 
important, a constant (or linear 
The errors of 7), in table 
explained in this way 

In the one-electron 
electron spectra considered in section 4 


manner similar to that 


third (on power terms were 
errotl would result 
| change too rapidly to be 
however 

three 
there are 
five experime ntally determinable differences between 
AA-values, and the five unknown constants in (60) 
sim 


two-electron and 


can be evaluated from these observed values 
larlv, the three constants in (55) can be evaluated 


from observed separations of terms in the two 
three-electron spectra The Bache 


and Goudsmit second approximation ts then equiva 


elec tron and 


lent to using the formulas obtained in this way to 
extrapolate values in the four-electron system 

This equivalence has been verified by using all 
available experimental data in titanium and = zit 
conitum to evaluate the radial parameters for all 


spectra with least squares The results confirmed 
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the conelusion already obtained from the limited 
amount of data utilized in table 1. In all spectra 
of both elements, the terms could be fitted with 
absolute errors less than 400 K. The parameters 
obtained could be adjusted to fit the relations con 
tained in (51 without significantly changing 
the original agreement. This verifies the conclusion 
that the separations of terms in the same configura 
tion are given by the Bacher and Goudsmit method 
with an absolute error less than 400 K in titanium and 
zirconium. The discrepancy in calculating A(¢/) 
A(d@ s) from relation (58) was found be 6560 Kk 
in titanium and 2086 K in zirconium. These values 
are approximately the same as the errors found for 
D, in table 1, which verifies the conclusion that 
calculated separations of terms in different configura 
tions may be considerably in error 

If the ionization potentials are accuratley known 


to (57) 


lo 


the aboslute values of the A’s can be calculated 
(rather than just the differences (58) and (59 

L(d/' bAtd HAC LA(d) +> 6 (61 
Lid’ s Lid Ad BAC 

,A(ds | Ad 6 62 

Lid? s 2A 2 A(ds 1(d?) —4A(ds A(s*) 

2A(d) } 2A(s) — 26,(ds 6), (ds) + 6 63) 
where 

6= /,—31,7 31,—1, 


The ionization potentials appear only in the combi- 
nation 6; /, is the principal ionization potention of 
the four-electron spectrum (71 1, for instance); J, 
is that for the three-clectron spectrum (as 7'i 1), ete 
If (61) to (63) are utilized, errors in the theory will 
not cancel as they may in the differences (58) and 
(59). In titantum the error of (61) in calculating 
A(d*) is found to be 24130 K, and the error of (62) 
in A(d@* 8) is 1757 K; on the other hand, the corre- 
sponding errors in zirconium are — 2100 K and — 4186 
K. The inconsistent behavior of the errors in these 
two calculations probably arises from the well-known 
need for better ionization potentials, particularly 
in the one- and two-electron spectra 


6. Weak Configuration Interaction 


The proof of the method of Bacher and Goudsmit 
is established with second-order perturbation theory, 
so that the theory is not directly applicable if con- 
figuration interaction is too strong to be evaluated 
with perturbation theory. However, it may not be 
directly applicable even if the spectrum has weak 
configuration interaction that can be evaluated with 
perturbation theory 

For example, consider interaction between terms 
of d* and those of d"~' s and assume that second-order 
perturbation theory can be used. The perturbation 
of, say, a d° *G term by a d? s *G term is evaluated 
in terms of a characteristic interaction parameter 
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s configurations: t| r 


and d 
s) in agreement with the 1 
For arbitrary n,. 
parameter defined between d" and d"™! s. Pay fy 
with the help of table XX of IIL the following ) 


perturbations are evaluated 


associated with the d 
designated Id d 


tion in (SO) of III there is a sin 


— 5OULAd d? s)| 
A(0"G) = W(as2G)— Wd? «*G) 
140/71 (d?- ds)|? 
(4 ) ; 
ae Wid?! D)— Wids 'D) 


By considering terms marked with asterisks in (2s 
it follows that configuration interaction between 


and d"-' s will be correctly accounted for by thy 

Bacher and Goudsmit theory if the following relation 

exists betyv een the two perturbations (64) and (65 
(> 5A d (j 3 7A -P 'T)) 8] Hh 


According to the perturbation calculation used by 
Bacher and Goudsmit, the energy difference W 
Gi Wie? s *G) would equal the difference WU 
'D)— Wids 'D) and the parameter //,(d°—d? s) would 
equal /7,(d?—ds If these two conditions held, then 
(66) would clearly be satisfied. Relations similar to 
(66) can be found from other formulas (such as (29 
(36), and (37 By assuming that all energy differ- 
ences Wid" aSl)— Wd"! se’SL) are equal, and all 
parameters //,(d"—d""' s) are identical, matrix ele- 
ments of configuration interaction can be calculated 
with the Bacher and Goudsmit theory by solving a 
set of simultaneous equations, although such a pro- 
cedure would not be the most practical. Similar 
considerations would apply to any pair of configura- 
tions; in particular, the interaction between d" s and 
d"~* s* is interesting because the parameter defined 
between these two configurations would be equal to 
the parameter defined between d" and d""! 
ing to this direct application of the Bacher and Goud- 
smit theory 

Alternatively, (66) could be satisfied if the para 
meters 11d ?? s) and H1,(d? ds) accidentally had 
the same ratio as the square roots of the energy 
denominators in (64) and (65). To get consistency 
for all formulas similar to (66) that could be derived 
it is necessary that a representative average energy 
difference exist that approximates separations of 
terms with the same S- and /-values in the two 
configurations (in many cases no average value will 
adequately represent the data). An easily defined 
value for this average is the difference in the A-values 
of the two configurations. In table 3 calculated 
H,-values and differences between A-values have 
been collected. In a given atom, the differences 
between A-values in this table are greater the 
greater the ionization; /7/, is also larger in the spectra 


SN accord- 


_ of higher ionization, according to the results obtained 


for vanadium. 

An extension of this argument indicates that the 
polarization energy arising from interactions with 
configurations in which only one electron is excited 
may be partially included in the second approxima 
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on of BG. However, it is unlikely that the inte: 
ction parameters would behave exactly as required 
«0 that it might not be well accounted for 
It is clear from table 3 that in the titanium spectra 
of moderate ionization, the Bacher and Goudsmit 
method would not aecount for the configuration 
interaction Even though the interaction is weak 
enough to be evaluated with second-order perturba 
tion theory the calculations discussed in section 5 
show that errors of 300 to 700 K would be introduced 
into formulas (28), (29), (36), and (37 Whenever 
the average term separation is small or negative 
which is the case in Titand Ti mu, no set of values for 
the interaction parameters can produce agreement 
When the configuration interaction is so strong that 
perturbation theory cannot be used, the errors are 
even larger. For zirconium, errors in excess of 1000 
K result when the formulas are applied directly 
The differences in A-values are all large and positive 
in chromium, so that the configuration interaction 
would probably be well accounted for in the second 
approximation; in this case, however, the interaction 
is weak and could be neglected in the linear theory 


7. Irrational Roots 


This section deals with one of the simplest exam 
ples in which irrationality of the energy levels plays a 
part As a preliminary, the two *D terms of the d 
configuration are discussed 

The determination of the energy levels of the two 
D) terms of d* was first carried out by Ufford and 
Shortlev.“ They chose two mutually orthogonal *?D 
eigenfunctions, which they characterized as a *D and 
6 *D; in the more recent notation of IT] these same 
two eigenfunctions are designated as *;D and *,D, 
respectively In terms of this choice of eigenfunc- 
tions, the electrostatic energy matrix was evaluated 


C. W. Ufford and G,. H. Shortley, Phys. Rev. 42, 167 (1982 See alsc, EF. t 


Condon and G. H, Shortley Theory of atomic spectra, p. 233 (Cambridge, 1951 
In ILI the classification of terms in d* is carried out by introducing the senior 
imber; this is given as the prefixed subscript in the notation above 


/ lie iif lia 18 / 
4 
oun 
and found to have the form 
v ae oe. ae 7 
D\3,21B A+3B+30C ° 


The energy eigenvalues of the two *D states of d° are 
the characteristic values of this matrix: these values 
as well as the observed energies for the 7D terms 
are specified by WeD*) and WeD 

The matrix S that diagonalizes (67 
form 


will have the 


S ( i yo ( | (GS) 


SMS ( 


) 69) 


The elements of S can be specified by the condition 
for the vanishing of the nondiagonal elements of 
SMS 


l4B(d (dP) JC\C ( C')38y 21 B(d 0. (70) 
Although WGD) and WD) have no physical reality, 
they will be subsequently to indicate the 
energies that would be associated with the diagonal 


elements of the matrix (67 The following formulas 
are obtamed from (69 


used 


WGD)=—Cc? WeD cou 


If the polarization energy can be accounted for with 
the L(L+-1) correction, then all the preceding dis 
cussion will still apply even when polarization energy 
is considered if the constant A in (67) is suitably 
redefined. 

By direct application of the methods in section 3, 
the energy for, say, a d‘°G term, can be expressed in 








the form 


Wd? °G)* = 10/21 W(d* {D)* + 3/4 Wd? PRY 

1/3 W(d* *F)4+- 99/140 Wd? ?G)" 

11/15 W(d* 7H 11/10 W(d? *P 

0/14 W(d? 'D)*— 29/10 W(d? °F 

19/14 Wd? 'G tWid 73) 
If the linear theory is valid, WG@D) in (73) can be 
replaced by the value (71 the approximation 
obtained in this way will be specified by a prime. In 


the exact application of the Bacher and Cioudsmit 


theory, (71) can still be used to replace Wid 2D 
but C) and Cy are not necessarily specified by (70 
Also, instead of evaluating Wid 2D?) which was 
done in deriving (73), it is necessary to evaluate 
Wihic I) ( ) | d The following is the 
correct formula whi hous obtained in this way 

Wald *G)=— Wad *G e(d'* °G 74) 
e(d'°G C? C38/21 Wid? 'S 2/7 W(d? §P 

50/147W(@2'D 2/7 Wd? °F) 

('¢.¢ (' 
2/49 Wd? 'G) 21W(d? *P) 
l4y21 

15W(d? 'D 21W (d? °F 15 W(d? 'G) 
Replacing the @ energies by their Slater formulas 
(including the L(/ | correction leads to the 
following approxunate form for (75 
e(d* °G)= 10/210, C,{ [4 Bid tC) CC (C3 

C3)3y 21 Bid’) (76) 


It follows that when ©, and ©, are determined from 
(70) then (74) will be equivalent to (73) if the ratio 
of B to Cis the same in @* as it is in d To establish 
the general correctness of the method of applying 
the Bacher and Goudsmit theory which was devel- 
oped in section 5, it is necessary that in general the 
ratios of similarly defined radial parameters be the 
same in all spectra of a given atom 

lt is well known that the ratio of 7? to ¢ 
nearly constant Values of this ratio, 
mainly from the calculations discussed in section 5, 
are given in table 4 Because of the large errors that 
may be present im the theory (+400 K), the values 
of r,, obtained by fitting the observed data exactly, 
are not well defined; to a lesser extent the d° and d x 
configurations also have too few terms to define ry 
precisely The data as a whole (giving strong weight 
to the two well determined r;-values) support the 
conclusion that r, is a slowly decreasing function of 
n. The exception apparent in titanium 
noticeable if the irrational roots are omitted from 
the least-squares calculations for 3d° and 3d° 4s 
parameters (the error in fitting the irrational terms 
is then excessive, however 

To estimate roughly the effect 
from inequality of this ratio, 


is very 
obtained 


is less 


that could result 
can be put in the 


48 


following form with the help of (70), and by use n 
approximate value of r equal to 0.29 


e(d* °G)= 40/21 (" : 2) 02 C3 Cd?) 


1.1607 r)C(d*), 7 
In Ti wm, €(d 2505 so that if r;—r,—0.05 r 
€ 145 hk This is likelv to overestimate the 1 


imum error that can result, vet within the accuracy 
of the theory it probably could be neglected HH 

the same treatment for the d' 'F term s vs 
that « will have exactly 3 times the value (75 
for the d' *D term the error will be 2.4 times is 
value, so it is not certain that this effect will always 
be negligible The ratio GC 
(it varies considerably, from 0.4 to 0.8 usually 


ever, 


need not be const nt 
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Though the method of BG introduces new non- 
linear effects if the parameter ratios are not equal, 
these effects would probably represent errors in the 
method. Because of this, and the ratios 
are so nearly equal, it is not expected that a corree- 
tion with this variation will produce improved 
agreement when added to the linear theory How- 
ever, irrational terms often cannot be fitted by the 
linear theory with the accuracy that is obtainable 
for rational terms * and some additional correction 
of this form may be needed. This correction would 
not change the agreement obtainable for the rational 
terms of d', but it would introduce corrections into 
the irrational terms. These considerations may be 
less academic for spectra with p- and d-electrons 
where there are more parameter ratios that may 
differ, and where most of the terms are irrational.” 


because 
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* Meshkov has introduced a method which combine ‘ later theory with 
+ Bacher and Goudsmit first approximation, and has obtained improved agree 
nent for the 3d? 4p conficuration of Ti u (footnote 1s A simple extension of 
the linear theory to inel correlations between p- and d-electrons may explain 
this improvement, but the possibility that new types of crorections are needed 
to explain it, cannot vet be excluded. His theory also leads to better agreeme 
than is obtainable with the Slater theory alon 1 the @ configuration of V 
und the d* configuration of V u (Phys. Rev. 88, 270 (1954)) In the d* configura 
tion his final mean deviation of +663 K is still rather large; this may be becat 
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Self-Ignition Temperatures of Combustible Liquids 
Nicholas P. Setchkin* 





The methods pre vuusly used for determination of the ignition temperature f liquid 
ire briefly reviewed, and the large discrepancies in reported values noted Various faetor 
‘ ( affect the letermined value of twnitior temperature are disecu i A practical 

tion apparat ‘ ting part of ar sulated spherical flask, designed to provide 

‘ ditions favoral low twhition temperature values, w= described, and data obtained 

ith this apparat on a considerable number of combustible liquids are presented The 

isefulme tt perature curves of a thermocouple junetion placed in mixture of 

‘ nt tibsle ia lair, aS Aa mea f studving the preignitio reactions in such mixture 
istrat ! 


1. Introduction 


Che fire hazards associated with the use, storage 


nd shipment of combustible liquids are so familiar 


to need no discussion, and the tremendous and 
varied consumption of such liquids carries the 
roblem of protection against these hazards into 


For the develop 
vent of effective protective measures, as well as for 


everv element of modern economy, 


eflicient use of these liquids, extensive and accurate 
nowledge of their ignition characteristics is highly 
important 

Although standard methods of 
determination have for many of 
the physical and chemical properties of materials, 
such as density, viscosity, solubility, there has 
the definition or 
determination — of temperature C'on- 
sequently, determinations of this quantity by 
different investigators have varied widely both in 
numerical value and significance, with much result 
ant confusion and difficulty in the interpretation 
and practical application of the information. 

Preparatory to the experimental work reported 
in this paper a study made of the 
methods previo sly employed for the determination 
of ignition temperature lt appeared that, although 
the ignition characteristics of a material were 
doubtless a function of some actual property of the 
material, the measurement of that property by any 
means readily available would be creatly affected 
by the ambient conditions Hence, any practical 
definition of ignition temperature which might be 
adopted must continue to be on a careful 
definition and standardization of the test apparatus 
and procedure It was further concluded that a 
criterion that indicated the lowest temperature 
producing self-ignition of the material under condi- 
tions as favorable as might normally occur, would 
be the useful in the practical analysis and 
elimination of hazards, and would also provide a 
suitable basis for data of research interest. 

Accordingly, the apparatus and test procedure 
deseribed in this paper were designed to determine 
values in the lower temperature range in’ which 
ignition may be possible. The method developed 
appears well suited to the purpose, and in tests of a 
large number and variety of materials it has given 


and 
established 


definitions 
been 


etc., 


been no general agreement on 


iwhnition 


was Various 


based 


most 
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results of satisfactory reproducibility The present 
paper covers the work carried out over a period of 
about 15 vears a considerable number of 


on conm- 


monly used combustible liquids 


2. General Considerations 


2.1. 


Previous Determinations of Ignition 
Temperature 


have been used in the 
may be 


Most of the methods that 


determination of ignition temperature 


grouped into two general types, the “‘rising-tem- 
perature” method, and the ‘“constant-temperature”’ 
method. In the former, a combustible mixture is 
introduced into the ignition apparatus, and its 
temperature is then continuously increased until 
positive evidence of ignition is observed The 
self-ignition temperature is usually taken as the 


point at which the time-temperature curve for the 
material definite 
the point of intersection between the time-lempera- 
ture curve for the material and that for the heating 
medium {1].' In the constant-temperature method 
of determination, the ignition apparatus is preheated 
to a selected initial constant temperature before the 
specimen is introduced, and no adjustment of heat 
input to the apparatus is made after introduction of 
the specimen The repeated with fresh 
specimens at different initial temperatures until the 
temperature at which ignition is produced 
has been established It has been determined as 
the temperature of the heating bath surrounding the 
reaction chamber, as the temperature of the com- 
bustion chamber wall, or as the temperature of the 
air within the chamber 

Within two the 
problem, numerous and diverse methods of determi 


shows a iInerease nm slope or as 


test is 


lowest 


these veneral approaches to 


nation have been developed (C. R. Brown [2] gives 
references to 16 different methods of the rising 
temperature type alone Following are brief 


outlines of a number of procedures that have been 
used, and the data in tables 1 and 2 illustrate the 
wide variation im values obtained by different 
investigators, 

Oil-drop method. Early determinations of ignition 
temperature with the so-called ‘‘oil-drop method” 
were made by Holm [3] in 1913 with a heated 


is paper. 








Tarte 1. Jgnition temperature (in degrees centigrade) of 
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Investigator A) Callendar; (B) Lewis; (C) Estrader D) Prettre; (E 
porcelain surface Moore [4] used a_ platinum 


crucible submerged in molten metal or embedded in a 
steel block. Sortman and Beatty [5] made a small 
cylindrical cavity of 43-ml capacity in a stainless- 
steel block and provided for the circulation of 
preheated air. Jentzsch and Von Zerbe [6] also 
used a steel block with four cylinders (7.1 ml) and 
circulation of oxygen or air. A modified type of the 
Moore apparatus, heated electrically, has been 
widely used by many investigators [7]. The initial 
temperature of the bath, or wall of the ignition vessel, 
leading to ignition is usually considered as the self- 
ignition temperature. The values given in table 1 
are reproduced from reference [7] 

Dynamic tube method. In this method, a Pyrex 
tube is heated in a conventional electric furnace to 
the desired temperature indicated by a thermocouple 
in contact with the outside surface of the tube. A 
combustible substance, being preheated separately 
or together with air or oxygen, is allowed to pass 
through the tube. The tests are repeated at differ- 
ent temperatures of the tube until evidence of igni- 
tion is observed. Apparatus of this type was used 
by Dixon and Coward in 1909 and 1934 [8, 9], by 
Mason and Wheeler [10], and by many others. 
Dixon and Coward limited their test to 15 seconds. 
Callendar [11] and Moureu [12] made a study of slow 
combustion and determined the ‘““Temperature of 
Initial Combustion” (T. I. C.) in an improved ap- 
paratus. This temperature was defined as the tube 
temperature necessary to provide the first appear- 
ance of products of combustion (steam, carbon 
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sids determined with different methods [15] 
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dioxide, aldehyde, ete and determined by 
chemical analysis of samples taken periodically 
during various phases of the reaction. 

This temperature has some relation to the “Tem- 
perature of Thermal Reaction” (T. T. R.) intro- 
duced by Estradére [13], and the ‘“‘Critieal Inflection 
Temperature” (C. I. T.) introduced by Lewis [14 
All these temperatures are related to the ““Tempera- 
ture of Initial Oxidation” [15] or, in other words, to 
the temperature of incipient self-heating. Prettre 
defined the ignition temperature by means of two 
criteria, the appearance of “cold flame” (T. C. F 
and the appearance of normal flame at “Minimum 
Ignition Temperature” (M. I. T.). A tabulation of 
ignition temperatures measured by this method is 
included in table 2, which is reproduced from 
reference [15] 

Evacuated-vessel method. In this method a com- 
bustible mixture is suddenly admitted into an 
evacuated vessel preheated to the desired tempera- 
ture. This method was originated by Mallard and 
Le Chatelier in 1880 and was used subsequently 
by Mason and Wheeler [10] and Townend [16] 
Ignition was indicated by rapid rise in pressure 
Lewis and also Pidgeon and Egerton [17] modified 
this method by admission of the cold combustible 
mixture into an evacuated vessel heated at a stead) 
rate (rising-temperature method). 

Bomb method. The bomb method was used by 
Neumann [18] and other investigators for determina- 
tion of the ignition temperatures of combustible 
mixtures at high pressures. 
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Adiabatic method. Heating of a 
ymbustible mixture to the ignition point by adia- 
atic compression was used by Tizard and Pye [19] 
nd Pignot [20]. It is hardly likely that ignition 
mperatures measured by this method, which 
volves the use of rising temperature and pressure, 
‘ directly comparable to measured at 
onstant pressure 
ASTM _ standard method. Committee D-2 on 
troleum Products and Lubricants of ASTM in 
30 adopted a standard method and apparatus 
D286—30) [21] for the determination of the ‘“Auto- 
neous Ignition Temperature.’ The standard ap- 
iratus consists of a conical Pyrex flask of 160-ml 
pacity, submerged in molten solder heated by a 
s burner. A few drops of the liquid to be tested 
allowed to fall into the flask, previously heated 
» different initial temperatures. The ignition 
mperature is defined as the lowest temperature of 
solder bath, indicated by a thermocouple placed 
a silica protection tube submerged in the bath 
ar the flask, necessary to cause ignition of the 
pecimen 


com pre SSi0n 


those 


2.2. Selection of Method and Definition of Self- 
Ignition Temperature 


With the rising-temperature procedure of deter- 

nation, the ignition temperatures vary with the 

te of temperature rise and with the criterion of 
gnition Elaborate control equipment Is necessary 
to obtain reproducible results either on the basis of 

fixed rate of heating or on that of maintaining a 
fixed temperature difference between any part of the 
specimen and its surroundings. Furthermore, the 
omparatively slow heating of the specimen in this 
method is not generally well suited to the study of 
liquids, many of which apparently produce vapors 
that form relatively stable compounds with oxygen 
at low temperatures. If the system has time in a 
low-temperature range to permit the formation of 
such compounds in appreciable amount, the observed 
ignition temperature may be well above that which 
might have been obtained with the constant-temper- 
ature method. In this latter method, on the other 
hand, precise control of the initial temperature 
may be accomplished with relatively simple equip- 
ment, and the rapid heating of the specimen limits 
the time in which low-temperature reactions may 
take place. The constant-temperature method of 
determination was selected for the present work, 
and the liquid charge was sprayed into the heated 
chamber as a fine dispersion, so that it rapidly 
became vaporized and mixed throughout the air in 
the chamber. 

The characteristic thermal effects when a combus- 
tible liquid is thus introduced into the chamber at 
different initial temperatures are illustrated by the 
time-temperature curves of figure 1. After injection 
of the dispersed charge, the temperature within the 
chamber first falls below the initial temperature be- 
cause heat is absorbed in vaporizing and heating the 
ample. In a short time, heat supplied mainly from 
he walls of the chamber reestablishes the initial 
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and nonignition reactions of comb istible liq 41ds 
temperature. At relatively low initial temperatures, 
there may be no further thermal effects, and the 
temperature of the sample-air mixture may remain 
constant at the initial value, as at to. At higher 
initial temperatures, exothermic reactions within the 
sample-air mixture may be induced and the tempera- 
ture of the reaction mixture raised above the initial 
value. The extent of this temperature rise depends 
upon the difference between the rate of heat genera- 
tion within the mixture and the rate of heat loss to 
the chamber walls. At initial temperatures that 
produce relatively mild reaction within the mixture, 
the heat generated is eventually lost through the 
chamber walls without visible effect on the mixture, 
and the temperature returns to the initial value after 
a period of time dependent on the intensity of the 
exothermic reaction. This situation is illustrated by 
initial temperatures ¢;, fe, ts, and ¢, in figure | How- 
ever, at temperature, characteristic of the 
material under study and of the particular test cham- 
ber emploved, the reaction within the mixture will 
become sufficiently intense to raise the temperature 
to the ignition point, as evidenced by flame or ex- 
plesion. The initial temperatures ¢; and f, represent 
such ignition reactions and further illustrate the 
reduction in ignition lag, which usually results as the 
initial temperature is increased ee ail 

For the present study, the following definition of 
self-ignition temperature was adopted: 


some 


The self-ignition temperature (SIT) of a liquid is the 
minimum temperature from which a combustible mixture 
with air will pass into a condition of observable combus 
tion in the experimental apparatus emploved 











2.3. Factors Affecting the Observed Values of SIT 
a. Thermal Characteristics of the Apparatus 


Inasmuch as ignition depends upon the tempera- 
ture reached by the reacting mixture, the governing 
factor in its development is the amount of heat re- 
tained within the mixture and available for increas- 
ing the temperature. Hence, ignition will result only 
when the reaction is intense enough to generate heat 
more rapidly than it is lost from the reaction mixture 
to its surroundings, and the initial temperature of the 
chamber must be high enough to induce such a 
reaction Therefore, the higher the rate of heat loss, 
the higher will be the initial temperature required 
to accelerate the reaction to ignition, that is, the 
higher will be the determined values of the SIT 

Thus, the design and construction of the reaction 
chamber affect the SIT values that will be obtained 
In order for the determinations to approach the low- 
est temperatures from which there may be danger of 
ignition, the chamber walls should be constructed 
with as low heat capacity and heat conductivity as 
may be practicable. The size and shape of the 
chamber, as well as the structural materials, must 


also be considered. 
b. Pressure and Uniformity of Temperature 


The relation between the pressure within the 
ignition chamber and the observed value of the SIT 
has not been determined However, for the develop- 
ment of a practical test method, it is reasonable to 
assume that atmospheric-pressure variations will 
have little effect on the determined values of the 
SIT. 

It is probable that the ignition reaction between 
the combustible substance and the oxidizing agent 
initiates with the most highly energized molecules of 
the mixture, and supplies the heat necessary to raise 
the energy level of adjacent molecules to the ignition 
point, thus propagating itself throughout the mixture 
(22). If all of the molecules present are very close 
to the required energy level at the same instant, the 
spread of ignition through the mass will be extremely 
rapid. Hence, the maintenance of uniform tempera- 
ture throughout the reaction chamber will favor rapid 
propagation of the ignition reaction 


c. Reaction Mixture Ratio 


The ratio of the combustible substance to air in 
the reaction mixture has an appreciable effect upon 
the SIT. Mixtures that are excessively rich or 
excessively lean usually ignite at higher temperatures 
than those close to some intermediate composition 


characteristic of the material. Carbon and soot, 
which may have been deposited from previous 
determinations with over-rich mixtures, may also 


influence the SIT value obtained. 


d. Method of Temperature Measurement 


Temperature measurements are usually made by 
means of one or more thermocouples either embedded 











in the heating bath or medium outside the cha 
walls or installed directly inside the reaction 
ber. Reproducible results may be obtained 
any of these methods, but because heat loss fro: 
chamber cannot be conveniently eliminated, 
values of SIT are usually observed with th 
arrangement Furthermore, a thermocouple p| 
inside the chamber may also be used to ind 
temperature changes during the reaction. 


e. Criteria of Ignition 


The intensity of the visual evidence of ign 
varies decidedly with the initial temperature { 
which the reaction proceeds. At high initial 
peratures, self-ignition usually occurs in more or less 
violent form, such as a bright flash of flame and 
audible expulsion of gases (explosion At somewhat 
lower initial temperatures the flaming becomes r 
tively mild, whereas at the lowest temperature from 
which ignition is detected, the flame may be very 
dull violet or blue, often visible only under dark- 
room conditions. The latter may be confused With 
cool flame reaction, which is not generally recognized 
as evidence of ignition even when, occasionally, it is 
accompanied by a rapid rise in the temperature of 
the reaction mixture. 

The time lag before the development of ignition is 
also affected by the initial temperature of the reac- 
tion. It is an inverse function of the initial temper- 
ature and may be only a fraction of a second at high 
initial temperatures and may amount to a number of 
minutes at the lowest temperature from which 
ignition occurs. For example, methyl cyclohexane 
showed an ignition lag of 30 seconds at 275° C, and 
using the same test equipment, the lag was 20 min 
10 sec at 248° C. By limiting the tests to 10 or 15 
seconds, as some investigators have done, the self- 
ignition temperature would have been determined as 
well above 275° C. 

In addition to the visual evidence of ignition, the 
time-temperature curves for many liquids provide 
supporting indication of its occurrence. The curves 
at ¢; and ¢, in figure 1 illustrate typical differences in 
the curves for ignition and nonignition reactions 
proceeding from only slightly differing initial tem- 
peratures. Although the rate of temperature rise 
was at first about the same in the two cases, it soon 
became much more rapid and pronounced in the 
ignition reaction and was followed by a sharper drop 
from the peak temperature: This decided difference 
in the time-temperature curves for ignition and non- 
ignition reactions is often helpful in the final deter- 
mination of the SIT, particularly when the detection 
of flame becomes difficult and somewhat uncertain 


t. Extraneous Stimuli of Ignition 


The ignition reaction may be induced by a number 
of means extraneous to the reaction mixture itself 
and devices used by various investigators have in- 
cluded pilot flames, sparks, detonators, overheated 
surfaces, and catalytic agents. The introduction of 
such stimuli is often accompanied by complex effects 
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her characteristics of the equipment, and in 
il it tends to complicate the problem of accurate 


yrretation and evaluation of the results Until 


extensive research has been carried out on 
ems related to then use it would appear 
sable to carefully avoid their introduction, 


intended or incidental 
3. Test Apparatus and Procedure 


3.1. Development of Apparatus 


liminary experiments were conducted with the 


iM standard ignition apparatus [21] with solde1 
and gas burner Because of difficulties in 
olling the temperatures and poor reproducibility 
e results, this apparatus was modified to include 
ic heating elements, a 0.2-liter spherical Pyrex 
a stirring arrangement in the bath, and direct 
surement of temperature in the flask. In order 
educe the heat loss from the reacting mixture to 
surroundings and improve uniformity of temper- 

re control in the flask, the solder bath was later 
laced with two sets of heating elements made of 
nichrome wire, 1n asbestos sleeving, wound around 
the upper and lower halves of the spherical flask 
and protected with a refractory cement. The final 
design of the ignition apparatus consisted of a 1-liter 
spherical Pyrex flask, as shown in figure 2, placed 
meentrically inside a refractory that a 
n. air space was provided between the flask and 
the casing for further reduction of heat loss. The 
heating wires were wound around the upper and 
lower halves of the casing and the neck, with inde- 
pendent temperature controls. ‘Temperatures of 
the three heaters were indicated by three thermo- 
made of chromel-alumel (26 B&S 
wire placed behind each group of heating 
Two additional thermocouples were placed in con- 
tact with the upper and lower surfaces of the flask. 
The sixth thermocouple, insulated by glass sleeving 
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Ficure 2 One liter spherical pyrex flask ignition apparatus 


for liquids with air space and refractory casing 


and also of 26 (B&S gage) wire, was placed in the 
lower part of the flask for indicating the tempera- 
ture inside the flask The flask neck was 
with a cement-asbestos stopper having two passages, 
one for the thermocouple wires and the other for 
injection of the charge 

The effect of chamber size on the observed self- 
ignition temperature was studied by means of addi- 
tional test units of the following capacities: 0.008, 
0.035, 0.2, 12, and 15 liters. ‘These units were con- 
structed in substantially the same manner as shown 
in figure 2, with the exception of the 15-liter appara- 
tus, the chamber walls of which were made of opaque 
in. thick, with no air space between 
the walls and the heating coils. From the 
shown in table 3 and figure 5 it is evident that, over 


closed 


silica about 
results 


; ; 


TABLE 3 Effect of flask size on ignition characteristics 
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self-ignition 
temperature becomes lower with increasing chamber 


the size range studied, the observed 


This 


rate of 


size in units of similar construction effect 
should not be unexpected as the heat 
loss per unit volume from the reacting mixture will 
be smaller with the larger chambers because of the 
smaller ratio of surface to volume. The _ higher 
values of SIT obtained with the 15-liter unit than 
with the 12-liter unit indicate a greater heat loss to 
the thicker opaque silica walls. Although slightly 


lower SIT values were obtained with the 12-liter 
Pyrex unit, it was unwieldy and fragile. Conse- 
quently, the standard, easily replaceable 1-liter 


flask was used as the more practical test apparatus, 
in making most of the determinations 

An investigation was also made of the comparative 
SIT values obtained with different types of casing 
around the 1-liter flask In addition to the refrac- 
tory casing separated from the flask by a ‘:-in. air 
space shown in figure 2, two other units were con- 
structed. In one of these a \-in.-thick copper casing 
was fitted as closely as possible around the flask, and 
in the other, a refractory casing with mirrored inner 
surface was separated from the flask by a ‘-in. air 
space The results obtained with several different 
liquids tested in these three units are shown in table 4 
In each case the highest value of SIT was obtained 
when using the metal-encased flask, whereas the 
lowest value was obtained with the mirrored refrac- 
tory casing. However, it appeared that the mirror 
coating on the refractory casing would not hold up 
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well at high temperatures and would be sub 
variable deterioration with repeated use. A . 
important that the thermal characteristics of t 
paratus remain as constant as possible, the 
refractory casing was selected for the major 

of the work 


TABLE 4 Effect of chamber wall construction on self 


femperature with the I-liter flask 
x 
‘ 
Liquid 

Meta Refractory 
r 
Acetic acid (1 i8N 
Acetome “MO 401 
Ethyl alcohol 40. 91 
CGiasoline (motor 249 243 
Isooctane 441 1") 
Noluene ”y 19 
Varsol 237 23 
Ethyl ether 172 “0 


The possibility that the thermocouple wires might 
exert) some catalytic effect on the reaction 
seriously considered. A test apparatus was carefull) 
prepared so that heat losses by conduction 
practically eliminated. | Comparative 
conducted, with the thermocouple remaining in its 
usual position in the chamber threughout th 
reaction, and with the thermocouple removed from 
the chamber just before injection of the liquid 
charge. The self-ignition temperature was a litth 
higher with the thermocoup!e present than without 
it, and increased slightly as the length of the wir 
in the chamber was increased, but no evidence of 
cataly tic action was detected 

Callendar [23] conducted experiments on th 
possible catalytic influence of a variety of metals 
placed in a test tube and, of those studied, onl) 
platinum black showed catalytic action. Thompson 
(24], and Matson and Dufour [25], used different 
metals for construction of the entire wall of the 
reaction chamber. Definite variation of the self- 
ignition temperature with the metal used was found 
but the effects were evidently related to the heat 
losses to and through the container walls 
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3.2. Test Procedure 


The procedure followed in establishing the self- 
ignition temperature of a combustible material con- 
sisted of a series of experiments at different tempera- 
tures and with different proportions of combustible 
material and air. In conducting an experiment, tly 
power supplied to the individual heating elements 
of the ignition apparatus was adjusted until the 
temperature inside the flask was uniform to 2 deg C 
and constant at some value estimated to be nea! 
the SIT. A charge of 0.2 ml of the combustible 
liquid per liter of flask volume was injected into the 
flask by means of a hypodermic syringe fitted with 
a fine needle. This provided the dispersion necessary 
for rapid evaporation and nearly simultaneous 
initiation of the reaction throughout the flask. T! 
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was permitted to proceed without any 
ons of the power supplied to the apparatus 
rder to follow the nature of the reaction, 
itions of the electromotive foree of the 
wouple inside the flask, at intervals of 15 to 
ynds, were plotted against time as the reaction 
ded. The temperatures of the heating coils 
isk surface were checked at intervals of 3 to 
ites 
experiment was repeated at different initial 
ratures until the lowest temperature from 
ignition occurred was determined for the 
of 0.2 ml of combustible material per liter of 
Additional experiments were then made with 
nt charges to establish the lowest temperature 
vhich the optimum ratio of combustible mate- 
»air would proceed to iwhnition 
Generally, ignition was clearly indicated by audible 
sible evidences (explosion, puff, flame), and it 
lso shown by characteristic changes in the time- 
erature curves, which served as dependable 
li when the visible evidence became uncertam 
»bseured As an aid to visual detection in the 
determinations, the interior of the reaction 
was viewed under dark room conditions by 
is of a omiurror placed above the flask opening 
er these conditions, even verv weak flames were 
netly visible, and flames obscured by the produc- 
on of dense smoke and fumes could be detected 
the illumination of the fumes 
\laterials such as gums, waxes, and tars, which 
too viscous at room temperature to be handled 
by syringe injection, were tested by placing them in 
small crucible suspended in the lower part of the 
flask by a thin wire This crucible arrangement was 
ilso used for at least the preliminary tests on certain 
oils, laequers, and paints, which produced especially 
sooty residues, because it helped to prevent an 
excessive accumulation of soot on the walls of the 
flask. It was necessary, however, to make a correc- 
tion for the cooling effect of the crucible, on any 
SIT values obtained in this way 
\fter each test, the products of decomposition 
were blown out of the flask by a jet of air in order to 
reduce the amount of carbon or soot deposited on 
the walls. The interior flask surface was cleaned at 
regular intervals during the preliminary tests for 
approximate determinations of the SIT, and before 
each test in the final determinations 


4. Test Results 


lhe following ignition characteristics were deter- 
mined under the specified test conditions for a 
imber of combustible liquids 
|. Self-ignition temperature 

Character and intensity of the evidence of 
lution 
3. Favorable fuel-air mixture ratio for minimum 
nition temperature. 
+. Time-temperature curves and ignition lags for 
ferent initial temperatures 
Values for the SIT given in tables 3 to 9 are the 


minimum values obtained under the most favorable 
conditions for ignition of the particular material 
Representative characteristics and time-temperature 
curves for the materials under study are shown in 
figures 3 to 22, but they do not in all cases include 
the determination that gave the lowest SIT value 
Some minor discrepancies will therefore be observed 
between the SIT values in the tables and those shown 
by the curves 

For convenience, the combustible materials studied 
in this investigation have been grouped as follows 

1. Liquids of low ignition temperature 

2. Common petroleum products 

3. Petroleum compounds 

1. Reference fuels 

5. Miscellaneous compounds 

6. Miscellaneous commercial materials 


4.1. Liquids of Low Ignition Temperature 


The ignition characteristics of two liquids, carbon 
disulfide and ethyl ether, ignitible at comparatively 


low temperatures are given in table 3 
a. Carbon Disulfide 


The self-ignition of carbon disulfide was usually 
indicated by a violent explosion accompanied by 
flame and a very rapid rise in temperature (see fig. 3 
The preignition reaction, however, was of a very 
mild form, indicated by a gradual rise in tempera- 
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Ficure 4 Effect of charge on wanition temperature of carbon 
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1-liter flask 











ture, even for temperatures just below the ignition 
Minor fluctuations in the time-temperature 
curves sometimes resulted from secondary reactions 
caused by rich mixtures or by from 
deposits on the flask surface if not cleaned before 
each test Both of these factors may have a retard- 
ing effect on the ignition reaction. Normally the 
lowest self-ignition temperatures for combustible 
liquids were observed when using charges of 0.2 to 
0.25 ml per liter of ai However, self-ignition 
temperatures of carbon disulfide showed a rapid drop 
as the fuel charge was from 0.2 to 0.05 
mi/liter of air and remained constant for 
charges below 0.03 ml, as shown in figure 4 

The carbon disulfide-air mixture had one eritical 
ratio for each initial temperature, above which the 
mixture remained nonignitible. For example, the 
critical charge for the initial temperature of 150° C 


region 


USCS « volved 


reduced 
almost 


was about 0.18 ml. With a charge of 0.25 ml, the 
reaction was not accelerated to ignition during 
the 2%-min test period If a fresh charge of air 


was admitted at the end of this nonignition reaction, 
This ignition probably 
concentration of the 
0.18 ml per 


ignition occurred 
from reducing the 
below the critical 


a violent 
resulted 
mixture 
liter 


value of 


b. Ethyl Ether 


Self-ignition of ethyl ether was characteristically 
indicated by a violent explosion accompanied by a 
large flame and rapid rise in temperature. Typical 
time-temperature curves are shown in figure 5 for 
different sizes of ignition chambers 


The nonignition reaction of this liquid was 
characterized by a slow rise of temperature for a 
considerable time after injection of the charge. 


During the first period of the exothermic reaction, 
extending from 5 to 10 seconds, the pressure inside 
the flask increased slightly and gases were expelled 
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from the neck with a barely audible sound 
temperature of the mixture continued to ris 
min, reached its maximu 
deg C the initial temperature Ho 
at slightly higher initial temperatures ig 
occurred at a time when the thermocouple ind 
a lower temperature of the mixture than the 
mum observed in the preceding test. This 
seemed to indicate that either the developm: if 
ignition is not a simple function of gas temper 
or the thermocouple and temperature-indi: 
system could not accurately follow the rapid 
perature changes at the point within the gas 


in about 10 or 15 
above 


ignition actually took place 


4.2. Common Petroleum Products 

The ignition characteristics of 18 common petro 
leum products were determined, and the final SIT 
values obtained are given in table 5 


PARLI > Tanitior characteristics of common pe 
nlucts determined in the 1-l flash 
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a. Motor Gasoline 


The time-temperature curves for motor gasolines 
were characteristic of a number of refined petroleum 
products and differed from those of combustible 
liquids of nonpetroleum origin. The exothermi 
reaction appeared to proceed in two stages. It was 
characterized by a consistently reproducible tem- 
porary drop in temperature as the reaction pro- 
gressed, probably caused by an interruption of the 
exothermic reaction or the occurrence of some endo- 
thermic The time-temperature curves of 
figure 6 show this two-stage reaction, the charac- 
teristics of which were not appreciably altered by 
varying the charge or the size of the flask. Ignition 


pre TCOCSS,. 


of gasoline usually occurred during the last stage of 
the reaction, whereas that of kerosine and cleaner’s 
solvent occurred in the first stage. 

A comparison of adjacent ignition and nonignition 
curves suggests that the rate of formation of stable 
compounds was highly sensitive to slight variations 
in temperature. 


Not only did the temperature in 
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tI nonignition reaction rise well above the SIT, 


but the reaction mixture had become so stable that 

vas not ignitible either by spark or flame, or by 
additional heating up to 400° C, unless diluted with 
fresh air 


lhe intensity of the evidence of ignition varied 
with the initial temperature. From initial tempera- 
tures close to the SIT, gasolines of high octane 
number tended to ignite with a dull blue flame 
observable under dark-room conditions, and gaso- 
lines of low octane number generally ignited with a 
mild orange flame accompanied by audible expulsion 
of gases (“puff”). Ignition at higher initial tem- 
peratures was usually indicated by a more violent 
flash of flame or by audible expulsion of vases 
It was noted also that ignition from high initial 
temperatures left the surface of the flask more 
heavily coated with soot and combustion products 
than did ignition from lower initial temperatures 

lenition lag was found to be closely correlated 
with the ignition temperature, as shown in figure 
7. The lowest self-ignition temperature was usually 
determined from tests conducted with different 
charges, as shown in figure 8 


b. Other Gasolines and Varsol 


The ignition characteristics of other gasolines and 
Varsol were substantially the same as those of motor 
and aviation gasolines (supply of 1937), with a 
two-stage reaction and with ignition occurring in 
the second stage of reaction 
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c. Kerosine and Cleaner’s Solvent 


Kerosine and cleaner’s solvent belong to another 
group of petroleum products, the ignition of which 
seemed to occur during the first stage of reaction in 
spite of the fact that the nonignition reaction pro- 






































aw Ficure 10 T'ime-temperature curves for cleaner's sol 
ceeded with a two-stage reaction ypu al time- F : 


temperature curves are shown in figures 9 and 10 
‘ After-burning”’ of cleaner’s solvent in the region of 
high ignition temperatures was more pronounced e. Diesel Fuels 
than with other fuels, and ignition at all tempera- 
tures was indicated by mild flame and a_ barely 
audible expulsion of CASES 


flask, charge equa 


The physical properties of four Diesel fuels, and 
their ignition characteristics determined in the 
liter flask, are given in table 6 


d. Crude Oil The time-temperature curves of figure 11 are ' 
representative for these fuels. The two upper | 
Lagunillas crude oil of gravity 17 API, flashpoint | curves were obtained for a straight-run fuel (55 — | 


15° ©, and sulfur content 1.91 percent, was tested | cetane); the two lower curves were obtained with 
in connection with a study of the explosion of an oil | the same fuel, to which 1 percent of ethyl nitrat 


tanker. ‘The SIT of this oil was found to be 262° C, | had been added. A small addition of the accelera 
with an ignition lag of 35 seconds. The ignition | tor reduced the ignition temperature from 230° to 


and nonignition reactions proceeded in a one-stage | 226° C, and the ignition lag from 2% min to 1% min 
pattern, in contrast with the two-stage reactions of | and the intensity of ignition changed from explosion 
the gasoline and kerosine groups to normal flame 





Tarnue 6 Physical properties and gnition characteristics of Diesel fuels 









* Key for evidence of ignition Mild flame; 3, flame; 4, mild explosion; 5, explosion 





58 











ppel 
»? 
with ' 
rate 
lera : 
A | 
min 


sion 














TEMPERATURE 





f. Furnace Oil No. 2 


lhe general ignition characteristics of furnace oil 
No. 2 were substantially the same as those of Diesel 
fuels of 41 and 55 cetane number, i. e., single-stage 
eaction with medium self-heating before ignition, 
vhich occurred at 233° C, with mild flame after 
> min 40 sec. 


g. Lubricating Oils 


Lubricating oils generally represent a complex 
mixture of different hydrocarbon components hav- 
ing different affinity for oxygen. Because of this, 
the self-ignition of oil vapor does not always proceed 
n a definite order, or as consistently as in the case 
of other combustible liquids. Oils having a com- 
paratively low ignition temperature, such as U.S 
Navy Nos. 2075 and 8190, ignited in a more or less 
efinite and reproducible manner. However, chlori- 
nated paraffin oils and motor lubricating oils used 
for internal-combustion engines usually showed 
rregular ignition in the zone of minimum ignition 
temperatures. Above this zone, ignition was ob- 
served regularly, and below it the oil vapors were 
definitely not ignitible. For chlorinated paraffin 
oil this doubtful zone covered a range of about 
0 deg C, but for motor oil it amounted to about 
) deg C, as indicated in figure 12. 
rhe reasons for this irregularity of ignition were 
established. However, from observation of the 
reacting mixture under a_ strong light directed 
rough the pyrex cover over the neck of the flask, 
appeared that poor distribution of the charge may 
ive been largely responsible. As these oils were 
aved into the reaction chamber, they formed 
use clouds that required an appreciable time to 
ome well mixed with the air throughout the 
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Figure 12 T'ime-temperature curves fo 


chamber. In this nonuniform mixture there was a 
tendency for ignition to take place in small local 
areas. At relatively high initial temperatures, such 
local ignition was commonly followed by multiple 
ignitions in other parts of the chamber, but at lower 
initial temperatures, a partial ignition frequently 
occurred without sufficient violence to affect the 
remainder of the mixture. Moreover, the 
clouds of unburned oil tended to obscure the evidence 
of a mild ignition, when viewed under dark 
room conditions, and detection was sometimes un 
certain. These variable conditions could be ex 
pected to produce irregular ignition results in the 
temperature range near the SIT. It was noted, 
furthermore, that the irregularities more 
pronounced as richer reaction mixtures or larger 
flasks were used, a situation particularly susceptible 
to nonuniform distribution of the charge 


dense 


even 


became 


h. Benzene 


The time-temperature curves obtained for pure 
benzene (C,H,) are shown in figure 13. Appreciable 
self-heating of the reacting mixture was observed 
50 deg C below the self-ignition temperature. At 
the initial temperature of 557° C, the mixture self- 
heated to 690° ©. However, at the initial tem- 
perature of 565° C, actual ignition was observed 
at 603° C, or about 87 deg C below the maximum 
temperature of the preceding test. Here, again, 
evidence 1s presented of the fact that the initial 
temperature from which the exothermic reaction 
proceeds to ignition is more significant than the 
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® The time-temperature curve the above liquids are ¢ 
hue flat 


temperature of the reacting mixture at the 
of whitton It is interesting to note that thy 
mum temperatures for each reaction fall on a 
Curve The shape of the nonignition reaction 
temperature curve is simular to but not 
with, that of singel stage reactions of other rhs 
Although both the tenition and nonmignition rv 
were characterized by rapid increases in te nper 
inition was indicated by a barely audible exp 
of dense smoke and mild flame 


1 Te iuene 


I ha elf-heating and inition characteri 
toluene were intilar to those of benzene. an 
temperature curves of substantially the 


were obtamed 
4.3. Petroleum Compound 


Mh groups of hydrocarbon compounds d 
from petroleum, were supplied by the divisio 
thermochemistry Phe irhhition characterist 
these hydrocarbons were determined im the 
apparatus, and are given in table 7 


a. Normal Paraffin 


Kigure 14 gives representative time-temper 
curves for fractions of the chain type, which sho 
substantially the same behavior, with the except 
of dodecane and cetane the noniwnition reaction 


which had a retarding portion similar to a_ ty 


stave reaction The rate of the mduc ine react 
preceding ignition gradually decreased, and 
whition lag increased with the carbon content 


iven in figures 
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elf-ignition temperature of liquids of the 
| 


type consistently diminished with an increasing 
ber of carbon atoms up to 16 C,H), as shown 
ie 15 The curve for fractions of larger 
bon content was not determined because of 
quantities of rare fractions. However, tests 


ery small amounts of dotriacontane (Cy 1H 

ated the self-ignition temperature to be some 
re above 300° C In other words, it would 
aur that the curve must turn upward after 
mon 16 Additional tests will be necessary to 


nfirm the shape of this curve 


b. Cycloparaffins 


The time temperature curves for evelohexane and 
hvlevelohexane, the two ring type fractions 
d, were similar. Curves for methvlevelohexane 
hown in figure 16 It will be noted that the 
tion appeared to be of a two stage type, with an 
nded period of constant temperature before the 
a heating of the second stage began lenition 

place during the second stage, and, as the 
al temperature was raised, the delay between 
es was gradually eliminated. The two-stage 
ern of reaction for evelohexane is nearly identical 
i the patterns obtamed for motor and aviation 
lines in the same region of initial temperatures 




















c. Branched-Chain Paraffins 


the minimum 
compound was 
evidences of 


Trime thylbutane 
self-ignition temperature of 
approached, the visible and 
ignition became weak, and precise determination of 
the SIT was difficult. Furthermore, the time- 
temperature curves for ignition and nonignition re- 
actions had substantially the same shape. However, 
initial temperature of 421° C the reaction 
definitely proceeded to ignition with a blue-violet 
flame. At lower light produced by 
the reacting mixture appeared to be that of a cool 


(7) iptane As 
r 6this 
audible 


at an 


temperatures 


flame reaction 

Tsooctane (2,2,4-T) imethyl pe ntane The 
characteristics of isooctane (and difficulties of ignition- 
temperature determination) were similar to those for 
trimethvlbutane. The observed reproducibility of 
blue flame, or cool flame, was very poor and was 
affected by deposits on the flask wall, by direction of 
the charge stream, and by relative humidity of the 
Injection of the charge along the side of the 


iwnition 


air 

neck wall usually stimulated ignition at slightly 
lower temperatures than injection directly to the 
bottom of the flask Self-ignition temperatures 


obtained at different times ranged from 408° to 
125° ©. The value of 420° C (with corresponding 
ignition lag of | min 20 sec) recorded in table 7 was 
fairly reproducible 

Although these difficulties were encountered with 
all sizes of flasks used, reproducibility of the deter- 
minations was slightly better in the 0.2 and 12-liter 
flasks than in the I-liter flask. The band of un- 
certainty in the 0.2-liter flask was no more than 
12 deg C, and in the 12-liter flask it was reduced to 
} deg ( ; 


4.4. Reference Fuels 

The reference fuels commonly used in determining 
the knock ratings of motor fuels comprise variously 
proportioned mixtures of isooctane and n-heptane 
Essentially the rating procedure consists in deter- 
mining what mixture most closely matches the per- 
formance of the fuel under study, using a standard 
CFR Cooperative Fuel Research C‘ommittee engine 
operated under specified conditions The percentage 
of isooctane in this mixture is taken as the octane 
number of the fuel being tested Details of the 
procedure differ with the purpose of the fuel, and the 
ASTM has established four different methods (motor 
research, aviation, and supercharged) for gasolines, 
and one for Diesel fuels The rating of a particular 
fuel may vary appreciably with the method used, 
and it is therefore necessary to consider the method 
of determination when knock ratings are used as a 
basis of comparison 

The ignition characteristics of these primary refer- 
ence-fuel mixtures were studied in the regular 1-liter 
test apparatus, and time-temperature curves for a 
number of them are shown in figure 17. In addition 
two other types of mixtures, formerly used as 
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secondary reference fuels were also tested 
secondary reference fuels were prepared frot 
specific fuels, and consisted of (1) mixtures 
A, with fuel Cy, and (2) mixtures of fuel ( 
fuel F Although improved methods of m 
ture and reduced cost of the primary referen: Is 
have now eliminated the need for intern 
standards, investigation of the well-defined se« 
reference fuels were of interest at the time th ate 
were performed 

In figures 18 and 19 the data obtained for pr 
and secondary reference-fuel mixtures, respect 


have been plotted to show the effect of char , 
composition on different ignition character Ss 
With increasing isooctane content, up to aboy 
70 percent, the SIT of the primary reference-fy 


mixtures increased gradually and continuously, an 
the ignition was marked by the characteristic pe 
flame of n-heptane. The intensity of ignition changed 
gradually from explosion to mild flame as the jso 
octane content increased With isooctane concer 
trations of 50 percent and higher an audible expu! 
sion of gases without flame was observed at tempers 
tures below ignition, and this ‘“‘puff” also diminishe: 
in intensity as the concentration of isooctane in 
creased, At 70 percent tsooctane, the SIT-compos 
tion curve showed a break; mixtures containing mor 
than 70 percent of isooctane igniting only at decid 
edly higher temperatures, and generally with the blue 
flame and considerable variability noted previously) 
for pure isooctane. Ignition was somewhat errati 
also in mixtures of slightly less than 70 percent iso 
octane content, forming a so-called 
nonignition zone’? similar to 
Prettre [15] with pentane and hexane by Townen 
[16] with hexane, by Sortman and Beatty [5] wit 
heptane and cetane, and by many other 
gators. The fact that nonignition zones with thes: 
latter liquids were not observed in the present 
investigation suggests that the appearance of suc! 
zones may be a function of the test apparatus an 
conditions rather than of the material under study) 
The series of secondary reference-fuel mixtures 
prepared from fuels Ay and C,, showed a continuous 
linear increase in ignition temperature with increasing 
C,. content as indicated in curve II of figure 19 
The series of mixtures of fuels Cy. and F;, however 
gave a pattern of ignition behavior (curve III 
similar to that obtained for the primary reference- 
fuel mixtures, including a secondary nonignition 
zone and the “puff without flame at temperatures 
below ignition for mixtures rich in the component 
having the higher ignition temperature 
In order to indicate the degree of 
between self-ignition temperature and motor-method 
octane number, SIT 
plotted for both the primary and secondary refer 
ence-fuel mixtures. In addition, self-ignition tem 
perature determinations were made on three othet 
gasolines of known octane number. 
ary reference fuels, only the end fuels A, and F 


“secondat 


those observ ed b 


invest 


correlation 


showed substantially the same SIT values as the 
reference fuels of corresponding octane 


primary 














octane-number curves were 


Of the second- 
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refe rence fuels 


1-liter flask 








temperature curves of fuel C, 


1-liter flask, charge equals 0.2 ml 











number ~ Onlv one of the three additional cvasolines 
tested gave an SIT value approximately the same 
as that of the primary reference fuel of corresponding 
octane number Similar poor correlation between 
ignition characteristics and octane-number ratings 
have been noted by other investigators with other 
test apparatus Jentzsch Eckert Zerbe and Carhart 

The effects of a small amount of ethyl fluid addi 
tive on the ignition characteristics of a fuel were also 
studied A typical time-temperature curve for the 
ignition of fuel Cy. (78.9 octane number) is shown in 
the lower part of figure 20.) The broken-line curves in 
the upper part of the figure are time-temperature 
curves for nonignition reactions, at increasing initial 
temperatures, of the same fuel with 1.25 ml of tetra 
ethyl lead added per gallon. The shape of these 
curves changed entirely, and manv additional stages 
of reaction developed before the reaction finally 
accelerated to iwhition 


4.5. Miscellaneous Compounds 


the 1-liter flask for 
XS 


The SIT values obtained i 
seven miscellaneous liquids are given in table 
The effect of varving the flask size is shown for several 
of them in table 3 

For all except methyl alcohol, both the ignition and 
nonignition reactions were characterized by an im- 
mediate rapid acceleration, as illustrated by the 
curves for acetone shown in figure 21. With methyl 


aleohol (fig. 22), the acceleration was somewhat 
delaved, and the curves for ignition and nonignition 
reactions differed distinctly in shape. Additional 


comments on the reaction characteristics of the 
individual liquids follow 
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a. Acetone 





At the minimum self-ignition temperature, acetone 
ignited with a blue-violet light, often initiated at and 
spreading from a small, red, glowing center of carbon- 
ized deposit of acetone on the flask wall. The light 
usually appeared without audible manifestation and 
continued to be visible for about 20 seconds at 491° C 
and about 15 seconds at 505° C The nature of 
this light gradually changed from dull blue-violet at 
191° © to blue-pink at 515° C and to a pink-violet 
flash with an audible explosion at higher tempera- 
tures. The initial blue-violet light was similar to 
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Fiaeure 21 Time-temperature curves for acetone 
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Figure 22 Time-temperature curves for methyl ak 









1-liter flask, charge equal 21 






that of a cool flame, and when first observed it was 
questioned whether this could be taken as an indica 
tion of ignition. However, in consideration of the 
very sharp increase in the time-temperature curve 
fig. 21) at 491° C, this temperature was taken as 
the self-ignition temperature 


















b. Acetic Acid 






tion of acetic acid occurred so irregularly that 
difficult to define the lowest self-ignition tem 

One possible explanation of such irregu- 
nav lie in the ability of acetie acid to exist in 
erchangeable states, “dimer and ‘“‘monomer”’, 
sed by the relationship 





CH,COOH MCH,COOH 
o forms may have different ignition tempera 
ind uncontrolled variables might have caused 
from one form to the other 
tion of acetic acid, in spite of the irregularity 
irrenee, Was distinetly defined by audible and 
evidence, and usually occurred about 8 seconds 
njection of the charge. Ignition was initiated 
blue flame, which, at the end of the ignition 
changed to a pink color. This pink color 
ho vel might have been due to reflection of the 
d furnace walls from fog formed within the 
bet Although the mixture remained clear and 
sparent throughout most of the reaction, there 
some evidence of the development of fog toward 
nd of the exothermic reaction 


c. Ethyl Alcohol 


Due to the rapid increase of pressure at the time 
evaporation, an audible expulsion of gases was 
before ignition. The actual ignition, at 
1° C after 40 seconds, was accompanied by blue 
flame with a second expulsion of vases 
Ethvl aleohol-water mixtures were also found to 
readily ignitible. The self-ignition temperature 
reased from 391° to 444° C as the water content 
if the mixture approached 70 percent. Very weak 
even with 80 percent 


obse rved 





lame was observed at 467° C 


d. Ethylene Glycol 


Che ignition of ethylene glycol in the 1-liter flask 
was indicated at 398° C by orange-blue flame and a 
rily audible expulsion of gases. However, the 
same liquid containing 3 percent of an antirusting 
ngredient (triethanolamine phosphate) ignited at 
5° C with a violent explosion. ‘Tests on ethylene 
© viveol-water solutions showed ignition of mixtures 
ntaining up to 90 percent of water (for which an 
enition temperature of 460° C was obtained), 
tl percent of water the 








hough with more than 40 
ne became very weak. 


wo other antifreeze liquids, R-19 and R-19d, 

it was \ tested, and the results indicated that these 
ndica- ds were highly ignitible. Both liquids were 
of the troleum distillates (kerosine type), liquid R-—19 
curve r the product without additive and R—-19d being 
ken as xture of the product with an ingredient reported 
prevent the softening of rubber hose in engine 

nstallations. The self-ignition temperature of R-19 

Was 211° C, with an ignition lag of 3 min 50 sec; 


R-19d ignited at 223° C after 1 min 30 sec. Ignition 








of both liquids occurred suddenly with a violent 
explosion when the thermocouple indicated a temper 
ature only 3 deg C above the initial temperature 


e. Methyl Formate 


Ignition of methyl formate indicated regu 
larly by a mild explosion and rapid rise in tempera- 
ture similar to that for acetone (fig. 21 


was 


{. Glycerine 


Glycerine ignited with a violent explosion accom 
panied by long flame and smoke expelled from the 
neck of the ignition apparatus. The rise in tempera- 
ture of the gas mixture immediately after injection 
of the charge was particularly extreme for both igni 
tion and nonignition reactions 

Observation of the flask interior 
immediately after injection of the charge, 
glycerine droplets accumulated on the bottom sur 
face of the flask and rapidly produced a dense white 
fog. After about 20 seconds, this fog became clear 
either due to superheating or partial oxidation 
This was followed, usually in about 5 seconds, by 
the appearance of a small blue flame with a pink 
center. Due to rapid propagation of flame, gases 
were expelled from the neck of the flask with con- 
siderable force 


that, 


Viscous 


showed 


g. Methyl Alcohol 


Methyl alcohol showed a mild form of ignition 
characterized by a dull blue flame and a weak expul- 
sion of clear gas, but the difference in shape of the 
time-temperature curves for ignition and nonigni- 
tion reactions (fig. 22) provided supplemental evi- 
dence of the definite location of the SIT 


4.6. Miscellaneous Commercial Materials 


A number of commercially produced combustible 
materials were tested in the 1-liter ignition apparatus 
and their ignition characteristics are given in table 9 
The time-temperature curves typical of many of 


these materials were similar to those for acetone 
fig. 21 
Parte 9 Ignition characteristics of miscellaneous commercial 
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a. Black Enamel 


Noticeable decomposition of synthetic black auto- 
mobile enamel (601) was indicated at 270° C by the 
appearance of smoke. At a temperature of about 
145° © carbonized residue of the charge self heated 
to the glow state. However, self-ignition occurred 
at 478° C with a large flame and audible expulsion of 
gas and smoke As a result of ignition, all the carbon 
residue was burned out and the flask surface became 
cleat 

b. Pine Tar 

Pine tar of commerical grade, due to its low 
‘fluiditv’’, was preheated in boiling water before in 
jection with a hypodermic syringe. The exothermic 
reaction began far below the ignition temperature 
and ignition indicated by flame and a mild 
expulsion of smoke. The flask after igni 
tion, was usually covered by a brownish coating and 
tar residue, baked to a hard carbon deposit removable 
only with a However, with pine oil the 
residue after ignition was soft and easily removable 
from the surface 


was 
surface, 


scrape 


c. Linseed Oil 


The ignition characteristics of linseed oil were 
similar to those of lubricating oils, and appreciable 
self-heating of the vapor-air mixture 


at temperatures far below the ignition temperature 


Wiis obser ved 


d. Varnish and Lacquer 


The reaction characteristics of varnish and lac- 
quers were substantially the same, particularly in 
the nonignition region. Self-heating of the vapor-air 
mixture initiated far below the self-ignition 
temperature, and often amounted to 20 to 25 deg C 
in nonignition reactions near the self-ignition tem 
perature, Ignition of these liquids, with the excep 
tion of interior varnish and lacquer thinner, was very 
mild; both the laequer thinner and interior varnish 
ignited with an audible explosion 


Wiis 


e. Turpentine 


The time-temperature curves for turpentine were 
similar to those for kerosine, having the same two- 
stage nonignition reaction with ignition definitely in 
the first stage 


5. Summary 


The apparatus and method described for the 
determination of self-ignition temperatures is prac- 
tical and readily adaptable to combustible liquids of 
widely different physical and chemical properties 
Data on a considerable number of liquids are pre- 
sented 

Time-temperature curves of a thermocouple junc- 
tion in mixtures of combustible liquid and air, offer 
a convenient means of studying the preignition 
reactions in such mixtures. Typical time-tempera- 
ture curves are given, 


66 


The importance of careful standardizati 
equipment and procedure in obtaining reprod 
results has been demonstrated. Thus, for rey 
vessels of similar construction, the self-ignitioy 
perature was found to decrease with increase jy 
of the chamber, over the capacity range st) 
with differences as high as 200 deg C being obt 
in vessels of 0.008-and 12-liter capacity Wit! 
rex reaction chambers of the Same size but 
different exterior casings, Variations up to near 
deg C' were observed in the determined self-ign 
temperature In view of the effects produces 
these modifications of equipment within a ¢ 
method, the large discrepancies in ignition tem 
tures reported by different investigators 
various methods are not surprising 

The temperatures measured with 
apparatus described are lower than those repo 
by most other workers 


ivnition 


The author is indebted to \l Sandholzer for; 
ance in the preparation of this paper 
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